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Preface 


This  research  effort  concentrates  on  the  develop-  , 
ment  of  flight  control  laws  for  the  X-29  forward  swept 
wing  demonstrator  aircraft  by  applying  recent  advances  In 
multivariable  control  law  theory.  Such  an  approach  Is  of 
great  Importance  to  the  Air  Force  because  past  design  tech¬ 
niques  are  very  tedious,  time  consvunlng,  complex,  and 
expensive  when  applied  to  the  design  of  the  Intercoupled 
multiple  Input,  multiple  output  flight  control  systems  now 
required  on  new  advanced  technology  aircraft. 

The  author  wishes  to  express  his  gratitude  to 
Captain  David  Potts  and  others  at  the  Plight  Dynamics 
Laboratory  for  sponsoring  this  thesis. 

The  technical  expertise,  assistance,  suggestions, 
and  stimulating  technical  discussions  provided  by  my  thesis 
advisor  Professor  John  J.  D'Azzo  had  a  profound  effect  on 
the  work.  Also  appreciation  Is  expressed  to  Professor 
D'Azzo  and  my  reader,  Professor  Constantine  H.  Houpls,  for 
carefully  and  thoroughly  reviewing  the  manuscript.  A  very 
special  thanks  to  Captain  John  West  at  NASA  Dryden,  Edwards 
AFB,  California  for  his  timely  and  Important  contributions 
in  obtaining  appropriate  models  of  the  aircraft.  Most 
importantly,  I  am  deeply  grateful  to  my  wife,  Patricia , 
whose  assistance,  patience,  understanding,  and  constant 
supportiveness  made  the  success  of  this  effort  possible. 

Major  Terry  L.  Courtheyn 
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AFIT/GE/ENG/84D-21 


Abstract 

Flight  control  laws  are  designed  for  the  X’’29  for¬ 
ward  swept  wing  demonstrator  aircraft  using  a  design  tech¬ 
nique  based  on  multivariable  control  law  theory  developed 
by  Professor  Brian  Porter  of  the  University  of  Salford, 
England.  The  computer-aided  design  program  called  MULTI 
is  used  to  develop  and  refine  the  control  laws.  MULTI  also 
simulates  the  complete  closed-loop  control  system  and 
generates  appropriate  time  response  plots  for  analysis. 

Aircraft  dynamics  for  several  points  in  the  flight 
envelope  are  represented  by  linearized  state  space  equations 
obtained  from  NASA  Dryden,  an  agency  responsible  for  the 
development  and  testing  of  the  X-29.  Decoupled  longi¬ 
tudinal  and  lateral  equations  are  used  to  design  separate 
longitudinal  and  lateral  controllers. 

Control  laws  are  developed  to  stabilize  the  aircraft 
and  perform  lonaitud in?*!  inun?;i'vers  (direct  climb,  vertical 
translation,  and  beta  pointing)  at  three  different  flight 
conditions,  with  and  without  first  order-actuator  dynamics 
and  computational  time  delay  added  to  the  simulation.  The 
responses  are  compared.  Although  some  degradation  of  per¬ 
formance  is  observed  when  the  actuator  dynamics  and  time 
delay  are  added,  the  results  show  fast  and  well-behaved 
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control  of  the  aircraft  motion  parameters  with  quick  and 
smooth  tracking  of  pilot  input  commands.  Similarly/ 
lateral  control  laws  are  developed  for  the  coordinated 
turn  and  beta  pointing  maneuvers  at  five  different  flight 
conditions  with  comparable  results. 

Finally/  future  areas  of  research  are  recommended 
which  also  include  proposed  modifications  and  additions  to 


MULTI. 


MULTIVARIABLE  CONTROL  LAW  DESIGN 


FOR  THE  X-29 


I .  Introduction 


l-l  Background 

Fighter  pilots  will  not  be  able  to  accomplish  their 
missions  unless  future  tactical  aircraft  incorporate  the 
latest  in  technological  advances  to  counter  the  ever- 
increasing  lethal  capability  of  the  enemy's  air  defense 
systems.  Adversaries  are  dramatically  improving  their 
surface  to  air  missiles,  radar  guided  anti-aircraft  artil¬ 
lery,  fighter-interceptor  aircraft,  air  to  air  missiles, 
and  other  air  defense  weapons. 

In  order  to  neutralize  these  deadly  threats,  tac¬ 
tical  aircraft  of  the  future  need  the  capability  to  fly 
higher,  farther,  faster,  and  to  fly  with  greater  maneuver¬ 
ability  and  survivability  at  all  altitudes  and  airspeeds. 

To  accomplish  these  goals,  future  vehicles  must  be  smaller, 
lighter,  more  fuel  efficJent,  and  have  better  aeronautical 
designs . 

New  airplanes  like  the  X-29  advanced  technology 
figh+  ^;r  ire  expected  tc  »,\eet  these,  challenges  with  innova¬ 
tive  .f.c  ..lU'M!;  such  as  forwa.T.d  swept  wings,  static  insta¬ 
bility,  ai.''l  hiultiple  functioti  ->1  surfaces. 


However,  the  aerodynamic  designs  which  provide 
this  enhanced  capability  cause  two  serious  problems. 

Because  of  the  statically  unstable  airframe,  It  Is  Impos¬ 
sible  for  a  pilot  to  fly  the  aircraft  without  automatic 
stability  augmentation  provided  by  a  flight  control  com¬ 
puter.  In  addition.  It  la  not  feasible  for  a  pilot  to 
independently  control  all  of  the  added  flight  control 
surfaces.  But  a  properly  designed  flight  control  system 
can  accomplish  this  difficult  task  by  optimally  deflecting 
the  appropriate  control  surfaces  when  the  pilot  commands 
a  specific  maneuver.  Thus,  the  performance  requirements 
for  future  fighter  aircraft  dictate  the  need  for  sophisti¬ 
cated  flight  control  systems. 

Digital  flight  control  systems  will  be  standard 
equipment  on  most  future  high  performance  aircraft.  Rapidly 
advancing  technology  In  digital  electronics  and  computa¬ 
tional  proficiency  of  microprocessors  make  them  Ideal  for 
solving  the  complicated  algorithms  required  in  state-of- 
':he-art  flight  control  systems. 

The  question  arises  as  to  the  best  method  for 
designing  the  control  laws  and  algorithms  for  the  flight 
control  digital  computer.  Classical  control  theory  tech¬ 
niques  are  beat  suited  for  the  simple  single  Input,  single 
output  (SISO)  flight  control  systems  currently  used  on  con¬ 
ventional  aircraft.  These  older  methods  are  very  tedious 
and  time  consuming  for  designing  complex  intercoupled 
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multiple  Input,  multiple  output  (MIMO)  systems.  In 
response  to  this  problem,  experts'  In  flight  control  theory 
have  developed  several  new  MIMO  design  methods  such  as 
linear  quadratic  regulator  (LQR)  and  quantitative  feedback 
theory  (QFT) .  This  thesis  uses  another  technique  recently 

* 

developed  by  Professor  Brian  Porter  of  the  University  of 
Salford,  England  (10) . 

1-2  Problem 

The  static  instability  of  the  X-29  plus  the  added 
number  of  flight  control  surfaces  dictate  the  necessity 
for  a  full  authority  digital  flight  control  system. 

Engineers  are  having  difficulty  designing  optimum 

t 

control  laws  for  the  X>:29  digital  flight  controller  because 
of  the  highly  relaxed  static  stability  and  the  complexity 
of  the  aircraft's  aerodynamic  design. 

Using  the  Porter  technique  and  the  computer-aided 
design  package  called  MULTI,  this  thesis  designs  improved 
control  laws  for  the  longitudinal  and  lateral  flight  con¬ 
trollers  for  the  maneuvers  and  flight  conditions  listed 
below: 

Longitudinal  Maneuvers 

1.  Dij:eci.  climb 

2.  Vertical  i-radijlation 

3.  Pitch  pointing 


Flight  Conditions 


1.  0.4  Mach  at  sea  level 

2.  0.7  Mach  at  15,000  feet 

3.  1.2  Mach  at  15,000  feet 

Lateral  Maneuvers 

1 .  Coordinated  turn 

2.  Beta  pointing 

Flight  Conditions 


1. 

0.4 

Mach 

at 

sea  level 

2. 

0.9 

Mach 

at 

sea  level 

3. 

0.7 

Mach 

at 

15,000  feet 

4. 

1.2 

Mach 

at 

15,000  feet 

5. 

0.9 

Mach 

at 

50,000  feet 

1-3  Scope 

This  project  expands  on  the  work  completed  by 
Lieutenant  Scott  Feldman  on  a  digital  flight  controller 
for  the  X-29  (5) .  Those  designs  did  not  include  actuator 
dynamics  and  computational  time  delay  which  reduce  the 
phase  margin  and  make  it  a  harder  task  to  assure  stabiliza' 
tion  of  the  system.  Specifically,  this  research  effort 
accomplishes  the  following  goals: 

1.  Redesign  the  flight  control  system  with 
actuator  dynamics  included  in  the  closed-loop  simulation 
but  external  to  the  plant  state  equations. 


2.  Evaluate  the  aircraft  responses  with  and 


without  computational  time  delay  and  actuator  dynamics 
included  in  the  closed  loop  simulation. 

3.  Enhance  the  capability  of  the  computer  program 
MULTI  by  developing  additional  FORTRAN  source  code  to: 

a.  calculate  the  steady- state  G(0)  and  0(0)“^ 
matrices  for  use  in  calculating  steady  state  aircraft  con¬ 
trol  deflections;  and 

b.  convert  output  simulation  data  from  dimen¬ 
sions  of  radians  to  degrees  prior  to  plotting  the  time 
responses. 

1-4  Standards 

This  research  effort  evaluates  the  final  digital 
control  system  design  using  computer  simulations.  The  con¬ 
troller  should  meet  the  following  standards: 

1.  Provide  stability  augmentation. 

2.  Show  fast  and  well-behaved  control  of  aircraft 
angles,  rates,  and  velocity. 

3.  Quickly  and  smoothly  track  the  pilot's  input 
commands. 

4.  Feedback  motion  parameters  that  are  measured 
with  currently  available  sensors. 

1-5  Approach 

This  research  work  applies  the  design  technique 
of  Professor  Brian  Porter  to  develop  the  appropriate 
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control  laws  for  the  X-29  digital  flight  controller.  His 
method  uses  the  state  space  time  domain  approach  of  multi- 
variable  control  theory.  Specifically,  integral  plus  pro¬ 
portional  control  is  applied  to  all  forward  loops.  In 
addition,  his  procedure  employs  output  feedback  which  has 
a  great  practical  advantage  because  the  outputs  to  be  con¬ 
trolled  are  accessible  and  easily  measured  with  sensors. 

More  information  is  presented  in  Appendix  A  where  the  theory 
of  Professor  Porter's  method  is  discussed  in  greater  detail. 

1-6  Materials  and  Equipment 

This  research  effort  uses  the  Control  Data  Corpora¬ 
tion  "Cyber"  main  frame  computer  located  at  Aeronautical 
Systems  Division,  Wright-Patterson  AFB,  Ohio  to  run  the 
software  programs  called  TOTAL,  ZERO,  and  MDLTI.  The 
computer-aided  design  package  called  MULTI  is  needed  to 
apply  Porter's  control  law  theories.  The  program  also 
allows  on-line  design  changes  and  generates  output  response 
plots.  A  comprehensive  user's  manual  with  illustrative 
examples,  describes  the  available  options  and  capabilities 
of  the  program  (8) . 

1-7  Overview 

Chapter  2  begins  with  a  general  description  of  the 
X-29  aircraft  and  continues  with  a  discussion  of  its 
unique  design  features.  Then,  Chapter  3  presents  the 
origin,  assumptions,  limitations,  and  dimensions  for  the 


linearized  state  space  equations  of  motion  for  the  X-29. 
Next,  Chapter  4  describes  the  detailed  procedure  developed 
in  this  research  effort  to  design  a  proportional  plus  inte¬ 
gral  multivariable  discrete  time  tracking  controller  for 
the  X-29  in  the  longitudinal  mode.  Detailed  results  for 
the  longitudinal  control  system  designs  are  presented  in 
Chapter  5  for  three  maneuvers  from  a  sample  of  flight  con¬ 
ditions.  Chapter  6  highlights  significant  steps  related 
to  the  lateral  design  procedure.  This  is  followed  by  a 
^presentation  and  analysis  of  the  results  of  the  lateral 
controller  designs  for  two  maneuvers  at  one.  flight  condi¬ 
tion.  Conclusions  are  given  in  Chapter  7  with  recommenda¬ 
tions  for  future  research.  Appendix  A  provides  a  descrip¬ 
tion  of  multivariable  control  theory  as  it  applies  to  the 
Porter  technique.  Then,  Appendices  B  and  C  document  improve< 
ments  made  to  the  computer  program  MULTI.  Next,  Appendices 
D,  E,  and  F  present  additional  data  not  included  in  the 
discussions  of  Chapter  5.  This  additional  data  includes 
state  space  matrices,  design  parameters,  controller 
matrices,  figures  of  merit,  and  response  plots  for  the 
longitudinal  controllers.  Similarly,  Appendices  G  and  H 
present  the  remainder  of  the  lateral  information  not  pro¬ 
vided  in  Chapter  6.  Finally,  Appendix  I  overviews  some 
unsuccessful  design  attempts. 
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2.  The  X-29  Aircraft 


2-1  General  Description 

The  X-29,  illustrated  in  Fig.  2-1,  is  an  experi¬ 
mental  advanced  technology  demonstrator  aircraft  built  to 
evaluate  the  feaslbixxty  and  capability  of  several  new 
aeronautical  design  concepts.  Specifically,  the  vehicle 
is  designed  to  exercise  the  full  spectrum  of  forward  swept 
wing  (FSW)  technologies  in  a  cost  effective  and  timely  man 
ner.  For  example,  56  percent  of  the  aircraft  weight 
includes  components  common  to  existing  aircraft  114:11). 
However,  the  remainder  of  the  vehicle  incorporates  many 
new  and  unique  design  features  as  outlined  below: 

1.  statically  unstable  airframe 

2.  forward  swept  wings 

3.  thin  supercritical  airfoils 

4.  close  coupled  canards 

5.  strake  flaps 

6.  variable  camber  control  surfaces 

7.  control  configured  design 

The  remainder  of  this  chapter  describes  these  new  innova¬ 
tions  in  more  detail. 


o 


Figure  2-1.  Grumman  x-29  Aircraft 

(16) 

Dimensions,  External: 

Wing  span 

27 

ft 

2 

in 

Wing  chords  at  root 

9 

ft 

8 

in 

at  tip 

3 

ft 

11 

in 

Wing  aspect  ratio 

4 

Length  of  fuselage 

48 

ft 

1 

in 

Height 

14 

ft 

3 

in 

Wheel  base 

18 

ft 

0 

in 

Weights 

Weight  empty 

13 

,326 

lb 

Max  fuel  weight 

3 

,978 

lb 

Max  takeoff  weight 

17 

,303 

lb 

2-2 


Airframe 


Grumman  designed  the  small  single  engine  aircraft 
with  an  overall  shape  developed  by  blending  an  optimized 
forward  swept  wing  and  canard  arrangement  to  an  exsltlng 
F-5A  forward  fuselage  module  with  modified  engine  ducts  and 
a  new  mid  and  aft  fuselage  (14i3}. 

Engineers  purposely  designed  the  airframe  with  high 
static  Instability  to  maximize  the  vehicle's  maneuverabil¬ 
ity  and  to  minimize  trim  drag.  As  a  result,  the  short 
period  roots  are  unstable  at  subsonic  speed  flight  condi¬ 
tions.  The  time  to  double  amplitude  Is  loss  than  0.5 
seconds.  This  rate  of  divergence  of  the  unstable  mode 
Illustrates  the  need  for  a  highly  reliable  stability  aug¬ 
mentation  system  (14:24). 

2-3  Forward  Swept  Wings 

Background.  Aeronautical  engineers  Invented  both 
forward  and  aft  swept  wings  during  World  War  II  in  a  concen¬ 
trated  effort  to  Increase  the  maximum  and  cruise  speeds  of 
aircraft.  Both  aft  and  forward  swept  configurations  delay 
Mach  drag  divergence  and  produce  weaker  oblique  wing  shock 
waves  at  transonic  velocities.  Initial  designs  had  wings 
swept  In  both  directions  and,  In  fact,  the  Germans  flew 
the  first  swept  wing  aircraft  (the  Jungers  JU  287  bomber) 
during  World  War  II.  It  had  forward  swept  wings  (6:147). 

Uni tuna te ly ,  continued  development  of  the  FSW  was  short 
■'.ived.  After  World  War  II,  engineers  attempted  additional 
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forward  swept  wing  designs,  but  no  aircraft  prototypes 
were  ever  built  due  to  an  inherent  malady — aeroelastic 
divergence.  The  existence  of  the  problem  was  known  from 
the  beginning  but  the  severity  was  not  fully  appreciated 
until  the  late  19408  when  Diederick  demonstrated  a  drastic 
drop  in  the  divergence  speed  with  increased  forward  sweep 
angle  (6:147).  The  divergence  phenomenon  increased  the 
FSW  weight  penalty  to  impractical  values.  Thus,  FSW  tech¬ 
nology  was  abandoned  and  designs  concentrated  solely  on 
aft  swept  wings.  Thirty  years  later,  the  invention  of 
advanced  lightweight  composite  materials  suddenly  pro¬ 
vided  a  solution  to  the  divergence  problem. 

Aeroelastic  divergence  occurs  when  the  wing  tip 
twists  upwards  as  dynamic  pressure  under  the  wing  increases 
This  upwards  twisting  increases  angle  of  attack  and  gener¬ 
ates  additional  torque  which  is  greater  than  the  structural 
elastic  restoring  forces.  As  a  result,  the  deformation 
continues  until  the  structure  fails  (6:147). 

Aeroelastic  Tailoring.  Grumman  material  engineers 
overcame  aeroelastic  deformation  by  using  a  new  process 
called  aeroelastic  tailoring.  Boron  fiber,  graphite, 
plastic  adhesives,  and  other  composite  materials  produce 
non-metallic  wings  30  percent  lighter  than  metal  wings  of 
similar  strength.  By  varying  the  fiber  direction  and  thick 
ness,  the  designer  tailors  the  wing  stiffness  character¬ 
istics  to  suit  a  particular . loading  condition,  distributing 


the  twisting  load  (greatest  at  the  wing  tips)  over  the 
entire  wing. 

Wave  Drag.  An  aircraft  flying  through  the  atmos¬ 
phere  produces  a  "pressure  wave"  as  the  fuselage,  wings,  and 
other  parts  of  the  vehicle  push  the  air  out  of  the  way. 
Pressure  waves  cannot  travel  faster  than  the  speed  of 
sound.  When  an  aircraft  flying  at  a  supersonic  velocity 
overtakes  the  pressure  wave,  the  vehicle  experiences  a 
phenomenon  called  "wave  drag."  The  magnitude  of  the  wave 
drag  is  reduced  by  using  a  minimum  and  uniform  aircraft 
cross-sectional  area  perpendicular  to  the  alrstream.  There¬ 
fore,  a  correctly  designed  fuselagt^  ^  "pinched  in"  where 
the  wing  attaches  to  the  fuselage  ''unt  for  the  added 

cross-sectional  area  of  the  wing.  Aft  swept  wings,  by 
necessity,  attach  to  a  "pinched  in"  fuselage  near  the  center 
of  gravity.  Unfortunately,  the  payload  also  is  concentrated 
near  the  aircraft's  center  of  gravity  to  minimize  the  pitch¬ 
ing  moment.  This  lowers  the  available  cargo  space  in  the 
vehicle.  Advantageously,  forward  swept  wings  connect  to 
the  fuselage  farther  aft  of  the  center  of  gravity  outside 
the  area  of  payload  concentration  (6tl48). 

Stall  Progression.  The  forward  swept  wing  provides 
the  x-29  with  the  capability  of  flying  at  low  airspeeds  and 
at  angles  of  attack  as  high  as  70  degrees.  This  is  pos¬ 
sible  because  stalls  progress  opposite  to  the  direction  of 
airflow  from  the  trailing  edge  to  the  leading  edge  and  from 
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the  most  aft  to  the  most  forward  part  of  the  wing.  There¬ 
fore,  aft  swept  wings  stall  first  at  the  wing  tips  near 
the  ailerons.  As  a  result,  the  ailerons  become  ineffective 
and  lateral  control  is  lost  at  approximately  30  degrees 
angle  of  attack.  Moving  the  ailerons  away  from  the  wing 
tips  is  not  practical  because  it  lessens  the  roll  moment 
due  to  a  shorter  lever  arm.  In  contrast,  for  the  forward 
swept  wing,  the  stall  begins  at  the  root  of  the  wing  away 
from  the  ailerons  and  progresses  toward  the  wing  tips. 

This  phenomenon  allows  greater  maneuverability  at  low  air¬ 
speed— a  valuable  aerial  combat  advantage. 

Super  Critical  Airfoil.  The  X-29  wings  are  not 
only  swept  forward,  but  also  "supercritical,"  which 
decreases  drag  and  minimizes  loss  of  lift  at  supersonic 
airspeeds.  Geometrically,  the  leading  edge  of  the  wing 
has  a  greater  radius  of  curvature  and  a  flatter  top  surface 
than  the  conventional  airfoil.  The  curvature  causes  air 
flowing  over  the  supercritical  wing  to  quickly  reach  super¬ 
sonic  speed.  The  flat  top  surface  sustains  the  supersonic 
velocity  of  the  air  so  the  shock  wave  forms  at  the  trailing 
edge  of  the  wing,  where  drag  is  lower  and  loss  of  lift  is 
less  (5:26). 

2-4  Canard 

The  canard,  located  just  behind  the  cockpit, 
increases  the  capability  of  the  forward  swept  wing  and  is 
the  primary  airfoil  for  generating  the  control  moments  and 
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stabilizing  the  aircraft's  longitudinal  axis  (see  Fig.  2-2). 
The  canard  and  wing  Interact  to  minimize  trim  drag  at  all 
mach  numbers.  Also  the  canard  directs  additional  air  to 
the  root  of  the  wing  to  further  resist  stall  at  high  angles 
of  attack. 

Since  the  X-29  is  statically  unstable;  the  center 
of  gravity  of  the  aircraft  is  aft  of  the  center  of  pressure. 
Therefore;  the  canards  must  generate  a  negative  lift  during 
one  "g"  trimmed  steady  level  flight  to  cancel  the  positive 
moment  produced  by  the  wing.  Fortunately;  during  high  "g" 
maneuvers;  the  canard  creates  a  positive  pitching  moment 
with  a  positive  aerodynamic  force  that  adds  to  the  total 

lift  on  the  vehicle.  Additionally;  this  arrangement  allows 

< 

greater  movement  of  the  center  of  gravity  and  enhances 
flexibility  in  payload  distribution. 

The  range  of  canard  surface  deflection  is  -60 
degrees  to  '*■30  degrees  with  a  rate  limit  of  105  degrees  per 
second  (14 1 18). 

2-5  Strake  Flaps 

The  strakes  in  Fig.  2-2  are  horizontal  airfoils 
connected  to  the  fuselage  from  the  wings  to  the  tail.  Pri¬ 
marily;  the  strakes  function  to  produce  a  nose  down  pitch¬ 
ing  moment  at  high  angles  of  attack.  The  flaps  at  the  end 
of  the  strakes  provide  the  control  power. 
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The  daflactlon  limits  for  the  itraka  flaps  are  -30 
degress  to  +30  degrees.  The  rate  limit  is  a  slow  27  degrees 
per  second  (15:3-47), 

2-6  Variable  Camber  Control  Surfaces 

Dual  hinged  flaperons,  attached  to  the  trailing 
edge  of  the  forward  swept  wing>  provide  high  lift  flaps, 
lateral  control,  and  discrete  variable  camber  to  maximize 
lift/drag  ratio  over  all  flight  conditions  (see  Fig.  2-2) . 

Dividing  the  flaperons  into  three  separate  sections 
along  the  wing  combats  aerodynamic  flutter.  The  inboard 
flaperon,  consisting  of  two  parts,  generates  additional 
lift  similar  to  a  standard  flap.  The  outboard  section 
operates  differentially  like  ailerons  to  provide  roll  con¬ 
trol.  All  three  sections  move  together  as  symmetrical 
flaperons  for  variable  camber. 

variable  camber  allows  a  flight  control  computer 
to  automatically  optimize  the  wing  shape  for  maximum  effi¬ 
ciency  during  sueady  level  cruise  flight.  In  the  "normal 
mode"  of  operation,  a  full  authority  digital  flight  control 
computer  trinut  not  only  the  variable  camber  flaperons,  but 
also  the  cau^rd  and  strake  flap  positions  for  optimal  air¬ 
flow  over  all  airfoils  to  minimize  drag  for  the  desired 
flight  lltion.  In  the  "reversion  mode"  an  analog  backup 
system  liits  the  canard  and  strake  flaps  to  zero  position 
and  adjusts  only  the  forward  swept  wing  camber  to  trim  the 
rti„craft  for  a  particular  flight  condition. 
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The  flaperons  deflect  In  the  range  of  '•19.5  degrees 
to  4*13  degrees  with  a  rate  limit  of  68  degrees  per  second 
(15;  3-44) . 

2-7  Summary 

The  X-29  advanced  technology  demonstrator  with  for¬ 
ward  swept  wings  and  high  negative  static  margin  aero¬ 
dynamics  provides  a  significant  advantage  In  performance 
and  overall  flexibility  as  follows; 

1.  Improved  area  rule  configurations 

2.  higher  volumetric  efficiencies 

3.  more  versatile  weight  and  balance  arrangements 

4.  Improved  pilot  visibility 

5.  smaller  aircraft  size— reduced  radar  profile 

6.  simpler  and  lighter  Installations  for 
hydraulics 

7.  more  compatible  wing  pylon  store  carriage 

8.  superb  maneuvering  over  a  wide  range  of  air¬ 
speeds  and  angles  of  attack 

9.  enhanced  fuel  efficiency  at  transonic  and 
supersonic  velocities 

The  added  features  of  this  aircraft  are  expected  to  revolu¬ 
tionize  air  combat  tactics  and  provide  a  superior  fighting 
capability. 


3.  Mathematical  Model  of  X-29  Aerodynamics 


3-1  Introduction 

The  development  of  control  laws  using  the  Porter 
design  method  requires  an  accurate  state  space  matrix 
representation  of  the  equations  of  motion  of  the  aircraft. 
The  equations  of  motion  are  a  mathematical  model  of  the 
physical  system  (plant)  that  relate  external  forces  and 
moments  to  the  linear  and  angular  accelerations  of  the 
vehicle  in  accordance  with  Newton's  laws  of  motion.  To 
derive  the  equations,  certain  assumptions  must  be  made  and 
an  appropriate  axis  system  defined. 

This  chapter  describes  the  origin,  axis  system, 
assumptions,  limitations,  and  dimensions  of  the  aircraft 
equations.  Also  included  in  the  chapter  are  the  mathemati¬ 
cal  equations  that  model  the  engine  thrust  and  control 
surface  actuator  dynamics.  Understanding  these  details 
is  critical  to  the  successful  design  of  a  control  system. 

3-2  Origin  of  Data 

Grumman  Aerospace  Corporation,  manufacturer  of  the 
X-29,  calculates  all  aerodynamic  data  including  stability 
derivatives  for  the  aircraft  using  wind  tunnel  tests  and 
other  techniques.  As  with  any  development  project,  they 
continually  update  and  upgrade  the  data  as  new  information 
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and  aircraft  modifications  occur.  In  late  1983/  Grumman 
supplied  aerodynamic  data  set  number  four  to  NASA  Ames 
Research  Center  at  Edwards  APB,  California.  On  24  May  1984, 
NASA  used  their  own  computer  program  (2)  and  Grumman 's 
aerodynamic  data  set  four  to  derive  the  six  degree  of 
freedom  state  space  matrix  eqtiations  used  in  this  thesis. 

3-3  Assumptions 

The  complexity  of  formulating  and  solving  problems 
in  flight  dynamics  necessitates  the  need  to  make  commonly 
accepted  approximations  involving  the  equations  of  motion. 
Judiciously  selected  assumptions  reduce  the  design  effort 
with  minimum  loss  of  accuracy. 

The  equations  of  motion  in  this  research  effort 
include  the  following  assumptions: 

1.  The  earth  is  an  inertial  reference.  This  is  a 
commonly  accepted  approximation  because  aircraft  maneuvers 
have  short  durations  compared  to  earth  rate.  Also  the 
flight  control  system's  sensors  are  too  insensitive  to 
detect  the  earth's  rotational  rate  and  corriolis  accelera¬ 
tion, 

2.  The  atmosphere  is  assumed  fixed  with  respect  to 
the  earth.  This  assumption  limits  the  scope  for  a  pre¬ 
liminary  design.  Final  designs  should  include  the  effects 
of  random  gusts,  wind  shears,  and  steady  winds. 

3.  The  aircraft  is  a  rigid  body.  Initial  designs 
generally  incorporate  this  assumption.  A  bending  mode 
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analysis  may  be  necessary  in  the  final  design  depending  on 
the  type  and  size  of  the  aircraft  and  the  maneuver  per¬ 
formed.  However,  limiting  the  flight  control  system  band¬ 
width  reduces  interaction  with  bending  and  vibration  modes 
which  generally  occur  at  higher  frequencies. 

4.  The  aircraft  mass  is  constant.  This  is  a  good 
approximation  during  a  10  to  15  second  maneuver. 

5.  The  airstream  surrounding  the  airplane  changes 
instantaneously  when  the  vehicle  is  disturbed  from  equi¬ 
librium.  At  transonic  and  supersonic  speeds  the  compres¬ 
sibility  drag  reduces  the  validity  of  this  assumption,  but 
it  is  a  good  first  /icproximation. 

6.  Aerodynamics  of  the  vehicle  are  fixed  at  the 
trim  point  for  a  particular  Mach  number  and  altitude. 
Although  perturbations  in  velocity  and  altitude  occur  dur¬ 
ing  the  iiiaheuvers,  the  stability  derivatives  change  slowly 
(except  in  the  transonic  region).  Therefore,  the  approxi¬ 
mation  is  valid  for  small  variations  from  the  reference 
point. 

7.  The  equations  of  motion  can  be  decoupled  into 
longitudinal  and  lateral  sets  of  equations.  The  only  cross 
coupling  term  in  the  equations  of  motion  for  the  X-29  is  a 
slight  influence  by  beta  (side  slip  angle)  in  the  pitch 
rate  equation.  This  is  safely  ignored  for  the  prelimi¬ 
nary  design  attempt  of  this  thesis. 
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8.  Linear  perturbation  equations  about  an  operating 
point  (flight  condition)  can  be  derived  from  the  nonlinear 
aircraft  dyneunlc  equations. 

3-4  Limitations 

Knowledge  of  the  limitations  of  the  aircraft  model 
ensure  that  commanded  alrcreft  maneuvers  remain  within  the 
accuracy  boundaries  of  the  equations.  Limitations  are  as 
follows: 

1.  The  equations  of  motion  are  written  In  the  body 
fixed  reference  frame  because  sensors  providing  feedback 
signals  measure  body  accelerations.  Also,  the  pilot  prefers 
to.  control  quantities  that  his  senses  can  easily  detect 
such  as  body  accelerations  and  rates.  Figure  3-1  Illus¬ 
trates  the  geometry  of  the  reference  frame.  The  origin  of 
the  body  axes  Is  at  the  mass  center  of  the  vehicle.  The 
axes  are  mutually  perpendicular.  The  x  and  z  axes  lie  In 
the  longitudinal  plane  and  the  y  axis  Is  In  the  lateral 
plane. 

2.  The  linear  perturbation  equations  are  derived 
from  the  nonlinear  equations  of  motion  of  an  aircraft  which 
Is  In  steady  level  trimmed  flight  at  a  particular  flight 
condition  (operating  point) .  Therefore,  solutions  to  the 
linear  equations  .represent  deviations  from  the  trimmed 
operating  point. 

3.  This  thesis  uses  the  "normal  mode"  equations. 

The  "normal  mode"  equations  of  motion  for  the  X-29  represent 
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an  optimum  steady  level  trimmed  condition  at  each  flight 
condition  which  minimizes  drag  on  the  aircraft  using  the 
canards,  symmetrical  variable  camber  wing  flaps,  and  the 
strake  flaps.  The  "reversion  mode"  equations  of  motion 
represent  a  trimmed  condition  using  only  the  symmetrical 
variable  camber  wing  flaps  with  the  canard  and  strake  flaps 
set  to  zero  deflection.  See  Chapter  2  for  additional 
details  on  the  function  and  operation  of  all  the  control 
surfaces. 

4.  The  linear  equations  of  motion  are  valid  for 
bank  angles  of  30  degrees  or  less  and  yaw,  pitch,  and  angle 
of  attack  angles  of  15  degrees  or  less. 

3-5  Dimensions 

The  linear  perturbation  equations  of  motion  fur¬ 
nished  by  NASA  (2)  have  the  following  variables  and  associ¬ 
ated  dimensions! 

1.  Pitch  angle  (9) — radians 

2.  Angle  of  attack  (a) — radians 

3.  Pitch  rate  (q) — radians  per  second 

4.  Change  In  forward  velocity  (u)--feet  per  second 

5.  Normal  acceleration  at  the  pilot  station 
(Ngp)— g' s 

6.  Normal  acceleration,  aircraft  center  of 
gravity  (N^g)— g's 

7.  Side  slip  angle  (B) — radians 

8.  Roll  angle  ((|)) — radians 


9.  Yaw  rate  (r) — radians  per  second 

10.  Roll  rate  (p) — radians  per  second 

11.  Lateral  acceleration  (Ny) — g's 

12.  .Canard  surface  deflection  (6^) --minus  60  degrees 
to  positive  30  degrees  (positive  deflection  produces  pitch 
up  moment/  negative  deflection  produces  pitch  down  moment) 

13.  Canard  surface  deflection  rate  --maxlmtim 
of  105  degrees  per>second 

14.  Strake  flap  surface  deflection  ( 6 g) --minus  30 
degrees  to  positive  30  degrees  (positive  deflection  produces 
pitch  down  moment/  negative  deflection  produces  pitch  up 
moment) 

a 

15.  Strake  flaps  deflectloxi  rate  (6 g) —maximum  27 
degrees  per  second 

16.  Thrust  (dip)— differential  throttle  angle  posi¬ 
tion  setting/  values  from  0.04  to  0.791  represent  idle  to 
military  thrust/  values  from  0.792  to  1.0  represent  after¬ 
burner  thrust.  Actual  values  of  engine  thrust  in  pounds 
of  force  are  not  defined. 

17.  Differential  wing  flap  (flaperon)  surface 
deflection  (6p) — minus  19.5  degrees  to  plus  13  degrees 

18.  Flaperon  surface  deflection  rate  (dp) — maximum 
of  68  degrees  per  second 

19.  Rudder  surface  deflection  (6p) — minus  31  degrees 
to  plus  31  degrees 
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20.  Rudder  surface  deflection  rate  — maximum 

of  141  degrees  per  second 

21.  All  the  equations  In  the  thesis  use  an  air¬ 
craft  weight  of  15000  pounds 

22.  Numbers  In  the  aircraft  equation  of  motion 
matrices  are  accurate  to  four  significant  figures 

3-6  Actuator  Dynamics 

The  equations  of  motion  for  the  X-29  aircraft, 
furnished  by  NASA,  do  not  Include  actuator  dynamics.  How¬ 
ever,  an  objective  of  the  thesis  Is  the  design  and  evalua¬ 
tion  of  a  flight  control  system  for  the  aircraft  with 
engine  thrust  and  control  surface  actuator  dynamics 
Included  In  the  complete  closed-loop  system. 

To  approximate  the  dynamics  of  the  canard  actuator 
and  engine  thrust  response,  this  research  effort  uses  first 
order  differential  equations  with  a  time  constant  of  0.05 
seconds.  It  Is  the  same  time  constant  utilized  by  Grumman 
engineers  in  their  Forward  Swept  Wing  Demonstrator  Tech¬ 
nology  Integration  and  Evaluation  Study  (14il8). 

The  strake  flap  actuator  time  constant  la  set  at 
0.20  seconds,  which  Is  four  times  longer  than  the  time 
constant  of  the  canard  actuator.  This  simulates  the  rate 
limit  of  27  degrees  per  second  for  the  strake  flap  which 
is  approximately  four  times  slower  than  the  maximum  deflec¬ 
tion  rate  of  105  degrees  per  second  for  the  canard. 


Using  the  above  time  constants  results  In  Laplace 
transfer  functions  as  follows: 

6^/e^  -  20/(3+20)  (3-1) 

Sg/Sg  -  5/{s+5)  (3-2) 

6^/e^  -  20/(s+20)  (3-3) 

where 


6^  ■  canard  deflection 
6g  ■  strake  flap  deflection 
6^  ■  thrust 

Bq  "  signal  Input  to  canard  actuator 
eg  ■  signal  Input  to  strake  flap  actuator 
e^  ■  signal  Input  to  engine  throttle  aotuator 

The  computer-aided  design  program  MULTI  (11) 
requires  the  above  equations  In  state  space  representation 
where 


and 


X  •  A  X  fa  ii 

■  -20  +  20  e^ 

is  -  “5  is  +  5  Sc 

It  "  “20  ix  20  e^ 

■  C  X 


(3-4) 

(3-5) 

(3-6) 

(3-7) 

(3-8) 
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3-7  Surotnarv 

The  assumptions  and  limitations  associated  with 
the  linear  perturbation  equations  of  motion  for  the  X-29 
aircraft  are  numerous  but  still  provide  an  adequate  model 
of  the  vehicle  Cor  the  initial  design  of  an  aircraft  gontrol 
system. 

The  next  chapter  describes  in  detail  the  design  of 

% 

the  longitudinal  axis  flight  control  system  which  includes 
the  aircraft  model  equations  and  control  laws. 


27 


4.  Longitudinal  Controller  Design  Procedure 


4-1  Introduction 

This  chapter  deaoribes  in  detail  the  procedure 
developed  in  thie  research  effort  and  subsequently  followed 
to  design  a  proportional  plus  integral  multivariable 
discrete- time  tracking  system  for  the  X-29  aircraft,  imple¬ 
mented  with  and  without  computer  time  delay  and  actuator 
dynamics.  The  comprehensive  design  procedure  applies  the 
Porter  design  method  and  includes  the  following  major 
stepsi 

1.  Validate  the  aircraft  model 

2.  Analyze  the  aircraft  open-loop  transfer 
functions 

3.  Select  the  required  output  variables  and  feed¬ 
back  signals 

4.  Check  controllability  and  observability 

5.  Calculate  the  transmission  zeros 

6.  Determine  measurement  matrix  if  required 

7.  Estimate  maximum  maneuvers  using  the  G(0)  matrix 

8.  Determine  appropriate  control  input  magnitudes 
and  rates 

9.  Specify  sigma  weighting  matrix  elements 

10.  Choose  the  epsilon  and  alpha  proportional  plus 


integral  gain  factors 


11.  Use  MULTI  tot 

« 

a.  calculate  the  and  matrices 

b.  form  the  closed-loop  control  system  to 
check  for  stable  roots  in  the  characteristic  polynomial 

c.  run  a  discrete  simulation  of  the  entire 
control  system  and  display  the  time  responses 

d.  optimize  the  design  with  actuators  and 
computational  time  delay 

Throughout  the  remainder  of  this  chapter/  sample 
calculations  for  the  flight  condition  of  0.7  maoh  at  15,000 
feet  illustrate  specific  details  in  each  step  of  the  design 
procedure  outlined  above.  Chapter  S  and  Appendices  D 
through  F  present  the  data  and  results  fox  the  other  flight 
conditions. 

4-2  Validation  of  Aircraft  Model 

Chapter  3  describes  in  detail  the  origin,  assump¬ 
tions,  limitations,  and  dimensions  associated  with  the 
state  space  representation  of  the  linear  perturbation  aqua¬ 
tions  of  motion  of  the  aircraft. 

This  research  effort  does  not  derive  the  equations 
of  motion  for  the  X-29  aircraft.  The  equations  are  a  pro¬ 
duct  of  NASA  Ames  Research  Center's  computer  program  "X-29 
Trim  and  Linearized  Data  Routine,  Forward  Swept  Wing  Linear 
Six  Degree  of  Freedom  Simulation"  (2) .  The  data  represents 
the  fourth  update  and  refinement  and  includes  Grumman' s 
aerodynamic  data  set  number  four. 
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A  typical  axampla  of  tha  information  providad  by 
NASA  ii  tha  A  and  B  matrix  longitudinal  axis  aquations 
for  0.7  maoh  at  15,000  fast  as  shown  balowt 


X 


s  « 
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• 

u 
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a 

a 

q 

■  a 

mm 

0 

0 

0 
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-32.13 

-0.01199 
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u 

-0.3246E-06 

-0.1298E-03 

-1.138 

0.9917 
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0.298621-03 
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18.92 

-0.4928 
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q 

*  ■ 
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-1 

r 

0 
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0 

-0.05895 

0.02621 

27.24 

-0.001161 

-0.7808B-03 

-0,002099 

«s 

0.1720 

-0.05630 

-0.1929 

- 1 

(4-1) 


This  thasis  assumas  tha  numarioal  valuas  in  tha  matrioas 
ara  oorraot.  Obviously,  tha  aquations  ralating  klnamatio 
ralationships  ara  aoourata.  In  addition,  a  ohaok  of  all 
tha  signs  on  tha  opan-loop  transfar  functions  for  tha 
unaugmantad  aircraft  indioata  tha  oorraot  diraotion  of 
movamant  of  tha  vahiola  for  a  glvan  control  daflaotion. 

4-3  Analysis  of  Aircraft  Modal 

An  analysis  of  tha  opan-loop  transfar  functions  of 
tha  unaugmantad  aircraft  provida  insight  into  tha  possible 
difficulty  in  designing  a  satisfactory  control  system.  As 
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described  In  Chapter  2-2  "Airframe,"  the  X-29  aircraft 
Is  statically  unstable  at  subsonic  mach  numbers.  The  fol¬ 
lowing  open-loop  transfer  functions  for  0.7  mach  at  15000 


feet  clearly  Illustrate  this  phenomenal 
Characteristic  Polynomial  (C.P.) 

-  (S  +  5.16)  (8  -  3.529)  (8  +  0.005907  1  j  0.06082)  (4-2) 
q/6Q  -  0.172  s  (8  +  1.01)  (s  +  0.0125) /C.P.  (4-3) 
q/6g  -  -0.0563  s  (s  +  1.4)  (s  +  0. 01254) /C. P.  (4-4) 
q/6^  -0.1929  s  (s  0.9985)  (s  -t-  0.3152)/C.P.  (4-5) 


u/6,.  -  -0.05895  (s  -H  0.7165  t  j  8.906)  (s  1.186) /C.P. 

^  (4-6) 

u/6o  ■  -0.02621  (s  +  0.3287  1  j  7.367)  (s  +  1.777) /C.P. 

®  (4-7) 

u/.;^  -  27.24  (8  •«-  5.139)  (s  -  3.489)  (s  -  0. 01705) /C.P.  (4-8) 

a/6r  ”  "0. 001161  (s  -  146.4)  (s  0.005991  i  j  0. 06417/C. P. 

(4-9) 

a/6c  "  -0.0007808  (s  +  72.0)  (s  +  0.005934  1  j  0. 065380) /C.P. 
®  (4-10) 

a/6m  ■  -0.02099  (s  +  93.35)  (a  -  0.01025  ±  j  0. 06406) /C. P. 

^  (4-11) 

There  is  an  unstable  root  at  3.529  in  the  characteristic 
polynomial.  Also,  several  of  the  transfer  functions  have 
zeros  in  the  right  half  plane.  As  an  illustration, 
examine  the  transfer  function  for  a/6^.  From  a  single¬ 
input-single-output  (SISO)  classical  root  locus  point  of 
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vlew/  the  zero  at  146.4  and  the  pole  at  3.529  depict  a 
system  that  Is  highly  unstable  for  all  positive  values  of 
feedback  gain.  Without  compensation,  the  closed-loop  sys¬ 
tem  characteristic  polynomial  has  a  root  between  3.529  and 
146.4  depending  on  the  value  of  the  feedback  gain.  This 
degree  of  Instability  in  an  aircraft  is  very  unusual.  The 
plant  is  a  non-minimum  phase  system.  Therefore,  the  suc¬ 
cessful  design  of  a  control  system  is  difficult.  This  is 
especially  true  when  actuator /sensor  dynamics  and  computa¬ 
tional  time  delay  are  added  to  the  basic  equations.  These 
additions  may  reduce  the  phase  margin  and  make  it  a  harder 
task  to  assure  stabilization  of  the  system. 

4-4  Controlled  Variables 

The  output  vector  '’j'”  represents  the  controlled 
variables.  The  outputs  can  be  states  of  the  aircraft  plant 
or  linear  combinations  of  states.  The  inputs  to  the  con¬ 
trol  system,  represented  by  the  vector  " v" ,  are  the  com¬ 
mands  to  be  tracked  by  the  output  vector  y  (see, Fig.  A-1). 
Examples  include  pitch  rate  and  angle  of  attack.  Inputs 
to  the  aerodynamic  model  (plant) ,  represented  by  the  vector 
"u" ,  are  the  controllable  mechanisms  on  the  vehicle  that 
can  modify  its  motion.  For  the  x-29  longitudinal  axis, 
the  inputs  include  the  canards,  strake  flaps,  thrust,  and 
symmetrical  flaperons.  For  the  lateral  axes,  the  inputs 
are  the  rudder  and  flaperons. 
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The  Porter  design  technique  requires  an  equal 
number  of  Input  and  output  (controlled)  variables.  Addi¬ 
tionally,  the  Porter  method  incorporates  output  feedback 
and  all  outputs  are  fed  back  to  the  Input  of  the  control 
system.  See  Appendix  A  for  more  theory  and  block  diagrams 
of  the  high-gain  error-actuated  closed-loop  control  system 
configuration. 

Each  Independently  controlled  variable  requires 
a  separate  control  input  for  a  decoupled  and  completely 
controllable  vehicle  response.  This  Is  a  basic  require¬ 
ment  which  is  independent  of  the  design  method.  For 
example,  controlling  pitch  angle,  velocity,  and  angle  of 
attack  requires  three  control  mechanisms  in  the  longitudi- 
nal  axis  of  the  aircraft  (canards,  strake  flaps,  and  engine 
thrust  for  the  X-29). 

Selecting  the  correct  control  output  variables  is 
very  important  and  depends  on  several  conditions  outlined 
below: 

1.  the  available  aircraft  control  inputs 

2.  practicality  and  reliability  of  measuring  the 
controlled  variables 

3.  the  desired  maneuvers 

4.  restraints  of  the  design  method 

Aircraft  Controls.  Canards,  symmetrical  flaperons, 
strake  flaps,  and  thrust  are  the  longitudinal  axis  aircraft 
control  inputs  on  the  X-29.  These  four  con:  i:ol  inputs 
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ovide  the  capability  to  command  four  independent  long!** 
tudinal  outputs.  It  is  difficult  to  command  four  separate 
outputs  that  do  not  conflict  with  each  other.  Therefore, 
the  designs  in  this  thesis  control  three  longitudinal  vari¬ 
ables  using  the  canards,  strakes  flaps,  and  thrust  as  the 
three  control  inputs.  The  sytmnetrlcal  flaperons  are  not 
used  as  a  control  input  because  their  primary  function  is 
trimming  the  aircraft  to  minimize  drag  at  each  cruise 
flight  condition  by  optimizing  the  air  flow  over  the 
canards,  wings,  and  strakes.  Since  this  research  effort 
uses  the  "normal  mode"  equations  of  motion,  the  symmetrical 
flaperons  are  already  set  to  the  optimum  position  for 
minimum  drag.  Also,  the  maneuvers  are  small  perturbations 
about  the  trimmed  operating  point,  thereby  minimizing  the 
need  to  move  the  symmetrical  flaperons. 

Measurement  of  Output. Variables.  The  pilot  feels 
angular  rates  and  accelerations  and,  in  actuality,  he  con¬ 
trols  those  quantities  when  flying  conventional  aircraft. 

In  addition,  it  is  more  reliable  to  measure  rates  and 
accelerations  using  rate  gyros  and  accelerometers  than  to 
measure  pitch  and  bank  angles.  The  pitch  and  bank  angle 
measurements  are  generally  obtained  from  an  ir.crtial  navi¬ 
gation  system  with  significantly  lower  hardware  reliability 
than  the  flight  control  system.  Angle  of  attack  measure¬ 
ments  are  reasonably  reliable  but  are  subject  to  wind  gusts 
and  turbulence. 


34 


This  thesis  Is  successful  in  performing  several 
Important  longitudinal  maneuvers  using  pitch  rate^  velocity, 
and  angle  of  attack  as  the  output  variables.  An  attempt  to 
control  normal  acceleration  was  unsuccessful  due  to  current 
limitations  with  the  design  method. 

Desired  Longitudinal  Maneuvers.  The  multiple  con¬ 
trol  surfaces  In  the  longitudinal  axis  allow  the  X-29  to 
generate  "direct  forces,"  l.e.,  aircraft  movements  uncoupled 
or  not  related  to  aircraft  rotation.  "Direct  force"  in  the 
longitudinal  axis  makes  possible  three  new  maneuvers  called 
"direct  lift,"  "pitch  pointing,"  and  "vertical  translation." 

The  "direct  lift"  mode  produces  a  longitudinal 
acceleration  while  holding  angle  of  attack  constant.  The 
aircraft  flight  path  angle  changes  due  to  additional  lift 
generated  without  changing  angle  of  attack.  To  accomplish 
this  maneuver  the  flight  control  system  commands  pitch  rate 
and  velocity  to  a  desired  value  and  angle  of  attack  to 
zero.  A  variation  to  the  "direct  lift"  mode  is  the  "direct 
climb"  whereby  pitch  i.^te  is  commanded  as  a  pulse  of  length 
sufficient  to  generate  a  desired  pitch  angle.  The  flight 
path  angle  equals  the  new  pitch  angle  since  alpha  is  com¬ 
manded  to  zero. 

The  "pitch  pointing"  mode  provides  the  pilot  the 
ability  to  command  changes  in  angle  of  attack  while  main¬ 
taining  a  constant  aircraft  flight  path  angle  and  velocity. 
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Commanding  a  pitch  rate  pulse  that  creates  a  pitch  angle 
equal  to  angle  of  attack  accomplishes  the  maneuver. 

"Vertical  translation"  allows  the  aircraft  to 
change  altitude,  i.e.,  translate  vertically,  without 
rotating  the  aircraft  or  changing  its  forward  velocity. 

This  maneuver  requires  pitch  rate  commanded  to  zero,  no 
change  in  forward  velocity,  and  a  negative  angle  of  attack 
to  translate  the  aircraft  upwards  and  positive  angle  of 
attack  to  translate  the  aircraft  downwards. 

The  ability  of  the  aircraft  to  perform  any  of  the 
"direct  force"  maneuvers  depends  on  the  force  generating 
capability  of  the  control  surfaces  and  their  deflection 
limits. 

Design  Method  Restraints.  Selecting  the  proper 
output  variables  depends  not  only  on  the  desired  aircraft 
maneuvers  but  also  the  restraints  of  the  design  method. 

For  example,  the  design  method  requires  the  number  of  out¬ 
puts  to  equal  the  number  of  inputs.  Furthermore,  the  tech¬ 
nique  mandates  that  all  input  variables  be  commanded  to  a 
specific  value,  including  zero.  Consequently,  some 
maneuvers  require  precomputations  to  determine  the  correct 
input  commands.  For  example,  a  constant  "g"  normal  accelera¬ 
tion  maneuver  requires  a  precomputed  pitch  rate  and  angle  of 
attack  command.  Still  other  maneuvers  are  difficult  to 
command.  As  an  example,  if  pitch  rate,  velocity,  and  angle 
of  attack  are  the  output  variables  and  the  desired  maneuver 


is  a  constant  thrust  climb;  the  Porter  control  system  needs 
a  velocity  input  command.  Precomputing  the  change  in 
velocity  is  difficult.  In  this  case/  changing  the  output 
matrix  to  remove  velocity  is  more  efficient  and  eliminates 
the  need  for  precomputations.  However,  modifying  the  out** 
put  matrix  does  introduce  additional  problems  with  imple¬ 
mentation  of  the  system. 

Careful  selection  of  output  variables  is  necessary 
t(.  design  an  efficient  flight  control  system. 

4-5  Controllability  and  Observability 

Properties  of  the  state  space  representation  of 
the  system  determine  controllability  and  observability.  A 
completely  controllable  system  implies  the  numerical  value 
of  each  state  variable  x(t)  at  any  initial  time,  tQ,  can 
be  changed  to  any  other  numerical  value  in  a  finite  time  by 
means  of  an  unconstrained  control  input  vector  u(t)  (3:442). 
For  a  time  variant  linear  system,  complete  controllability 
requires  the  controllability  matrix  to  have  full  rank 
where 

Rank  -  Rank  [B  M  ....  -  n  (4-12) 

A  system  is  completely  observable  if  every  state 
affects  the  output  y  (t) .  This  implies  that  the  numerical 
value  of  every  state  x(t)  can  be  determined  from  the  mea¬ 
surements  of  the  output  vector  y(t)  (3:442).  For  a  time 
invariant  system,  the  observability  matrix  Mq  must  have 


full  rank  where 

rank  Mq  -  Rank  [C^  aV  ...  -  n  (4-13) 

Also  a  system  is  completely  controllable  and  observ¬ 
able  if  there  are  no  input  and  no  output  decoupling  zeros. 

In  lieu  of  directly  determining  the  rank  of  the  and  Mq 
matrices#  this  thesis  uses  an  alternative  computer  program 
called  "Zero"  to  calculate  the  transmission#  input 
decoupling#  and  output  decoupling  zeros  of  the  state  and 
output  aircraft  plant  equations  (7). 

In  this  research  effort#  the  X-29  aircraft  open- 
loop  plant  is  completely  controllable  and  observable#  i.e.# 
there  are  no  input  or  output  decoupling  zeros#  for  all  the 
flight  conditions  and  output  variable  configurations  (C 
matrices) . 

Likewise#  a  check  of  the  controllability  and  observ- 
eddility  of  the  closed-loop  tracking  system  augmented  with 
integrators  is  very  important.  A  sampled  data  system 
augmented  with  a  "Porter  configured"  proportional  plus 
integral  digital  controller  is  completely  controllable  and 
observable  when  the  open-loop  plant  is  controllable  and 
observable  and 

■  number  of  states  (4-14) 

4-  number  of  outputs 


rank 


4)(T)-I  i|)(T) 
-Tr  0 


Multiplying  this  matrix  by  the  frequency  "f"  does  not 
change  the  rank  of  the  matrix.  For  the  case  when  "f" 
approaches  Infinity,  this  becomes 


^(T)-I  i|j(T) 

A  B 

rank  fi 

■  rank 

-Tr  0 

-C  0 

as  shown  by  Bradshaw  and  Porter  (10i822). 

This  satisfies  the  definition  of  total  controllabil¬ 
ity  whereby  every  bounded  Input  yields  bounded  outputs  and 
states. 

However,  when  pitch  rate  Is  an  output  variable,  the 
above  matrix  Is  rank  deficient  and  the  closed-loop  tracking 
system  does  not  satisfy  the  strict  definition  of  controlla¬ 
bility.  But  the  system  Is  still  "conditionally”  control¬ 
lable.  The  clroumstanoe  often  occurs  when  rates  and  accel¬ 
erations  are  the  output  variables.  The  cause  Is  a  transmis¬ 
sion  zero  at  the  origin  that  makes  pitch  angle  unbounded 
with  a  constant  pitch  rate  Input.  This  is  desirable  for 
fighter  aircraft  as  explained  In  more  detail  In  the  next 
topic . 

4-6  Transmission  Zeros 

By  definition,  a  transmission  zero  of  a  transfer 
function  Gis)  blocks  the  transmission  of  an  exponential 
signal  from  the  Input  to  the  state  variable  X(s) .  For 
example  If 
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(4-16) 


u(t)  - 

then 

x(t)  #  (4-17) 

As  discussed  in  Appendix  A»  Multivariable  Control  Law 
Theory,  the  tranamission  zeros  of  a  system  are  regions  to 
which  certain  finite  (slow)  roots  of  the  closed  loop  char¬ 
acteristic  equation  migrate  with  increasing  gain.  As  the 
system' s  closed  loop  gains  approach  infinity  or  sampling 
time  approaches  zero  for  discrete  systems,  the  affected 
roots  asymptotically  approach  the  transmission  zero  loca¬ 
tions.  Output  feedback  does  not  alter  the  number  or  loca¬ 
tion  of  transmission  zeros. 

Multivariable  control  theory  requires  transmission 
zeros  in  the  left  half  Laplace  s-plane  to  insure  stability 
at  high  gains.  However,  this  does  not  guarantee  stability 
for  gains  leas  than  infinity  because  the  locus  of  roots  may 
travel  into  and  out  of  the  unstable  region  before  approach¬ 
ing  the  locality  of  the  transmission  zeros.  Therefore, 
systems  with  stable  transmission  zeros  are  theoretically 
guaranteed  stable  only  for  infinite  gains.  However,  for 
most  systems,  especially  minimum  phase  plants,  gains  much 
lower  than  infinity  result  in  stable  systems  with  excellent 
performance.  Thus,  control  is  achieved  with  gains  of 
reasonable  magnitude. 


A  tucoeeiful  controller  design  for  a  non-mlnlmum 
phase  plant  with  unstable  transmission  zeros  may  be  pos¬ 
sible  with  the  Porter  method.  In  such  cases  gain  bound¬ 
aries  must  be  established  where  the  locus  of  roots  are  in 
the  stable  region.  Unfortunately,  there  is  no  theoretical 
method  to  determine  the  gain  boundaries  and  the  system  may 
not  be  stable  for  any  value  of  gain.-  Currently,  trial  and 
error  is  the  only  technique  available.  Although  there  may 
be  a  range  of  gains  for  stable  operation,  finding  the  gain 
values  can  be  very  difficult  and  time  consuming. 

The  X-29  is  a  non-minimum  phase  plant  with  open- 
loop  zeros  and  poles  in  the  right  half  s  plane  as  detailed 
earlier  in  this  chapter.  For  example,  the  unaugmented 
open-loop  aircraft  plant  without  control  surface  actuators 
for  0.7  maoh  at  15000  feet  has  a  transmission  zero  at  the 
origin.  Therefore,  unstable  roots  in  the  closed-loop  charac¬ 
teristic  equation  are  easily  possible  depending  on  the 
values  of  the  gains. 

Theoretically,  a  transmission  zero  at  the  origin 
is  undesirable  if  the  control  system  requires  "bounded 
input-bounded  output"  controllability.  The  structure  of 
the  X-29  aircraft  state  equations  and  the  desire  to  control 
pitch  and  roll  rates  cause  transmission  zeros  to  occur  at 
the  origin.  This  means  thab  a  bounded  input  may  result  in 
an  unbound  output.  For  example,  in  this  investigation, 
the  output  variables  are  pitch  rate,  velocity,  and  angle  of 
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attaalc)  and  the  atate  variable!  are  pitch  angle >  velocity, 
angle  of  attack,  and  pitch  rate.  A  pitch  rate  step  command 
result!  in  the  pitch  angle  ramping  to  infinity.  This  Is 
predictable  and  normal  since  the  pitch  angle  Is  simply  the 
Integral  of  the  pitch  rate.  Hence,  In  fighter  aircraft  with 
"rates"  as  Inputs  and  outputs,  the  pilot  achieves  the 
desired  pitch  angle  by  simply  commanding  a  pitch  rate  pulse. 

4-7  Measurement  Matrix 

If  the  ^  does  not  have  full  rank,  the  plant  is 
"Irregular"  and  the  control  system  requires  a  measurement 
matrix.  See  Appendix  A  for  an  explanation  of  Irregular 
designs.  The  measurement  matrix  Introduces  additional 
transmission  zeros,  which  can  be  changed  by  altering  the 
numerical  values  of  the  measurement  matrix  elements.  Fur¬ 
thermore,  a  sparse  measurement  matrix  generates  transmis¬ 
sion  zeros  In  the  s-plane  which  are  the  Inverse  of  the 
measurement  matrix  elements.  As  an  example,  a  0.2  element 
In  the  measurement  matrix  results  In  a  transmission  zero 
at  -5.0.  But  the  relationship  Is  true  only  for  a  sparse 
measurement  matrix. 

The  control  engineer  should  make  an  effort  to 
eliminate  the  need  for  a  measurement  matrix  by  reordering 
the  equations  or  selecting  different  outputs,  if  possible. 

A  measurement  matrix  is  always  necessary  if  control  surface 
actuator  equations  are  Incorporated  as  part  of  the  A  matrix 
because  of  the  structure  of  the  equations. 


The  longitudinal  deiigne  in  thlo  reiearoh  effort 
require  no  measurement  matrix  beoause  the  outputs  are  pitch 
rate,  velocity,  and  angle  of  attack.  The  first  Markov 
parameter  is  satisfied  since  the  matrix  CB  has  full  rank. 

In  preliminary  designs,  not  presented  in  this  report,  the 
output  variables  include  pitch  angle.  This  results  in  a 
rank  defioienoy  of  one  for  CB  and  a  need  for  a  measurement 
matrix  with  a  minimum  of  one  non-zero  element. 

4-8  Steady-State  Transfer 
Funotlon 

The  maximum  magnitude  of  any  steady- state  maneuver 
is  dependent  on  the  aircraft  plant  and  is  independent  of 
the  control  system.  The  largest  possible  maneuver  depends 
on  the  aircraft  dynamios,  structural  limitations,  and  con¬ 
trol  surface  deflection  limits.  The  G(0)  matrix  for  the 
open-loop  aircraft  plant  is  useful  for  evaluating  maximum 
commanded  inputs  based  on  the  control  surface  deflections 
limits.  The  derivation  of  the  G(0)  matrix  begins  with  the 
state  spaas  matrix  representation  of  the  aircraft  open  loop 
plant  where 

i  ■  ^  +  Bu  (4-18) 

y  “  Cx  (4-19) 

Solving  for  y(s)  gives 

Y(8)  -  (C[8l  -  U(s)  (4-20) 
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wh«r« 


■\V 


G(B)  -  (C[«I  -  A]"^B) 


(4-21) 


Solving  for  U(b)  ronulti  In 


U(«)  ■  G{i)’^y(i) 


(4-22) 


Using  the  Lapleoe  final  value  theorem 


u(t),,  -  lim  [s  U(s)]  -  lim  [•  G(s)“^Yli)]  (4-23) 


§■+•0 


when  2^(t)  is  a  constant  vector 


V(b)  -  K/s 


(4-24) 


and  the  free  s  in  thci  numerator  and  denominator  oanosl, 
thus  the  final  result  is 


u(t),,  ■  G(0)“^K 


(4-25) 


For  the  X-29  at  0.7  mach  and  ISOOO  feet  with  the  pitch 
angle#  velocity#  and  angle  of  attack  as  the  output  vari¬ 
ables;  the  steady  state  transfer  function  isi 


G(0) 


-1 


0.1915 

>3.469 

1.183 


-0.03771 

-0.1150 

-0.0004662 


-395.0 

-870.7 

-0.0496 


(4-26) 


Calculations  using  the  above  relationship  show  that  the 
maximum  angle  of  attack  command  is  two  degrees  for  vertical 
translation  without  exceeding  control  deflection  limits. 
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For  the  pitch  pointing  maneuver  the  maximum  pitch  angle  and 
angle  of  attack  command  la  also  two  degrees.  Furthermore, 
the  GiO)"^  matrix  for  this  aircraft  demonstrates  that  a 
very  small  change  in  steady- state  angle  of  attack  commands 
very  large  changes  in  the  control  surface  deflections. 

This  is  easily  seen  by  observing  the  relatively  large 
numbers  in  the  third  column  of  the  0(0)**^  matrix.  The 
relationship  is  valuable  in  the  design  phase  when  selecting 
the  sigma  element  that  controls  the  angle  of  attack.  The 
element  value  must  be  small  to  prevent  excessive  control 
surface  deflections. 

Another  use  of  the  6(0)“^  matrix  allows  the 

designer  to  oross-oheok  the  validity  of  the  steady  state 

» 

closed- loop  control  system  time  plots.  After  the  comple¬ 
tion  of  the  transients,  the  final  values  of  the  control 
surface  positions  and  the  thrust  magnitude  on  the  response 
plots  should  agree  with  those  calculated  with  the  G(0)~^ 
matrix  in  equation  (4-25) . 

Appendix  B  documents  an  improvement  to  MULTI  that 
provides  the  capability  to  calculate  the  G(0)  and  GCO)**^ 
matrices  directly  from  the  A,  B,  and  C  matrices  of  the 
plant.  Also,  the  user  is  given  the  option  to  enter  the 
plant  output  vector  "y"  and  have  MUI^TI  automatically  calcu¬ 
late  the  input  vector  "u"  using  the  relationship 

a(t»3,  Z<t)„  (4-27) 


‘.y 

f>y. 
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where 


(4-28) 


^(t)  ■  a  constant  vector 

8  8 

The  user  may  also  enter  the  plant  input  vector  "u"  and  have 
MULTI  automatically  calculate  the  plant  output  vector  "y" 
using  the  equation 


where 


i:(t)„ 


G(0)  U(t)^g 


u(t)  ■  a  constant  vector 

8  8 


For  example,  the  direct  climb  maneuver  for  0.7  mach 
at  15000  feet  has  command  inputs  which  generate  7.5  degree 
pitch  and  flight  path  angles  with  zero  change  in  velocity 
and  angle  of  attack.  For  these  conditions,  the  0(0)'^ 
matrix  calculates  steady  state  values  of 


U''  » 

i*# 


canard  deflection  ••  0.02  degrees 
strake  flaps  deflection  ■  -0,44  degrees 

thrust  ■  0.15  throttle  angle  ratio 


i 


These  values  agree  with  the  time  response  plots  of  the 
closed-loop  tracking  system. 

Often  the  G(0)”^  matrix  is  singular  when  rates  and 
accelerations  are  output  variables.  As  an  example,  when 
pitch  rate  is  an  output  variable  the  G(0)  matrix  is  singular 
(no  inverse  is  possible)  because  pitch  rate  is  the  deriva¬ 
tive  of  pitch  angle.  An  extra  "free"  s  occurs  in  the 


,■  * 
Ik.. 


W- 
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numerator  of  the  corresponding  row  of  the  G(s)  matrix. 

That  s  does  not  cancel  when  applying  the  Laplace  final 
value  theorem  and  the  outcome  is  an  all  zero  row.  As  a 
consequence  of  the  structure  of  the  state  equations ,  the 
alrqraft  control  surfaces  do  not  In  this  case  directly  con¬ 
trol  pitch  rate  because  It  Is  kinematically  related  to 
pitch  angle,  l.e. , 

i-2+0*i^+0*ig+0*iT 

4-9  Control  Input  Magnitudes 

and  Rates 

The  Porter  design  method  creates  a  controller  which 
tracks  the  Input  commands.  For  step  Inputs,  1:he  controller 
rapidly  drives  the  control  surfaces  to  maximum  deflection 
at  very  high  rates.  Actuators  may  reduce  the  surface  rates, 
but  not  to  magnitudes  below  the  specification  limits.  How¬ 
ever,  ramping  the  Input  command  to  the  desired  value  slows 
down  the  rate  of  control  surface  deflection  to  acceptable 
limits.  The  MULTI  computer-aided  design  program  provides 
the  capability  to  represent  command  Inputs  as  ramp-to-hold 
signals.  For  many  of  the  maneuvers  In  this  research  effort 
the  Input  command  ramp  time  durations  are  calculated  by 
dividing  the  maximum  surface  deflection  position  by  the 
minlmvim  surface  rate.  For  example,  most  of  the  X-29  input 
signals  are  ramped  over  a  1.1  second  interval  as  estab¬ 
lished  by  the  strake  flaps  position  and  rate  limits. 
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In  addition  to  the  rate  limit  restrictions,  the 
maximum  magnitude  and  duration  of  the  Input  commands  must 
keep  the  aircraft  motion  within  the  linear  region  of  the 
aircraft  equations  of  motion.  For  example,  in  this  Invest! 
gatlon,  pitch  rate  is  an  Input  command  and  Its  magnitude 
and  the  pulse  duration  determine  the  final  pitch  angle  for 
the  direct  climb  maneuver. 

Furthermore,  the  pitch  rate  pulse  rotates  the  nose 
of  the  aircraft  and  produces  normal  acceleration  of  the 
vehicle  during  the  pull  up.  The  relationship 

q  -  (1845/U)  *  (4-30 

where 

q  *  pitch  rate  In  degrees  per  second 
U  ■  forward  velocity  in  feet  per  second 
■  normal  acceleration  in  "g's" 

calculates  the  correct  pitch  rate  for  a  desired  "g"  level 
(1:107).  For  the  direct  climb  maneuvers  in  this  thesis, 
the  peak  magnitude  of  the  pitch  rate  pulse  generates  a 
change  of  one  "g"  for  the  flight  condition  of  0.7  Mach  at 
15000  feet  as  follows: 

q  ■  (1845  /  740  ft  per  sec)  *  1  g  »  2.5  deg  per  sec  (4-31 

Since  the  aircraft  is  already  at  one  "g”  in  steady  level 
flight,  an  increase  of  one  "g"  creates  the  equivalent  of  a 


two  g"  normal  pull  up  acceleration  during  the  rotation 
phase  of  the  direct  climb  maneuver.  The  two  "g”  pull  up 
is  sufficient  to  evaluate  the  controller's  stability  aug¬ 
mentation  capability  in  the  longitudinal  mode.  For  com¬ 
parison  purposes,  the  same  input  command  is  used  for  the 
other  flight  conditions. 

For  the  vertical  translation  and  pitch  pointing 
maneuvers  the  limiting  factor  is  not  the  linear  range  of 
the  aircraft  equations  of  motion,  but  the  control  surface 
deflection  limits.  The  0(0)“^  matrix  is  used  to  calculate 
the  maximum  possible  control  inputs  for  the  maneuvers  at 
steady  state  without  exceeding  the  control  surface  position 
limits.  These  input  magnitudes  need  adjustment  only  if 
the  control  surface  deflections  overshoot  their  maximum 


steady  state  values  during  the  transient  response. 

4-10  Sigma  Weighting  Matrix 

The  elements  of  the  sigma  weighting  matrix  assign 
the  first-order  roots  which  approach  infinity  with  increas¬ 
ing  gain.  In  addition,  as  the  gain  nears  infinity,  the 
slow  modes  vanish  because  they  become  asymptotically 
uncontrollable  and  unobservable.  This  leaves  only  the  fast 
modes  (infinite  roots) .  In  addition,  if  the  sigma  weight¬ 
ing  matrix  is  also  diagonal,  the  system  exhibits  increas¬ 
ingly  tight  and  non-interacting  behavior  as  the  gain  of  the 
closed-loop  system  approaches  infinity.  See  Appendix  A 
for  more  details. 


From  a  practical  Implementation  point  of  view,  the 
values  of  the  sigma  weighting  matrix  elements  ultimately 
determine  the  time  responses  of  the  outputs.  Each  diagonal 
element  affects  the  time  response  of  the  corresponding  out¬ 
put.  As  an  example,  the  X-2d  longitudinal  mode  controller 
has  the  following  outputs: 

y^^  -  pitch  rate  (q) 

y2  "  forward  velocity  (u) 

y^  ■  angle  of  attack  (a) 

Therefore,  the  sigma  matrix  has  a  3  x  3  dimension.  Spe¬ 
cifically,  the  diagonal  element  in  row  one  of  the  sigma 
matrix  affects  the  time  response  of  the  pitch  rate,  the 
diagonal  element  in  row  two  modifies  the  forward  velocity 
response,  and  the  third  row  diagonal  element  influences 
the  angle  of  attack  response.  In  general,  Increasing  the 
value  of  a  sigma  element  causes  a  faster  time  response  for 
the  corresponding  output.  However,  for  the  aircraft  to 
respond  more  rapidly,  the  controller  commands  greater  con¬ 
trol  surface  deflection  positions  and  faster  rates  which 
may  exceed  the  design  specifications. 

For  a  minimum  phase  plant  which  is  statically 
stable,  the  initial  selection  of  the  weighting  matrix  ele¬ 
ments  for  a  preliminary  design  is  easy  and  most  designers 
begin  by  assigning  the  value  of  one  for  all  the  elements. 
After  completing  the  computer  simulation,  the  engineer 


emalyzes  the  transient  responses  and  Iterates  the  appropri¬ 
ate  sigma  element  magnitudes  to  Improve  undesirable  results. 

Unfortunately/  the  X-29  aircraft  Is  a.  non-mlnlmvim 
plase  plant  with  static  Instability.  In  fact,  starting 
with  a  value  of  one  for  all  the  sigma  weighting  matrix  ele¬ 
ments  gives  an  unstable  response  for  all  flight  conditions. 
Thus  a  lengthy  trial  and  error  process  begins  to  locate 
sigma  values  which  produce  stability.  An  analysis  of  the 
steady  state  G(0)~^  matrix  provides  some  Insight  on  the 
relative  magnitude  of  the  elements  that  may  produce  a  stable 
system.  As  discussed  earlier  in  this  chapter,  the  X-29 
control  surface  deflections  are  very  large  for  small 
angle  of  attack  changes.  Therefore,  the  sigma  element 
controlling  the  angle  of  attack  must  be  small  to  achieve 
a  stable  system  with  acceptable  control  surface  deflections. 
For  Illustration,  the  following  sigma  weighting  matrix  pro¬ 
vides  the  best  direct  climb  time  response  for  the  X-29  at 
0.7  Mach,  15000  feet  without  control  surface  actuators  or 
computer  time  delay: 


1.7  0  0 


L 


0 

0 


1.2  0 

0  0.01 


4-11  Proportional  Plus  Integral  Gains 

The  most  efficient  technique  to  determine  the  best 
proportional  and  Integral  gains  begins  by  setting  the  a 


and  e  gain  multipliers  to  magnitudes  of  one.  Next, 

Iterate  the  individual  weighting  matrix  elements  to  opti¬ 
mize  the  transient  response  for  each  output.  In  effect, 
this  sets  the  best  proportional  gain,  which  primarily 
affects  the  transient  response  of  the  system.  Finally, 
adjust  only  the  integral  gain  to  improve  steady  state  per- ‘ 
f ormance . 

However,  in  MULTI,  the  alpha  factor  varies  the  prov 
portional  gain  Instead  of  the  desired  integral  gain  as 
shown  below: 


Kq  "  ae[CB]“^I 

(4-32) 

-  e[CB]"^Z 

(4-33) 

where 

Kq  ■  proportional  gain  matrix 
•  Integral  gain  matrix 

To  efficiently  solve  this  problem,  first  set  the  c 
parameter  to  the  desired  integral  gain.  Next,  enter  an 
a  parameter  equal  to  the  reciprocal  of  e.  This  adjusts  the 
proportional  gain  back  to  its  original  value.  The  result 
is  t  multiplying  the  matrix  by  the  appropriate  integral 
gain. 

4-12  MULTI 

The  successful  design  of  a  digital  controller  for 
a  realistic  aircraft  plant  using  the  techniques  developed 
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by  Profssaor  Porter  requires  an  interactive  computer  pro** 
gram  to  perform  the  arithmetic  operations  necessairy  to 
calculate  the  control  matrices >  run  the  simulations/  and 
then  Iterate  the  process  many  times  to  achieve  the  optimum 
design.  The  interactive  computer  program  called  "MULTI" 
provides  that  capability  (I) . 

Kq  and  Matrices.  After  entering  the  aircraft 
plant,  the  desired  sigma  weighting  matrix,  epsilon,  alpha, 
and  measurement  matrix  for  irregular  designs,  MULTI  calcu¬ 
lates  the  Kq  andK-|^  control  matrices  using  Porter's  con¬ 
trol  law.  As  an  illustration  with  the  X-29  aircraft,  0.7 
Mach  at  15000  feet,  no  actuators,  and  no  time  delay;  the 
optimum  design  parameters  and  the  resulting  control 
matrices  calculated  using  MULTI  are  as  follows: 

a  -  0.71429 

e  -  1.4 


L- 


1.7  0  0 

0  1.2  0 

0  0  0.01 


6.652 

-9.955 


0.007181 

-0.1295 


-2,819 

-8.621 


0.02397 


0.04419 


-0.002194 


9.313 


0.01005 


-3.947 


13.94  -0.1813 

0.03356  0.6187 


-12.07 

0.003071 


Simulation.  In  addition  to  the  above  items ,  MULTI 
requires  the  sampling  time,  computational  delay  time,  total 
simulation  time,  calculation  step  size,  command  input  vec¬ 
tor,  and  actuator /sensor  state  space  equations  before 
running  the  simulation. 

This  thesis  uses  a  sample  period  of  0.025  seconds 
which  is  easily  within  the  capability  of  current  flight  con 
trol  microprocessors. 

Final  designs  in  this  research  effort  employ  a 
computational  time  delay  of  one  sample  period  to  allow  the 
digital  controller  time  to  calculate  the  control  signal. 

The  total  simulation  time  must  be  long  enough  to 
ensure  the  transient  response  is  completed  and  the  steady 
state  response  remains  stable.  This  report  uses  a  simula¬ 
tion  time  of  16  seconds. 

The  calculation  step  size  is  set  equal  to  the 
sampling  period  unless  there  is  a  special  need  to  observe 
the  response  between  sampling  times. 

Each  maneuver  is  entered  into  MULTI  via  a  command 
vector.  MULTI  has  the  capability  to  pulse  and  ramp-to- 
hold  all  input  variables. 


MULTI  has  the  versatility  to  add  actuator  and 
sensor  dynamics  external  to  the  aircraft  plant,  with  this 
Implementation,  the  actuators  and  sensors  are  not  used  In 
the  computatiO'n  of  Kq  and  control  matrices.  However, 
the  actuator  and  sensor  dynamics  are  Included  In  the  com¬ 
plete  closed- loop  system  simulation. 

In  this  paper,  the  system  Is  first  designed  without 
actuators.  Then,  using  the  same  design  parameters  as  a 
starting  point,  the  simulation  Is  rerun  with  actuator 
dynamics  added.  Based  on  an  analysis  of  the  time  response 
plots,  the  parameters  are  modified  and  the  simulation  is 
run  again.  The  process  Is  repeated  several  times  until 
the  best  possible  design  Is  achieved.  The  same  iterative 
procedure  is  employed  when  computer  time  delay  Is  added  to 
the  system. 

Appendix  C  documents  an  improvement  to  MULTI  that 
provides  the  capability  to  easily  convert  the  simulation 
output  from  units  of  radians  to  degrees  prior  to  plotting 
the  data  and  calculating  the  figures  of  merit.  At  the 
prompt,  the  user  enters  the  vectors  and  states  requiring 
conversion. 

4-13  Summary 

This  chapter  details  the  11  major  steps  necessary 
to  successfully  design  a  digital  controller  with  actuators 
and  time  delay  for  the  X-29  aircraft.  For  clarity. 


numerical  example*  illuatrate  many  of  the  step*  in  the  prO' 
cedure.  The  next  chapter  preaents  the  design  data  and 
evaluates  the  results. 


5.  ReaultB  for  Longitudinal  Control 


System  Designs 


5“1  Introduction 

This  chapter  details  the  design  data  and  presents 
time  history  response  plots  for  the  X-29  aircraft,  Imple-  . 
mented  with  a  "Porter  configured"  error  actuated  flight 
control  system.  This  thesis  Includes  results  for  the 
direct  climb,  pitch  pointing,  and  vertical  translation 
maneuvers  at  0.4  Mach,  sea  level;  0.7  Mach,  15000  feet; 
and  1.2  Mach,  15000  feet.  For  illustration  and  analysis, 
responses  are  presented  for  the  three  longitudinal  maneu¬ 
vers  with  samples  selected  from  the  three  flight  conditions. 
Design  data  and  plots  for  the  maneuvers  at  flight  condi¬ 
tions  not  detailed  In  this  chapter  are  Included  In  the 
appendices. 

This  research  effort  Investigates  the  effects  of 
adding  actuator  dynamics  and  computational  time  delay  to 
the  simulation.  For  each  maneuver  and  flight  condition,  a 
series  of  longitudinal  controllers  are  designed  with  the 
following  system  configurations! 

1.  basic  aircraft  equations  (plant)  only 

2.  actuator  dynamics  added  as  an  Input  to  the 

plant 


3.  batlc  plant  vlth  one  sample  period  computa¬ 
tional  time  delay 

4.  plant  with  both  actuator  dynamics  and  computa¬ 
tional  time  delay 

5-2  Direct  Climb  Maneuver  for 

6»)  Mach  at  15000  Feet 

The  longitudinal  state  space  matrices  are  tabu¬ 
lated  In  Table  5-1.  The  command  vector  v,  design  param¬ 
eters;  and  controller  matrices  for  this  maneuver  are  pre¬ 
sented  In  Table  5-2.  Table  5-3  lists  the  figures  of  merit. 
The  aircraft  time  response  plots  are  provided  In  Figs.  5-1 
through  5-24.  Note  that  all  angles  are  converted  from 
radians  to  degrees  before  plotting. 

In  this  maneuver  a  pitch  rate  pulse  command  gener¬ 
ates  a  steady  state  climb  angle  of  6  ■  7.5  degrees  with  a 
one  g  change  In  normal  acceleration  during  the  rotation 
phase  of  the  maneuver.  This  is  within  the  linear  range  of 
the  aircraft  equations  of  motion.  The  change  In  velocity 
and  angle  of  attack  are  commanded  to  zero. 

An  analysis  of  the  time  histories  for  all  configura¬ 
tions  Including  actuators  and  computational  time  delay  shows 
very  good  aircraft  responses. 

As  expected;  the  time  plots  with  no  actuators  and 
no  computational  time  delay  are  the  best.  For  example, 
the  aircraft  closely  and  smoothly  follows  the  commanded 
pitch  rate  pulse  with  only  a  slight  overshoot  (see  Pig.  5-1)  . 
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TABLE  5-1 


LONGITUDINAL  STATE  SPACE  MATRICES 


FLIGHT 

• 

CONDITION;  0.7  Mach  at  15000 

Feet 

A  (Plant  Matrix) 

.OOOOE+OO 

-.3213E+02 

-.3243E-06 

.OOOOE+00 

.OOOOE+00  .OOOOEtOO 

-.1199E-01  -.1718E+01 

-.1298E-03  -.1138E+01 

.2986E-03  .1892E+02 

.lOOOE+Ol 

-.3504E+00 

.9917E+00 

-.4928E-I-00 

B  (Control  Input  Matrix) 

.OOOOE+00 

-.5895E-01 

-.1161E-02 

.1720E+00 

.OOOOBtOO 
.2621B-01 
-.7808B-03 
- .  56  30B-0i 

.OOOOE-fOO 

.2724E+02 

-.2099E-02 

-.1929E+00 

C  (Output  Matrix) 

.OOOOE+00 

.OOOOE+00 

.OOOOB+00 

.OOOOE+00  .OOOOE+00 
.lOOOE+01  .OOOOE+00 
.OOOOE+00  .lOOOE+01 

.lOOOE+01 

.OOOOE+00 

.OOOOE+00 

Notes 

1.  States  (listed  in  order)  are  pitch  angle, 
change  in  forward  velocity,  angle-of-attack,  and  pitch 
rate. 


2.  Control  inputs  (listed  in  order)  are  canards, 
8tra](e  flaps,  and  thrust. 

3.  Outputs  (listed  in  order)  are  pitch  rate, 
change  in  forward  velocity,  and  angle  of  attack. 
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TABLE  5-2 


DESIGN  PARAMETERS  AND  CONTROLLER  MATRICES 
FOR  LONGITUDINAL  CONTROLLERS 


FLIGHT  CONDITION:  0.7  Mach  at  15000  Feet 

MANEUVER;  Direct  Climb 

COMMAND  VECTOR  v*:  1.1,  0.0436,  3,  4.1 


Alpha 

Epsilon 

Sigma** 

0,  6;  0, 
0,  0,  0, 

Plant  Only 

0 

0 

Kq*** 

1.7 

-.6652E+01 

.7181E-02 

-.2819E+01 

0.7143 

1.4 

1.2 

-.9956E+01 

-.1295E+00 

-.862lE-f01 

0.01 

.2397E-01 

.4419E-01 

.2194E-02 

Plant  and  Computational  Time  Delay 

Alpha 

Epsilon 

Sigma** 

V  *** 

-0 

1.9 

.7435E+01 

.7181E-02 

-. 5639E+01 

0.7143 

1.4 

1.2 

-.1113E+02 

-.1295E+00 

-.1724E+02 

0.02 

.2680E-01 

.4419E-01 

.4387E-02 

Plant  with  Actuator  Dynamics 

Alpha 

Epsilon 

Sigma** 

1.9 

.7435E+01 

.7181E-02 

-.2819E+01 

0.7143 

1.4 

1.2 

-.1113E+02 

-.1295E+00 

-.8621E+01 

0.01 

. 2680E-01 

.4419E-01 

.2194E-02 

Plant 

with  Actuator  Dynamics  and  Computational 

Time  Delay 

K  *** 

Alpha 

Epsilon 

Sigma** 

-0 

2.3 

.9000E+01 

.7181E-02 

-.4229E+01 

0.7143 

1.4 

1.2 

-.1347E+02 

-.1295E+00 

-.1293E+02 

0.015 

.3244E-01 

. 4419E-01 

. 3290E-02 

Notes 

^Each  pulse  entry  in  v  has  four  parts: 

1.  The  time  in  seconds  the  input  reaches  steady- state; 

2.  Steady-state  value  (angles  in  radians); 

3.  Time  the  input  leaves  steady-state;  and 

4.  Time  the  input  reaches  zero. 

**Diagonal  element  values  of  the  matrix. 

***Kq  =  alpha  *  K^. 
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TABLE  5-3 


LONGITUDINAL  FIGURES  OF  MERIT 

FLIGHT  CONDITION:  0.7  Mach  at  15000  Feet 

MANEUVER:  Direct  Climb 

COMMAND  VECTOR  v:*  1.1,  0.0436,  3,  4.1 

0 /  0 /  Of  0 

0 /  Of  Of  0 


Plant  Only 


Output 

Peak  Value** 

Peak  Time 

Sett line  Time 

Pitch  Rate 

2.66 

1.4 

*** 

Forward  Velocity 

-0.02 

2.8 

•khii 

Angle  of  Attack  ' 

1.22 

1.925 

*** 

Pitch  Angle 

•  7.7 

3.95 

4.0 

Flight  Path  Angle 

8.1 

5.0 

8.0 

Plant 

and  Computational  Time 

Delay 

Output 

Peak  Value** 

Peak  Time  Settling  Time 

Pitch  Rate 

2.78 

1.525 

Forward  Velocity 

-0.0414 

2.625 

*** 

Angle  of  Attack 

1.28 

1.925 

*** 

Pitch  Angle 

7.7 

3.9 

4.0 

Flight  Path  Angle 

8.2 

4.5 

8.0 

Plant  witn  Acruator  Dynamics 


Output 

Pitch  Rate 
Forward  Velocity 
Angle  of  Attack 
Pitch  Angle 
Flight  Path  Angle 


Peak  Value** 

Peak  Time 

Settling  Time 

2.71 

1.4 

*** 

-  .0199 

2.45 

*** 

1.29 

1.925 

*** 

7.6 

3.8 

7.6 

8.2 

5.0 

8.1 

Plant  with  Actuator  Dynamics  and  Computational  Time  Delay 


Output 

Pitch  Rate 
Forward  Velocity 
Angle  of  Attack 
Pitch  Angle 
Flight  Path  Angle 


ak  Value** 

Peak  Time 

Settling  Time 

2.85 

1.525 

h  hie 

-0.033 

2.0 

*  *  * 

1.4 

1.925 

*¥fk 

7.7 

4.0 

8.0 

8.2 

5.0 

8.1 

Notes 

1.  Figures  of  merit  for  the  pitch  angle  and  the 
flight  path  angle  are  estimated  from  the  response  plots. 

2.  Pitch  rate,  forward  velocity,  and  angle  of 
attack  are  components  of  the  output  vector  y. 

3.  Pitch  rate  and  angle  of  attack  are  commanded 
to  achieve  desired  responses  in  the  pitch  angle  and  the 
flight  path  angle  for  each  maneuver. 

‘Commanded  inputs  in  order: 

(1)  pitch  rate  in  radians  per  second; 

{2)  forward  velocity  in  feet  per  second;  and 
(3)  angle  of  attack  in  radians. 

“Angles  are  converted  from  radians  to  degrees  prior 
to  calculating  the  figures  of  merit. 

** ‘Final  value  of  the  output  is  zero.  Settling  time 
undefined  using  definition  of  2  percent  of  final  value. 
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Figure  5-3 


Figure  5-4 
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Figure  5-20 
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Figure  5-22 
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The  change  in  velocity  in  Fig.  5-5  is  very  slight.  The 
angle  of  attack  momentarily  peaks  at  1.25  degrees »  then 
decreases  to  -0.9  degrees  during  the  transient  before 
returning  to  zero  degrees  in  the  steady  state  (see  Fig. 

5-9) .  The  maximum  strake  flap  deflection  of  30  degrees  in 
Fig.  5-13  occurs  at  approximately  3.5  seconds.  The  strake 
flap  deflection  rate  peaks  at  16  degrees  per  second  which 
is  well  within  the  rate  limit  of  27  degrees  per  second. 
Because  the  deflection  and  rate  limits  of  the  strake  flap 
are  the  governing  factors  in  all  the  maneuvers »  the  canard 
deflection  position  and  rate  always  remain  within  limi,ts 
(see  Figs.  5-13  through  5-16).  An  analysis  of  the  time 
history  of  the  control  surface  movements  shows  the  canard 
initially  deflects  in  a  positive  direction  to  rotate  the 
nose  of  the  aircraft  upwards.  Because  the  aircraft  is 
statically  unstable^  the  aircraft  tends  to  continue 
rotating.  The  strake  flap  offsets  the  rotation  by  pro¬ 
ducing  a  negative  moment.  This  is  observed  by  the  positive 
deflection  and  resulting  positive  lift  generated  by  the 
strake  flap  which  is  located  aft  of  the  center  of  gravity. 
The  maneuver  is  effj.ciently  accomplished  since  the  lift 
created  by  both  the  canard  and  strake  flap  during  the 
maneuver  remains  positive  at  all  times.  In  contrast,  a 
pitch  up  maneuver  in  a  conventional  aircraft  requires  the 
stabilator  or  elevator  to  produce  a  negative  lift,  thereby 
reducing  the  total  lift  on  the  vehicle. 
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the  design  parameters  completely  eliminate  the  oscilla¬ 
tions. 

Responses  for  the  direct  climb  maneuver  are  also 
very  good  for  the  other  flight  conditions.  The  design 
parameters  and  plots  are  presented  in  Appendices  D  through 
•F. 


5-3  Vertical  Translation  Maneuver 

i^or  0.4'‘Mach  at  Sea  Level 

The  aircraft  longitudinal  state  space  matrices  are 
tabulated  in  Table  5-4.  The  command  vector  V/  design 
parameters,  and  controller  matrices  for  this  maneuver  are 
presented  in  Table  5-5.  The  figures  of  merit  are  listed 
in  Table  5-6.  The  aircraft  time  plots  are  provided  in 
Figs.  5-25  through  5-48.  Note  that  the  angles  in  all 
plots  are  in  degrees. 

The  objective  of  this  maneuver  is  to  change  the 
aircraft's  altitude  without  rotating  the  vehicle  while 
maintaining  a  constant  forward  velocity.  The  maneuver 
requires  that  pitch  rate  and  forward  velocity  be  commanded 
to  zero.  A  commanded  negative  angle  of  attack  increases 
the  aircraft's  altitude  and  a  positive  angle  of  attack 
decreases  the  vehicle's  altitude.  Thus,  the  aircraft  trans¬ 
lates  upwards  only  when  the  canard  and  strake  generate  a 
change  ir  positive  lift  that  is  greater  than  the  li.ft  lost 
on  the  main  wings. 


TABLE  5-4 


LONGITUDINAL  STATE  SPACE  MATRICES 


FLIGHT  CONDITION:  0.4  Mach  at  Sea  Level 


(Plant  Matrix) 


.OOOOE+00 

-.3217E+02 

.OOOOEtOO 

.OOOOE+00 


.OOOOE+00 

-.1660E-01 

-.3172E-03 

.2823E-03 


.OOOOE+00 
. 4880E-01 
-.lOOlE+01 
.1012E+02 


.lOOOE+01 

-.4273E+00 

.9880E+00 

-.4340E+00 


(Control  Input  Matrix) 


.OOOOE+00  .OOOOE+00  .OOOOE+00 

-.4891E-01  .2592E-01  .3779E+02 

-.9426E-03  -.8021E-03  -.6761E-02 

.9328E-01  -.3603E-01  -.2679E+00 


(Output  Matrix) 


.OOOOE+00 
.  OOOOE+00 
.  OOOOE+00 


.OOOOE+00 

.lOOOE+01 

.OOOOE+00 


.OOOOE+00 

.OOOOE+00 

.lOOOE+01 


.lOOOE+01 
.OOOOE+00 
.  OOOOE+00 


Notes 

1.  States  (listed  in  order)  are  pitch  angle « 
change  in  forward  velocity,  angle-of-attack,  and  pitch 
rate . 

2.  Control  inputs  (listed  in  order)  are  canards, 
St rake  flaps,  and  thrust. 

3.  Outputs  (listed  in  order)  are  pitch  rate, 
change  in  forward  velocity,  and  angle  of  attack. 


TABLE  5-5 


DESIGN  PARAMETERS  AND  CONTROLLER  MATRICES 
FOR  LONGITUDINAL  CONTROLLERS 


■  FLIGHT  CONDITION:  0.4  Mach  at  Sea  Level 
MANEUVER;  Vertical  Translation 


COMMAND  VECTOR 

V*;  0 ,  0,  0 

,  0 

0,  0,  0 

/  0 

1.1,  -0 

.0349,  50,  1 

50 

Plant  Only 

Alpha 

Epsilon 

‘Sigma** 

V  *** 

-0 

1.5 

.1105E+02 

-.'7022E-02 

-.3313E+01 

0.4545 

2.2 

1.0 

-.1319E+02 

-.2160E+00 

-.8588E+01 

0.01 

.2335E-01 

.2660E-01 

.1603E-02 

Plant  and  Computational 

Time  Delay 

Alpha 

Epsilon 

Sigma** 

V  At*  * 

-0 

1.5 

.1105E+02 

-. 7022E-02 

-.8281E+01 

0.4545 

2.2 

1.0 

-.1319E+02 

-.2160E+00 

-.2147E+02 

0.025 

.2335E-01 

.2660E-01 

.4008E-02 

Plant  with  Actuator  Dynamics 

Alpha 

Epsilon 

Sigma** 

Kq*** 

1.  5 

.1105E+02 

-. 7022E-02 

-.2484E+01 

0.4545 

2.2 

1.0 

-.1319E+02 

-. 2160E+00 

-.6441E+01 

0.0075 

.2335E-01 

. 2660E-01 

.1202E-02 

Plant 

with  Actuator  Dynamics  and  Computational 

Time  Delay 

Alpha 

Epsilon 

Sigma** 

V  *** 

-0 

1.5 

.1105E+02 

-.1053E-01 

-.6625E+01 

0.4545 

2.2 

1.5 

-.1319E+02 

-. 3240E+00 

-.1718E+02 

0.02 

.2335E-01 

. 3990E-01 

.  3206E-02 

Notes 

*Each  pulse  entry  in  v  has  four  parts; 

1.  The  time  in  seconds  the  input  reaches  steady- state; 

2.  Steady-state  value  (angles  in  radians); 

3.  Time  the  input  leaves  steady- state;  and 

4.  Time  the  input  reaches  zero. 

**Diagonal  element  values  of  the  matrix. 

***Kq  =  alpha  *  K^. 
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TABLE  5-6 


LONGITUDINAL  FIGURES  OF  MERIT 


FLIGHT  CONDITION:  0.4  Mach  at  Sea  Level 
MANEUVER:  Vertical  Translation 


TABLE  5-6  — Continued 


Plant  with  Actuator  Dynamics 


Output 

Peak  Value** 

Peak  Time 

Settling  Time 

Pitch  Rate 

-0.0772 

1.925 

Forward  Velocity 

-0.0009 

0.7 

*** 

Angle  of  Attack 

-2.03 

5.95 

4.375 

Pitch  Angle 

-0.17 

4.1 

4.1 

Plight  Path  Angle 

1.85 

5.0 

5.0 

Plant  with  Actuator  Dynamics  and  Computational  Time  Delay 


Output 

Peak  Value** 

Peak  Time 

Settling  Time 

Pitch  Rate 

-0.22 

1.75 

*  A* 

Forward  Velocity 

-0.0019 

0.525 

A** 

Angle  of  Attack 

-2.05 

4.2 

5.25 

Pitch  Angle 

-0.3 

4.0 

4.1 

Flight  Path  Angle 

1.8 

3.9 

4.5 

- 

Notes 


1.  Figures  of  merit  for  the  pitch  angle  and  the 
flight  path  angle  are  estimated  from  the  response  plots. 

2.  Pitch  rate,  forward  velocity,  and  angle  of 
attack  are  components  of  the  output  vector 

3.  Pitch  rate  and  angle  of  attack  are  commanded 
to  achieve  desired  responses  in  the  pitch  angle  and  the 
flight  path  angle  for  eech  maneuver. 

*Conunanded  inputs  in  order: 

(1)  pitch  rate  in  radians  per  second; 

(2)  forward  velocity  in  feet  per  second;  and 

(3)  angle  of  attack  in  radians. 

**Angles  are  converted  from  radians  to  degrees  prior 
to  calculating  the  figures  of  merit. 

***Final  value  of  the  output  is  zero.  Settling  time 
•'ndefined  using  definition  of  2  percent  of  final  value. 
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Figure  5-28 


VERTICAL  TRANSLATION  (C.4M  OK  FT' 

Figure  5-33 


CRNARO  RNO  STRRKE  FLRPS 


VERTtCRL  TRANSLATION  (0.4M  OK  FT) 


Figure  5-37 


VERTICAL  TRANSLAMON  10.411  OK  FT)  DELAf=0.C25  SEC 
Figure  5-38 


VERTir.qi  TRANSLATION  (O.^M  OK  FT 


Figure  5-41 


VERTICAL  TRANSLATION  (O.dM  OK  FT)  DELAr=0,C2S  SEC 

Figure  5-42 


12.00 


14.00 


te.oo 


In  this  manauver,  a  ramp-to-hold  Input  signal  com¬ 
mands  the  angle  of  attack  to  negative  two  degrees.  In  addi¬ 
tion,  the  pitch  rate  and  velocity  are  commanded  to  zero. 

The  result  Is  an  upward  translation  of  the  aircraft. 

All  the  time  plots  show  smooth  and  well  behaved 
aircraft  responses  with  and  without  actuators  and  computa¬ 
tional  time  delay. 

For  the  configuration  without  actuators  and  computa¬ 
tional  time  delay,  the  aircraft  smoothly  follows  the  com¬ 
manded  angle  of  attack  with  just  a  slight  overshoot  as 
shovm  In  Fig.  5-33.  The  velocity  changes  only  a  fraction 
of  a  foot  per  second  during  the  maneuver  (see  Fig.  5-29) . 

The  aircraft  rotates  slightly  during  the  transient  response 
and  correspondingly,  the  pitch  rate  of  Fig.  5-25  momentarily 
peaks  at  minus  0.06  degrees  per  second  before  returning  to 
zero.  Because  it  la  unrealistic  to  expect  the  pitch  rate 
to  remain  at  exactly  zero  during  the  maneuver,  a  small 
pitch  angle  error  results  from  the  Integral  of  the  pitch 
rate  response  and  the  resulting  pitch  angle  Is  about  minus 
0.2  degrees  (see  Fig.  5-45).  Therefore,  the  flight  path 
angle  la  slightly  less  than  the  commanded  angle  of  attack 
(see  Figs.  5-45  through  5-48).  This  occurs  because  the 
flight  path  angle  Is  Lhe  dl^rerence  oetween  the  pitch  angle 
and  the  angle  of  attack. 

When  rates  are  the  input  variables,  then  angles  are 
indirectly  controlled.  Therefore,  a  slight  error  always 


occurs  in  the  steady  state  response.  This  .error  Is  com¬ 
pletely  eliminated  only  i£  angles  are  the  Input  variables. 

The  canard  and  strake  deflections  are  both  positive 
and  within  the  deflection  limits.  The  waveshapes  are 
smooth  with  a  slight  overshoot  (Fig.  5-37) . 

Adding  computational  time  delay  requires  an 
Increase  In  the  angle  of  attack  weighting  element  In  the 
sigma  matrix.  With  the  additional  gain,  the  results  are 
very  similar  to  the  responses  for  the  basic  aircraft  con¬ 
figuration. 

When  actuator  dynamics  are  Incorporate.^,  as  an  Input 
to  the  aircraft  plant,  the  sigma  weighting  matrix  element 
controlling  the  angle  of  attack  Is  reduced  slightly  to 
minimize  the  overshoot.  As  a  result.  It  takes  the  aircraft 
approximately  0.5  seconds  longer  to  reach  the  commanded 
angle  of  attack  (see  Fig.  5-35).  The  pitch  rate  goes 
momentarily  positive,  then  negative,  making  the  area  under 
the  pitch  rate  curve  small  (see  Fig,  5-27).  The  resulting 
pitch  angle  err^or  of  0.15  degrees  In  Fig.  5-47  Is  less  than 
for  the  basic  aircraft.  The  canard  and  strake  deflections 
In  Fig.  5-39  are  within  limits  and  have  a  smooth  response. 

The  design  of  the  controller  with  both  actuator  and 
dynamics  and  computational  time  delay  has  very  good  time 
responses  which  are  very  similar  to  those  of  the  basic 
aircraft.  To  achieve  the  best  design,  the  weighting  matrix 
element  corresponding  to  the  angle  of  attack  Is  Increased 
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to  a  value  lower  than  for  the  configuration  with  only  time 
delays  but  higher  than  for  the  design  with  just  actuators. 

Responses  for  the  vertical  translation  maneuver 
are  also  very  good  for  the  other  flight  conditions.  The 
design  parameters  and  plots  are  presented  In  Appendices  D 
through  F. 


5-4  Pitch  Pointing  Maneuver  for 
0.4  Mach  at  Sea  Level 

The  longitudinal  state  space  matrices  are  listed 
in  Table  5-4.  The  command  vector  v,  design  parameters, 
and  controller  matrioes  for  this  maneuver  are  tabulated  in 
Table  5-7.  Table  5-8  presents  the  figures  of  merit.  The 

aircraft  time  response  plots  are  provided  in  Figs.  5-49 

» 

through  5-72.  Note  that  all  angles  are  converted  from 
radians  to  degrees  before  plotting. 

The  purpose  of  this  maneuver  is  to  change  the  air¬ 
craft  pitch  attitude,  i.e.,  point  the  nose  of  the  aircraft 
longl>.udinally ,  without  changing  the  vehicle's  flight  path 
a,nglQ  or  velocily<  To  accomplish  this  maneuver,  the  pitch 


rate  oOiiunami  is  pulsed  to  gene-.-te  a  desired  pitch  angle 
and  the  chang*'^  of  velocity  is  commanded  to  zero.  In  addi¬ 
tion,  the  angle  of  attack  J.s  commanded  to  equal  the  nevr 
pitch  angle  so  the  resulting  flight  path  angle  remair>s 
unchanged . 
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TABLE  5-7 


DESIGN  PARAMETERS  AND  CONTROLLER  MATRICES 
FOR  LONGITUDINAL  CONTROLLERS 


PLIGHT  CONDITION:  0.4  Mach  at  Sea  Level 

MANEUVER;  Pitch  Pointing 

COMMAND  VECTOR  v*s  0.8,  0.0349,  1.0,  1.8 

0,  0,  0,  0 

0.8,  0.0349,  50,  50 


Plant  Only 


Alpha 

Epsilon 

Sigma** 

Kq*** 

0.357 

2.8 

1.9  .1400E+02 
1.0  -.1670E+02 
0.015  .2957E-01 

-.7019E-02 

-.2159E+00 

.2659E-01 

-.4967E+01 

-.laSSEtOO 

.2404E-02 

Plant  and  Computational  Time  Delay 

Alpha 

Epsilon 

Sigma** 

jjj^lQQQIIIIIIIIIII^ 

0.357 

2.8 

1.9  .1400E-f02 
1.0  -.1670E+02 
0,015  .2957B-01 

-.7019E-02 

-.2159E+00 

.2659E-01 

-.4967E+01 

-.1288E4-02 

.2404E-C2 

Plant  with  Actuator  Dynamics 

Alpha 

Epsilon 

Sigma** 

V  *** 

io 

0.357 

2.8 

1.9  .1400Et02 
1.0  -.1670E+02 
0.105  .2957E-01 

-.7019E-02 

-.2159E+00 

.2659E-01 

-.4967E+01 

-.1266E-l‘02 

.2404E-02 

Plant 

with  Actuator  Dynamics  and  Computational 

Time  Delay 

Alpha 

Epsilon 

Sigma** 

Kq*** 

0.357 

2.8 

1.6  .1179E+02 
0.5  -.1406E+02 
0.015  .2490E-01 

-. 3510E-02 
-.1079E+00 
.1330E-01 

-.4967E+01 

-.1288Et02 

.2404E-02 

Notes 

*'Each  pulse  entry  in  v  has  four  parts: 

1.  The  time  in  seconHs  the  input  reaches  steady- state; 

2.  Steady-state  value  (angles  in  radians); 

3.  Time  the  input  leaves  steady- state;  and 

4.  Tinu  the  input  reaches  zero. 

**Diagonc.i  v' ^ment  values  of  the  matrix. 

***^0  “  alpha  *  Kj . 


97 


TABLE  5-8 


LONGITUDINAL  FIGURES  OF  MERIT 

FLIGHT  CONDITION!  0*4  Mach  at  Sea  Level 

MANEUVER!  Pitch  Pointing 

COMMAND  VECTOR  V!*  0.8,  0.0349,  1.0,  1.8 

■“  0,  0,  Of  0 

0.8,  0.0349,  50,  50 


Plant  Onl^ 


Output 


Peak  Value**  Peak  Time  Settling  Time 


Pitch  Rate 
Forward  Velocity 
Angle  of  Attack 
Pitch  Angle 
Flight  Path  Angle 


2.08 

0.875 

*** 

0.008 

1.225 

*** 

2.02 

5.425 

4.2 

2.05 

1.8 

1.8 

0.3 

2.5 

4.5 

Plant  and  Computational  Time  Delay 


Output 


Peak  Value**  Peak  Time  Settling  Time 


Pitch  Rate 
Forward  Velocity 
Angle  of  Attack 
Pitch  Angle 
Flight  Path  Angle 


2.19 

1.05 

il*i* 

0.018 

1.225 

*  *  * 

2.01 

7.175 

5.6 

2.05 

1.8 

1.9 

0.7 

2.5 

6.0 
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TABLE  5-8  — Continued 


Plant  with  Actuator  Dynamics 

Output 

Peak  Value** 

Peak  Time 

Settling  Time 

Pitch  Rate 

2.16 

1.05 

*ilt* 

Forward  velocity 

-0.0102 

1.05 

*** 

Angle  of  Attack 

2.036 

5.25 

4.2 

Pitch  Angle 

2.1 

1.6 

1.6 

Flight  Path  Angle 

0.3 

2.4 

5.6 

Plant  with  Actuator  Dynamics  and  Computational  Time  Delay 


Output 

Peak  Value* 

*  Peak  Time 

Settling  Time 

Pitch  Rate 

2.26 

0.875 

*** 

Forward  Velocity 

-0.045 

1.225 

*** 

Angle  of  Attack 

2.02 

6.65 

5.25 

Pitch  Angle 

2.01 

1.8 

2.5 

Flight  Path  Angle 

0.75 

2.3 

5.2 

Note  a 


‘  1,  Figures  of  merit  for  the  pitch  angle  and  the 

flight  path  angle  are  estimated  from  the  response  plots. 

2.  Pitch  rate>  forvard  velocity,  and  angle  of 
attack  are  components  of  the  output  vector 

3.  Pitch  rate  and  angle  of  attack  are  commanded 
to  achieve  desired  responses  in  the  pitch  angle  and  the 
flight  path  angle  for  each  maneuver. 

^Commanded  Inputs  In  order: 

(1)  pitch  rate  In  radians  per  second; 

(2)  forward  velocity  In  feet  per  second;  and 

(3)  angle  of  attack  In  radians. 

** Angles  are  converted  from  radians  to  degrees  prior 
to  calculating  the  figures  of  merit. 

***Final  value  of  the  output  is  zero.  Settling  time 
undefined  using  definition  of  2  percent  of  final  value. 
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PITCH  PCINTINO  <0.4f1  OK  F^)  OE;Lflr=0.025  '5EC  flC’-UqTORS 


Figure 


5-52 
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^ORt•RRO  VElOCITr  (FT/SEC) 
.'2  -0.C8  -O.C4  C.CQ 


Figure  5-58 


Figure  5-6A 
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An  analysis  of  the  time  history  plots  for  all  con¬ 
figurations  of  actuators  and  oomputational  time  delay 
show  smooth  and  well  behaved  responses. 

Without  actuator  dynamics  and  computational  time 
delay  in  the  aimulationi  Fig.  5-49  shows  the  pitch  rate 
response  of  the  aircraft  slightly  overshoots  the  pitch  rate 
pulse  input  command.  The  overshoot  is  greater  when  time 
delay  and/or  actuator  dynamics  are  added  to  the  simulation 
(see  Figs.  5-50  through  5-52).  As  a  result,  the  actual 
pitch  angle  response  has  a  steady  state  value  slightly 
greater  than  calculated  from  the  pitch  rate  pulse  input 
signal.  However,  the  angle  of  attack  control  input  is  set 
equal  to  the  desired  (not  actual)  pitch  angle  because  it 
is  not  possible  to  exactly  precompute  the  actual  pitch 
angle  response.  The  angle  of  attack  is  slightly  less  than 
the  resulting  pitch  angle  depending  on  the  overshoot  in 
the  actual  pitch  rate  response.  The  difference  between 
the  pitch  angle  and  angle  of  attack  determines  the  error  in 
the  flight  path  angle.  The  flight  path  angle  error  ranges 
from  0.1  degrees  for  the  basic  aircraft  configuration  in 
Fig.  5-69  to  0.3  degrees  with  oomputational  time  delay  and 
actuator  dynamics  added  to  the  simulation  as  illustrated 
in  Fig.  5-72. 

A  maximum  of  two  degrees  of  pitch  pointing  are 
possible  without  exceeding  the  canard  and  strake  deflection 
limits  for  0.5  Mach  at  sea  level.  The  limits  vary, 
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depending  on  the  flight  condition,  and  are  easily  calcu¬ 
lated  using  the  0(0)"^  steady  state  transfer  function 
matrix  as  described  In  Chapter  4. 

Analysis  of  the  canard  and  strake  flap  deflection 
positions  in  Fig.  5-61  through  5-64  show  the  canard  and 
strake  correctly  generate  a  positive  lift  to  offset  the 
additional  positive  lift  produced  by  the  increased  angle 
of  attack  on  the  main  wings. 

The  velocity  control  for  all  oases  is  excellent 
with  a  drop  of  less  than  0.1  feet  per  second  during  the 
transient  response  (see  Figs.  S-53  through  5-56) .  The 
thrust  increases  to  offset  the  increased  drag  induced  by 
the  higher  angle  of  attack  as  shown  in  Figs.  5-65  through 
5-68. 

The  pitch  angles  in  Figs.  5-69  through  5-72  increase 
smoothly  to  a  steady  state  value  with  a  small  overshoot. 

The  steady  state  value  of  the  pitch  angle  equals  the 
integral  of  the  pitch  rate  response. 

The  responses  for  the  pitch  pointing  maneuver  are 
also  good  for  the  other  flight  conditions.  The  design 
parameters  and  plots  are  presented  in  Appendices  D  through 
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Results  for  the  Lateral  Control  System  Designs 

6-1  Introduction 

This  chapter  presents  the  design  data  and  time 
history  response  plots  for  the  X-29  aircraft,  Implemented 
with  a  "Porter  configured"  error  actuated  flight  control 
system  for  the  lateral-directional  dynamics.  Results 
Include  control  laws  for  the  coordinated  turn  and  beta 
(sideslip)  pointing  for  five  flight  conditions  as  follows: 

1.  0.4  Mach  at  Sea  Level 

2.  0.9  Mach  at  Sea  Level 

3.  0.7  Mach  at  15,000  Feet 

4.  1.2  Mach  at  15,000  Feet 

5.  0.9  Mach  at  50,000  Feet 

For  illustration  and  analysis,  results  are  pre¬ 
sented  for  the  two  lateral  maneuvers  at  the  flight  condi¬ 
tion  of  0.7  Mach  at  15,000  feet.  Design  data  and  plots 
for  the  maneuvers  not  detailed  in  this  chapter  are  included 
in  Appendices  G  and  H. 

Because  the  lateral  control  system  is  designed  by 
primarily  following  the  design  procedure  developed  in 
Chapter  4  for  the  longitudinal  controller,  a  separate 
chapter  on  the  lateral  design  procedure  Is  not  necessary. 
However,  before  presenting  the  results  of  the  lateral 
control  system  designs,  this  chapter  first  highlights- 


significant  information  related  to  major  steps  in  the 
design  procedure  and  discusses  those  differences  between 
the  lateral  and  longitudinal  design  procedure. 

6-2  Analysis  of  Aircraft  Model 


when  actuator /sensor  dynamics  and  computational  time  delay 
are  added  to  the  basic  equations.  These  additions  may 
reduce  the  phase  margin  and  make  It  a  harder  task  to  assure 
stabilization  of  the  system.  Although  the  lateral  equa¬ 
tions  are  statically  unstable,  they  are  "more  stable"  than 
the  longitudinal  equations  which  have  poles  and  zeros 
farther  to  the  right  In  the  right  half  s-plane  (see  Chap¬ 
ter  4-3).  As  a  result.  It  Is  less  tedious  to  design  a 
controller  for  the  lateral  mode  of  the  aircraft. 

6-3  Controlled  Variables  and 
Aircraft  Controls 

For  this  research  effort,  the  bank  angle  and  the 
sideslip  angle  are  selected  as  the  controlled  variables 
which  are  represented  by  the  output  vector  y.  The  rudder 
and  flaperons  are  the  aircraft  control  surfaces  chosen  for 
the  control  input  vector  u. 

An  attempt  to  augment  the  aircraft  state  equations 
with  actuator  dynamics  and  then  control  lateral  accelera¬ 
tion  and  roll  rate  was  unsuccessful  due  to  current  limita¬ 
tions  with  the  design  method. 

6-4  Controllability,  Observability , 
and  Transmission  Zeros 

The  computer  program  "Zero"  (7)  is  used  to  calcu¬ 
late  the  transmission,  Inpu'c  decoupling,  and  output  decou¬ 
pling  zeros  of  the  state  and  output  aircraft  plant  equations. 
With  the  lateral  control  system  configured  as  described  in 
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the  previous  section#  the  aircraft  plant  is  completely 
controllable  and  observable;  i.e. #  there  are  no  input 
decoupling  or  output  decoupling  zeros#  for  all  the  flight 
conditions.  In  addition#  the  transmission  zeros  are 
located  in  the  left  half  s-plane.  For  example#  at  the 
flight  condition  of  0.7  Mach  at  1S#000  feet  the  transmis¬ 
sion  zeros  occur  at  -3.3431  and  -3.3047  with  a  measurement 
matrix  of: 


M 


0.3  0 

0  0.3 


6-5  Measurement  Matrix 

A  measurement  matrix  is  required  for  the  lateral 
control  system  because  ^  does  not  have  full  rank.  As 
explained  in  Chapter  4-7#  the  measurement  matrix  intro¬ 
duces  additional  transmission  zeros.  They  can  be  changed 
by  altering  the  measurement  matrix.  Since  the  measurement 
matrix  feeds  back  the  derivatives  of  the  states,  the 
transient  error  between  the  control  system  input  vector 
V  and  the  output  vector  y  is  minimized  by  reducing  the 
number  of  measurement  matrix  elements  to  the  smallest 
number  that  provides  a  satisfactory  response.  Also#  if  the 
input  vector  is  a  constant  value  in  the  steady  state#  then 
the  measurement  matrix  does  not  affect  the  steady  state 
error  since  the  derivatives  of  the  states  are  equal  to 


For  the  lateral  control  system  configuration  pre¬ 
sented  In  this  thesis,  a  measurement  matrix  with  diagonal 
elements  of  0.1  results  In  satisfactory  performance  of 
the  control  system  only  if  the  basic  aircraft  equations 
(without  actuator  dyneunlcs  and  computational  time  delay) 
are  Included  In  the  simulation.  However,  when  actuator 
dynamics  and/or  computational  time  delay  are  added  to  the 
closed-loop  simulation  of  control  systems,  it  is  necessary 
to  Increase  the  diagonal  measurement  matrix  elements  to  a 
value  of  0.3  to  prevent  unbounded  oscillations  in  the 
output . 

6-6  Control  Input  Magnitudes 

In  this  thesis,  the  maximum  magnitude  and  duration 
of  the  Input  command  vector  are  adjusted  to  maintain  air¬ 
craft  movement  within  the  linear  region  of  the  aircraft 
equations  of  motion.  This  makes  anlaysls  of  the  results 
easier  and  more  realistic.  For  example,  calculations  are 
made  to  determine  the  turn  rate  generated  as  a  function  of 
the  bank  angle  for  a  level  coordinated  turn  at  a  constant 
airspeed.  This  calculation  Is  Independent  of  the  aircraft 
and  depends  on  the  bank  angle  of  the  vehicle.  This  equa¬ 
tion  relating  the  airspeed  and  bank  angle  to  the  yaw  rate 
Is  derived  using  the  following  relationships.  First,  the 
lateral  g's  generated  by  an  aircraft  in  a  level  coordinated 
turn  is  equal  to  the  tangent  of  the  bank  angle.  Using  this 


relationship  and  the  equation  for  centrifugal  accelera¬ 
tion; 

a  -  v^/R  (6-6) 


where 

2 

a  •  acceleration  In  feet/second 

V  ■  tangential  velocity  In  feet  per  second 
R  ■■  turn  radius  in  feet 

the  circumference  of  a  circle: 

C  -  27rR  (6-7) 

where 

C  ■  circumference  of  circle  In  feet 
R  ">  turn  radius  of  the  aircraft  in  feet 

and  the  constants 

2 

1  g  ■  32.2  feet /second 
360  degrees  In  a  circle 

the  following  relationship  is  derived  relating  bank  angle 
and  velocity  to  yaw  rate: 

r  -  (360)  (32.2)  tan  (J>/27tv  '  (6-8) 

where 

r  =  yaw  rate  in  degrees /second 
(|i  =>  bank  angle  in  degrees 

V  “  tangential  velocity  in  feet /second 


119 


As  an  example »  for  the  flight  condition  of  0.7 
Mach  (740  feet/second)  at  15,000  feet,  a  15  degree  bank 
results  in  a  yaw  rate  of  0.668  degrees/second  using  the 
equation  above.  This  agrees  closely  with  the  results  of 
0.652  obtained  from  the  simulation  on  MULTI.  The  error 
Increases  with  greater  bank  angle  because  the  linear  equa¬ 
tions  of  motion  approximate  the  tangent  of  the  bank  angle 
as  equal  to  the  bank  angle  in  radians.  This  is  a  good 
approximation  only  for  small  angles. 

For  the  beta  pointing  maneuver  the  sideslip  angle 
is  limited  to  10  degrees  due  to  nonlinearities  in  the  drag. 

6-7  Results 

The  introduction  to  this  chapter  describes  the 
maneuvers  and  flight  conditions  used  in  the  design  of  the 
lateral  controllers  for  this  research  effort.  The  follow¬ 
ing  sections  present  the  results  for  the  coordinated  turn 
and  the  beta  pointing  maneuver  for  an  aircraft  with 
actuator  dynamics  and  computational  time  delay  added  to  the 
simulation  for  the  flight  condition  of  0.7  Mach  at  15,000 
feet.  Data  and  response  plots  for  the  remaining  flight 
conditions  are  in  Appendices  G  and  H. 

6-8  Coordinated  Turn  for  0.7  Mach 

at  15,000  Feet 

The  lateral  state  space  matrices  for  this  flight 
condition  are  tabulated  in  Table  6-1.  The  command  vector  v, 


TABLE  6-1 


LATERAL  STATE  SPACE  MATRICES 


FLIGHT 

CONDITION 1 

0.7  Mach  at  15000 

Feet 

A  (Plant  Matrix) 

.OOOOE4-00 

.OOOOEtOO 

.lOOOE+01 

.5710E-01 

.4334E-01 

-.2210E+00 

.5710E-01 

-.9984E4-00 

.OOOOEtOO 

-.1793E+02 

-.2989Et01 

.1090E-I-01 

.OOOOEtOO 

.7836EtQl 

-.9651E-01 

-. 9504E-01 

B  (Control  Input  Matrix) 

.QOOOEtOO 

.OOOOBtOO 

-.8361E-03 

.9071E-03 

.1493E-t-01 

.2670Et00 

.9608E-01 

-.8655B-01 

C  (Output  Matrix) 

.lOOOE+01 

.OOOOE+OO 

.OOOOE-t^OO 

.OOOOEtOO 

.OOOOE+OO 

.lOOOEtOl 

,00'00E+00 

.OOOOEtOO 

Notes 


1.  States  (listed  In  order)  are  bank  angle,  side¬ 
slip  angle,  roll  rate,  and  yaw  rate. 

2.  Control  inputs  (listed  in  order)  are  flaperon 
and  rudder. 

3.  Outputs  (listed  in  order)  are  bank  angle  and 
sideslip  angle. 
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design  parameters,  and  controller  matrices  for  this  maneuver 
are  presented  In  Table  6-2.  Table  6-3  lists  the  figures  of 
merit.  The  aircraft  time  response  plots  are  presented  In 
Figs.  6-1  through  6-4.  Note  that  all  angles  are  converted 
from  radians  to  degrees  before  plotting. 

In  this  maneuver  the  bank  angle  command  Is  ramped 
In  one  second  to  a  steady  state  value  of  15  degrees  (see 
Fig,  6-1) .  The  aircraft  response  follows  the  Input  command 
very  smoothly  with  a  negligible  overshoot  of  only  0.2 
degrees.  The  sideslip  angle  la  commanded  to  zero.  The 
time  response  plot  (see  Fig.  6-2)  shows  no  change  In  the 
sideslip  angle.  The  figures  of  merit  In  Table  6-3  show 
that  the  sideslip  angle  actually  peaks  at  0.5  seconds  to 

I 

a  very  small  -0.01  degrees  during  the  transient  portion  of 
the  maneuver  (too  small  to  see  on  the  plot)  before  returning 
to  zero  In  the  steady  state.  Figure  6-3  shows  a  peak 
flaperon  deflection  of  0.8  degrees  at  approximately  0.5 
seconds  and  a  peak  rudder  deflection  of  0.1  degrees  at 
0.9  seconds.  The  plots  for  the  flaperon  and  rudder  are 
relatively  smooth  and  easily  within  the  deflection  and  rate 
limits  for  both  control  surfaces.  The  roll  and  yaw  rates 
shown  In  Fig.  6-4  are  smooth  and  wall  controlled.  The  roll 
rate  peaks  at  15  degrees  per  second  at  0.5  seconds  before' 
returning  to  zero.  The  yaw  rate  peaks  at  0.7  degrees/ 
second  before  settling  at  0.652  degrees/second  which  is 
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TABLE  6-2 


DESIGN  PARAMETERS  AND  CONTROLLER  MATRICES 
FOR  LATERAL  CONTROLLERS 

MANEUVERi  Coordinated  Turn 

COMMAND  VECTOR  v*J  1,  0.2618,  50,  50 

“  0,  0,  0,  0 

ALPHAt  1.0  EPSILON:  1.0  SIGMA**:  0.5,  0.5 

MEASUREMENT  MATRIX**:  0.3,  0.3 

PLANT  WITH  ACTUATOR  DYNAMICS  AND 
COMPUTATIONAL  TIME  DELAY 


Flight  Condition  Kq  ■ 


0.4 

Mach 

at 

Sea  Level 

/  .1779E+01 

-.4041E+01 

-.5172E+00 

.2705E+02 

0.9 

Maoh 

at 

Sea  Level 

.5671E+00 

-.1557E+01 

.1405Et00 

•.7334E+01 

0.7 

Mach 

at 

15000  Feet 

.1088Et01 

-.2729E+01 

.1383E-(-00 

.1561E+02 

1.2 

Mach 

at 

15000  Feet 

.1012E+01 

-.2528E+01 

.1971E+01 

.1247E+02 

0.9 

Mach 

at 

50000  Feet 

.4174E+01 

-.1032Et02 

-.1387E+01 

.4071E+02 

Notes 

*Each  pulse  entry  in  v  has  four  parts: 

1.  The  time  In  seconds  the  Input  reaches  steady  state; 

2.  Steady  state  value  (angle  in  radians); 

3.  Time  the  input  leaves  steady  state;  and 

4.  Time  input  reaches  zero. 

**Diagonal  element  values  of  the  matrix. 


TABLE  6-3 


LATERAL  FIGURES  OP  MZ^IT 


MANEUVER:  Coordinated  Turn 

COMMAND  VECTOR  V*:  1,  0.2618,  50, 

”  0,  0,  0,  0 

PLANT  WITH  ACTUATOR  DYNAMICS  AND 
COMPUTATIONAL  TIME  DELAY 


50 


Flight  Condition:  0.4  Mach  at  Sea  Level 


Output 

Peak  Value** 

Peak  Time 

Settling  Time 

Bank  Angle 

15.2 

2.1 

1.5 

Sideslip  Angle 

-0.01 

0.5 

*** 

Roll  Rate 

15.0 

0.5 

*** 

Yaw  Rate 

2.0 

1.0 

1.5 

Flight  Condition:  0.9 

Mach  at  Sea 

Level 

Output 

Peak  Value** 

Peak  Time 

Settling  Time 

Bank  Angle 

16.08 

2.3 

3.9 

Sideslip  Angle 

-.02 

1.3 

*** 

Roll  Rate 

15.0 

1.0 

★  *  * 

Yaw  Rate 

1.0 

1.0 

1.8 

Flight  Condition:  0.7 

Mach  at  15000 

Feet 

Output 

Peak  Value** 

Peak  Time 

Settling  Time 

Bank  Angle 

15.3 

2.2 

1.5 

Sideslip  Angle 

-0.006 

0.8 

•k** 

Roll  Rate 

15.0 

0.  5 

*  ** 

Yaw  Rate 

0.7 

1.1 

1.5 

TABLE  6-3 — Continued 


Flight  Condi. tj-jii;  1,2  Mach  at  15000  Feet 


Output 

Peak  Value** 

Peak  Time 

Settling  Time 

Bank  Angle 

15.9 

2.3 

3.7 

Sideslip  Angle 

-0.04 

1.4 

*** 

Roll  Rate 

15.1 

1.0 

Aikflr 

Yaw  Rate 

1.0 

1.1 

1.9 

Flight  Condition:  0.9 

Mach  at  50000 

Feet 

Output 

Peak  Value** 

Peak  Time 

Settling  Time 

Bank  Angle 

15.08 

2.2 

1.9 

Sideslip  Angle 

.02 

0.7 

*** 

Roll  Rate 

17.0 

0.5 

*  *  * 

Yaw  Rate 

1.9 

0.4 

1.8 

Notes 


1.  Figures  of  merit  for  the  roll  rate  and  the  yaw 
rata  are  estimated  for  the  response  plots. 

2.  Bank  angle  and  sideslip  angle  are  components  of 
the  output  vector 

*Commanded  Inputs  In  order  are: 

(1)  bank  angle  In  radians;  and 

(2)  aide  slip  angle  In  radians. 

’*'*Angles  are  converted  from  radians  to  degrees  prior 
to  calculating  the  figures  of  merit. 

***Final  value  of  the  output  is  zero.  Settling  time 
Is  undefined  using  the  definition  of  2  percent  of  final 
value . 


ROI  L  «HO  rw  RATES  ( DEG/SEC » 
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very  close  to  the  computed  value  of  0.668  degrees /second 
as  discussed  In  Chapter  6-6. 

Analysis  of  the  time  histories  for  the  other  flight 
conditions  also  show  good  aircraft  responses.  The  data 
is  presented  in  Appendices  G  and  H. 

6-9  Beta  Pointing  Maneuver  for 

6.7  Mach  at  Is', OOP  Feet 

The  aircraft  lateral  equations  of  motion  are  tabu¬ 
lated  in  Table, 6-1.  The  command  vector  v,  design  parameters, 
and  controller  matrices  for  this  maneuver  are  presented  in 
Table  6-4.  The  figures  of  merit  are  listed  in  Table  6-5. 

The  object  of  this  maneuver  is  to  sideslip  the  air¬ 
craft  while  maintaining  a  wings  level  flight  attitude.  The 
sideslip  angle  is  the  angle  formed  between  the  velocity 
vector  and  the  nose  of  the  aircraft.  Pilots  often  use  the 
maneuver  during  crosswind  landings.  The  aircraft  is  side¬ 
slipped  into  the  wind  to  maintain  the  desired  ground  track 
while  simultaneously  keeping  the  nose  of  the  aircraft 
pointed  toward  the  runway,  the  fuselage  of  the  aircraft 
parallel  to  the  centerline,  and  the  wings  in  level  flight. 
Beta  pointing  is  also  helpful  for  making  minor  corrections 
during  formation  and  air  refueling  and  is  used  in  certain 
combat  maneuvers  and  tracking  tasks. 

In  this  maneuver,  a  ramp- to- hold  input  signal 
commands  the  sideslip  angle  to  10  degrees  (see  Fig.  6-5) . 

In  addition,  the  bank  angle  is  commanded  to  zero. 
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TABLE  6 '•4 


DESIGN  PARAMETERS  AND  CONTROLLER  MATRICES 
FOR  LATERAL  CONTROLLERS 


MANEUVER;  Beta  Pointing 

COMMAND  VECTOR  v*i  0,  0,  0,  0 

1‘,  0.17453,  50,  50 

ALPHA;  1.0  EPSILON;  1.0 

MEASUREMENT  MATRIX:  0.3,  0.3 


PLANT  WITH  ACTUATOR  DYNAMICS  AND 
COMPUTATIONAL  TIME  DELAY 


Flight  Condition 

Sl^a** 

0.4 

Mach 

at 

Sea  Level 

0.1, 

1.5 

— g-i 

.3557E+00 

X 

-.1212E-I-02 

-.1034E+00 

.8116E+02 

0.9 

Mach 

at 

Sea  Level 

1.0, 

2.4 

.1134E+01 

-.7473E+01 

.2809E+00 

.3520E-I-02 

0.7 

Mach 

at 

15000  Feet 

0.4, 

1.5 

.8704E+00 

-.8187E+01 

■ 

.1106E+00 

.4682E+02 

1.2 

Mach 

at 

15000  Feet 

1.0, 

2.4 

.2023Ef01 

-.1213E+02 

.3942E+01 

.5983E+02 

o 

• 

so 

Mach 

at 

50000  Feet 

1.0, 

2.4 

. 8349E+01 

-.4954E+02 

2775E+01 

.1954E-H03 

Notes 


*Each  pulse  entry  in  v  has  four  parts; 

1.  The  time  in  seconds  the  input  reaches  steady  state; 

2.  Steady  state  value  (angle  in  radians); 

3.  Time  the  input  leaves  steady  state;  and 

4.  Time  input  reaches  zero. 

**Diagcnal  element  values  of  the  matrix. 
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TABLE  6-5 


LATERAL  FIGURES  OF  MERIT 


MANEUVER;  Beta  Pointing 

COMMAND  VECTOR  v^:  0,  0,  0,  0 

1,  0.17453,  50,  50 


PLANT  WITH  ACTUATOR  DYNAMICS  AND 
COMPUTATIONAL  TIME  DELAY 


Plight 

Condition:  0.4 

Mach  at  Sea 

Level 

Output 

Peak  Value** 

Peak  Time 

Settling  Time 

Bank  Angle 

-7.9 

1.7 

*** 

Sideslip  Angle 

10.0 

10.0 

2.5 

Roll  Rate 

-9.5 

1.2 

imu 

Yaw  Rate 

-10.5 

1.0 

2.0 

Flight  Condition:  0.9  Mach  at  Sea 

Level 

Output 

Peak  Value**  Peak  Time 

Settling  Time 

Bank  Angle 

-4.5  1.4 

*** 

Sideslip  Angle 

10.0  10.0 

3.3 

Roll  Rate 

-5.0  1.2 

*** 

Yaw  Rate 

-11.0  1.0 

2.0 

Flight  Condition:  0.7 

Mach  at  15000 

Feet 

Output 

Peak  Value** 

Peak  Time 

Settling  Time 

Bank  Angle 

-3.45 

1.3 

Sideslip  Angle 

10.0 

10.0 

2.8 

Roll  Rate 

-4.5 

0.5 

Yaw  Rate 

-10.0 

0.9 

2.1 

TABLE  6-5""Contlnuecl 


Flight 

Condition:  1.2 

Mach  at  15000 

Feet 

Output 

Peak  Value** 

Peak  Time 

Settling  Time 

Bank  Angle 

-6.6 

1.4 

*** 

Sideslip  Angle 

10.0 

10.0 

3.3 

Roll  Rate 

-7.5 

1.2 

*** 

Yaw  Rate 

•  -11.0 

1.1 

2.0 

Flight  Condition:  0.9 

Mach  at  50000 

Feet 

Output 

Peak  Value** 

Peak  Time 

Settling  Time 

Bank  Angle 

-0.7 

1.0 

*** 

Sideslip  Angle 

10.0 

10.0 

2.1 

Roll  Rate 

-4.0 

0.5 

*** 

Yaw  Rate 

-10.0 

0.9 

2.1 

Notes 


1.  Figures  of  merit  for  the  roll  rate  and  the  yaw 
rate  are  estimated  for  the  response  plots. 

2.  Bank  angle  and  sideslip  angle  are  components  of 
the  output  vector 

*Commanded  Inputs  In  order  are: 

(1)  bank  angle  In  radians;  and 

(2)  side  slip  angle  in  radians. 

^i^Angles  are  converted  from  radians  to  degrees  prior 
to  calculating  the  figures  of  merit. 

***Pinal  value  of  the  output  is  zero.  Settling  time 
Is  undefined  using  the  definition  of  2  percent  of  final 
value . 
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Figure  6-6  shows  the  aircraft  smoothly  follows  the  commanded 
beta  angle  without  an  overshoot.  As  expected,  the  rudder 
deflection  required  to  establish  the  sideslip  angle  causes 
a  negative  moment  about  the  x  eixls.  This  causes  the  bank 
angle  to  momentarily  change  during  the  transient  response 
until  corrective  aileron  Is  applied  by  the  control  system. 
For  this  flight  condition,  the  bank  angle  peaks  at  -3.45 
degrees  (see  Table  6-5)  before  returning  to  a  steady  state 
value  of  zero.  When  the  bank  angle  feedback  element  In  the 
sigma  matrix  is  Increased  to  reduce  the  maximum  bank  angle, 

I 

the  wings  begin  to  rock  In  an  attempt  to  maintain  a  zero 
bank  angle.  Pilots  often  perform  the  beta  pointing  maneuver 

with  manual  Inputs  to  the  aileron  and  rudder  In  conven- 
•  ► 
tlonal  aircraft  and  expect  the  bank  angle  to  change  momen¬ 
tarily  during  the  transient  portion  of  the  maneuver.  There¬ 
fore,  It  Is  preferred  to  accept  a  slightly  larger  error 
during  the  transient  as  long  as  It  Is  non-oscillatory. 

When  the  wings  rock  or  flutter  at  relatively  high  frequen¬ 
cies,  pilots  may  falsely  perceive  that  the  aircraft  Is 
unstable  and  uncontrollable  which  is  not  desirable. 

Figure  6-7  shows  the  time  history  of  the  rudder  and 
flaperon  deflection  angles.  Although  the  plots  are  not  as 
smooth  as  those  for  the  coordinated  turn,  the  rudder  Is 
still  operating  within  the  rate  limit  of  141  degrees  per 
second.  The  resulting  roll  and  yaw  rates  are  well  con¬ 
trolled  as  shown  in  Fig.  6-8. 
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Analyais  of  the  plots  for  this  maneuver  at  the 
other  four  flight  conditions  shows  good  aircraft  responses. 
The  data  Is  presented  In  Appendices  G  and  H. 

6-10  Summary 

This  chapter  begins  with  a  description  of  the' design 
procedure  used  to  obtain  the  lateral  control  laws.  Empha¬ 
sis  Is  placed  on  an  analysis  of  the  aircraft  open-loop 
transfer  functions  Including  the  location  of  the  zeros  and 
poles  and  the  controllability  and  observability  of  the  sys¬ 
tem.  Also  discussed  Is  the  effect  of  the  measurement 
matrix  on  the  transmission  zeros  and  the  transient  response. 
Than  the  data  for  the  coordinated  turn  and  beta  pointing 
maneuvers  are  presented  for  the  flight  condition  of  0.7 
Mach  at  15,000  feet.  The  chapter  concludes  with  a  brief 
description  of  the  maneuvers  and  an  analysis  of  the  time 
history  plots. 
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7.  ConcluBloriB  and  Recommendationa 


7"1  CoOTtianta  on  the  Raaulta 

The  major  objective  of  thla  thbaia  la  aocompllahed 
with  the  aucceaaful  deaign  of  flight  control  lawa  for  the 
Btatlcally  unatable  longitudinal  and  lateral  modea  on  the 
X-‘29  aircraft.  The  control  lawa  are  Implemented  with  con¬ 
troller  matricea  that  have  oonflguratlona  and  galna  deter¬ 
mined  by  applying  the  multivariable  control  theory  of 
Profeaaor  Brian  Porter. 

Very  good  reaulta  are  obtained  with  and  without 
actuator  dynamic a  and  computational  delay  added  to  the  aya- 
tem.  The  time  reaponaea  are  amooth  and  well  controlled 
and  meet  eatabllahed  requirements  auoh  aa  control  aurfaoe 
deflection  and  rate  llmlta.  The  daalgna  are  not  optimum 
in  every  aenae,  but  the  Porter  technique  provides  a  rela¬ 
tively  straightforward  and  quick  procedure  for  developing 
the  control  laws.  The  method  has  many  advantages  over  other 
multivariable  deaign  techniques  from  a  practical  applica¬ 
tion  point  of  view.  These  advantages  are  listed  as  followat 

1.  Analog  and  digital  designs  are  developed  using 
eaaentially  the  same  design  procedure. 

2.  The  method  uses  output  feedback  without  requir¬ 
ing  full  state  information. 
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3.  A  high  degrae  of  decoupling  between  the  outputs 
is  easily  obtainable. 

4.  Computer  assisted  software  provides  the  capa¬ 
bility  to  iteratively  fine  tune  the  design  to  obtain  the 
best  possible  responses. 

Although  there  are  many  advantages  to  the  Porter 
method,  this  research  effort  encounters  some  limitations 
which  must  be  overcome  to  make  this  technique  a  truly 
superior  way  to  solve  actual  aircraft  flight  control  prob¬ 
lems.  These  limitations  lare  discussed  in  the  next  section. 

7-2  Improving  the  Design  Process 

Refinement  and  expansion  of  Porter's  control  law 
theory  is  desirable  to  overcome  some  problems  that  are 
encountered  in  this  research  effort  when  applying  the  Porter 
method  in  a  practical  way  to  an  aircraft  with  serious  flight 
control  problems. 

First,  theory  should  be  developed  to  determine  a 
method  of  selecting  initial  sigma  weighting  matrices  for  a 
non-minimum  phase  plant  that  will  result  in  a  stable 
response.  Currently,  trial  and  error  is  the  only  method 
available,  which  proves  to  be  very  time  consuming  and  not 
always  successful. 

Another  limitation  occurs  when  actuator  dynaini'- 
are  implemented  into  the  state  equations  of  the  aircraft 
plant.  Then  a  measurement  matrix  becomes  necessary  and 
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derivatives  of  the  states  must  be  fed  back.  This  is  often 
undesirable.  For  example,  in  this  research  effort  the  mea¬ 
surement  matrix  generates  additional  transmission  zeros 
which  are  in  the  right  half  s-plane  and  the  result  is  an 
inability  to  achieve  a  stable  design.  See  Appendix  I  for 
additional  details.  Until  theory  can  be  expanded  to  permit 
a  feedforward  D  matrix,  acceleration  cannot  be  assigned  as 
an  output  variable  without  implementing  actuators  into  the 
state  matrix. 

Third,  the  Porter  method  can  be  more  useful  in 
practical  application  if  some  of  the  inputs  can  float  for 
certain  maneuvers.  Currently,  a  Porter  "configured” 
flight  control  system  requires  every  input  to  have  a  spe¬ 
cific  command  irrespective  of  the  maneuver.  Thus,  for  some 
maneuvers,  nonessential  inputs  require  precomputation  or 
the  output  C  matrix  must  be  changed. 

Finally,  the  theory  needs  expanding  so  weighting 
can  be  assigned  to  control  surfaces  prior  to  calculating 
the  control  matrices.  This  will  eliminate  overdriving 
leas  effective  control  surfaces  when  other  control  surfaces 
still  have  unused  control  power  remaining. 

7-3  MULTI  Improvements 

Recommendations  to  the  MULTI  program  are  based  on 
its  extensive  use  during  this  research  effort. 


Ultimately  the  MULTI  program  should  be  rewritten 
using  structvired  programming  techniques  along  with  adequate 
documentation. 

The  ZERO  program,  which  calculates  the  input  decou¬ 
pling,  output  decoupling,  and  transmission  zeros,  should  be 
added  to  the  MULTI  program.  This  will  save  time  and 
eliminate  duplicate  entries  of  system  matrices  into  each 
prograun. 

An  option  should  be  added  to  MULTI  to  provide  the 
capability  to  plot  aircraft  control  input  rates  on  the 
terminal  screen  and  Calcomp  plotter.  Also,  control  input 
rate  limiting  needs  to  be  added  to  the  simulation. 

Option  28  needs  estpansion  so  figures  of  merit  can 
be  calculated  not  only  for  the  outputs  but  also  for  the 
states,  any  combination  of  states,  and  the  control  inputs. 

7-4  Proposed  Future  Work 

In  general,  it  is  highly  desirable  to  control 
rates  and  accelerations  as  opposed  to  angles  which  are 
often  difficult  to  measure.  Angle  information  is  generally 
provided  by  inertial  navigation  systems  which  have  lower 
hardware  reliability  than  the  flight  control  system.  Also 
pilots  control  rates  in  conventional  aircraft  and  feel 
accelerations. 

Therefore,  future  work  on  the  X-29  should  concen¬ 
trate  on  controlling  rates  and  accelerations.  Rates  as 
output  variables  generally  cause  transmission  zeros  to 


139 


ooaur  at  the  origin  which  uaually  makes  stabilizalton  o£  the 
system  more  tedious  but  not  Impossible. 

An  unsuccessful  attempt  to  control  acceleration  in 
this  research  effort  results  from  the  inability  to  find  a 
stable  design  when  actuators  are  implemented  in  the  state 
equations.  Future  research  efforts  should  attempt  to  over¬ 
come  this  problem  by  finding  an  easy  way  to  implement  a 
D  matrix  into  the  design  procedure. 

Limiting  the  bandwidth  of  the  control  system  reduces 
suacept;ibility  to  noise  and  other  objectionable  high  fre¬ 
quency  disturbances.  Therefore^  it  is  desirable  to  investi¬ 
gate  the  loop  transmission  frequency  of  the  final  control 
system  design. 


Appendix  A:  Multivariable  Control  Theory 
(Reproduced  from  Reference  4) 

This  thesis  utilizes  the  multivariable  design  method  of 
Professor  Brian  Porter  of  the  University  of  Salford;  England 
(10) .  The  design  method  uses  output  feedback  with  high-gain 
error-actuated  controllers.  Output  feedback  is  advantagous 
since  state  variables  may  be  difficult  to  measure  while 
system  response  data  are  more  readily  available. 

System  State  Equations 

Porter's  method  works  equally  well  for  either  contin¬ 
uous  or  discrete  systems;  but  it  is  often  easier  to  first 
examine  a  system  in  the  continous  time  domain.  This  is 
because  designs  can  be  more  easily  visualized  in  the  s-plane 
than  in  the  z-plane.  A  continuous  time  system  is  represent¬ 
ed  by  the  state  space  model t 

e 

X  ■  Ax  +  BU 

y  ■  Cx  (A-1) 

where 

A  ■  continuous  plant  matrix  (n  x  n) 

B  ■  continuous  input  control  matrix  (n  x  m) 

C  -  continuous  output  matrix  (  i  x  n) 

X  ■  state  variable  vector  with  n  states 
u  ■  input  vector  with  m  inputs 
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^  ■  output  vector  with  l  outputs 


The  system  Inputs  for  an  aircraft  are  the  control  surface 

deflections  or  actuator  input  commands,  and  the  system 

outputs  are  aircraft  responses  affected  by  the  inputs. 

» 

The  method  does  not  allow  for  a  feedforward,  D,  matrix. 
If  such  a  matrix  is  present  in  the  original  state  space 
model,  the  control  inputs  must  be  redefined  as  states  so 
that  the  D  matrix  is  absorbed  into  the  C  matrix.  This  can  be 
accomplished  by  incorporating  the  actuator  dynamics  into  the 
plant  model.  Actuator  inputs  then  become  control  inputs. 

To  employ  Porter's  method,  it  is  desirable  (but  not 
necessary)  to  partition  the  system  state  equations  as 
follows: 


(A-2) 


The  equations  are  partitioned  so  that  B2  and  C2  are  square 
(m  X  m)  and  (  i  x  i )  matrices,  respectively.  The  method 
requires  that  the  number  of  inputs  to  the  system  equals  the 
number  of  outputs  which  means  m  »  i  ,  and  therefore  the 
dimension  of  B2  equals  the  dimension  of  C2.  It  is  often 
possible  to  form  the  state  equations  so  that 
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For  the  discrete  case  the  system  equations  are  written 
as  follows: 

x[(k+l)T]  =  ^(kT)  +  _!j^{kT) 

y(kT)  -  rx(kT)  (A-3 

where 

^  ■  exp (AT)  ■  discrete  plant  matrix 

jfe,  ■  I  exp(AT)Bdt  ■  discrete  input  control  matrix 
0  -  - 

r  ■  c  »  discrete  output  matrix 

In  the  above  equations  T  is  the  sampling  period^  and  k 
takes  on  integer  values  from  zero  to  plus  infinity. 

System  With  Output  Feedback 

Figure  A-l  shows  the  block  diagram  for  a  continuous 
output  feedback  system,  where  v  is  the  command  input  vector 
and  y  is  the  desired  output  vector.  The  blocks  for  the 
plant  are  derived  directly  from  the  system  state  equations, 
Equation  (A-l) .  The  proportional  plus  integral  controller 
has  three  parameters.  Kg,  and  g,  which  must  be  deter¬ 
mined  by  the  designer.  The  output  signal  of  the  controller 
u,  is  given  in  the  following  control  law  equation: 

ii  *  g(Kge  +  Kj^/edt)  (A-4 

where 

u  is  the  output  signal  of  the  controller 


Figure  A-1.  System  block  Diagram  —  Continuous  Case 


e  is.  the  error  signal  at  the  input  of  the  controller 
Kq  is  the  proportional  gain  matrix' 

is  the  gain  matrix  for  the  integral  term 
g  is  the  scalar  forward  path  gain 

A  measurement  matrix  M  is  included  in  the  system  if  the 
plant  is  irregular.  Regular  and  irregular  plants  are  dis¬ 
cussed  later. 

The  discrete  system  block  diagram,  shown  in  Figure  A-2, 
is  similar  to  the  continuous  system,  but  Equation  (A-4) 
becomes 

u(kT)  -  (1/T) [KQe(kT)  +  K^z{kT)]  {A-5) 

where  the  forward  path  gain  g  equals  the  sampling  frequency, 
(1/T).  The  z,(kT)  matrix  is  derived  from  the  backward 
difference  equation, 

z[(k+l)T]  -  z(kT)  +  Te(kT)  (A-6) 

The  steps  to  be  taken  next  in  the  design  method  depend 
on  whether  or  not  [^]  has  full  rank,  i.e.,  does  it  have  an 
inverse.  If  the  matrix  [CB]  has  full  rank,  the  plant  is 
called  "regular"  and  no  measurement  matrix  M  is  needed. 
However,  if  [CB]  does  not  have  full  rank,  the  plant  is 
called  "irregular"  and  M  is  needed  to  form  a  new  matrix  [FB] 
(see  Equations  (A-12)  through  (A-14) )  which  does  have  an 
inverse.  This  is  explained  in  more  detail  in  the  next 
sections.  When  the  partitioned  B  matrix  in  Equation  (A-2) 


has  the  form 


then 


(A-7) 


[CB] 

(A-8) 

[PB] 

(A-9) 

Regular  Plant 

For  the  system  to  be  classified  as  ''regular''  the  first 
Markov  parameter  [^1  must  have  full  rank.  If  this  is  true 
the  gain  matricies  can  be  found  from 

Kq  "  (A-10) 

and 


(A-11) 


where 


a  is  a  constant  which  assigns  the  ratio  of  proportional 
to  inte^iral  control 

^  is  the  diagonal  weighting  matrix 

The  diagonal  weighting  matrix,  £_  »  diag  02'  •••  >  Of)  » 

is  specified  by  the  designer.  Each  (i  =  deter¬ 

mines  the  weighting  of  the  effect  of  a  particular  error 
signal  on  each  control  input. 
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Irregular  Plant 

If  the  first  Markov  parameter  [^]  is  rank  deficient/ 
then  the  plant  is  called  "irregular".  In  this  case,  the  C 
matrix  roust  be  replaced  by 

P  -  [Pj_  I  Pj]  (A-12) 

where 

£l  ■  +  Mill  <A-13) 

Pj  ■  [Cj  +  MA12I  (A-14) 

The  matrix  M  in  the  above  equations  is  a  measurement  matrix 
which  is  chosen  such  that  the  matrix  [PB]  has  full  rank. 

The  designer  chooses  the  measurement  matrix  so  that  it  is  as 
sparse  as  possible/  thus  the  smallest  number  of  additional 
measurements  are  required.  Reference  13  gives  an  approach 
for  selecting  the  measurement  matrix  to  achieve  optimal 
decoupling.  Once  M  is  formed/  Kq  and  are  computed  by 


^0 

•  [PB]“^£. 

(A-15) 

h 

-  otPBl'^i 

(A-16) 

which  are  similar  to  Equations  (A-10)  and  (A-11) . 

Por  irregular  plants  the' error  vector  e  is  defined  as 

eav-w  (A- 17) 
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where 


(A-18) 


For  step  Inputs  the  values  of  the  rates r  become  zero  In 
the  steady  state  because  they  represent  kinematic  variables. 

AsvnMPtotio  Characteristics 

As  the  gain  factor  of  the  system,  9  (or  1/T  for  the 
discrete  case) ,  approaches  Infinity,  the  system  transfer 
function  matrix  G(s)  assumes  the  asymptotic  form 

r(X)  -  Ed)  +  (A-19) 


where 

r(X)  Is  the  slow  transfer  function  matrix 

r(X)  Is  the  fast  transfer  function  matrix 

The  roots  of  the  asymptotic  closed-loop  transfer  func¬ 
tion  may  be  grouped  into  three  setst  ^2'  I3'  Table 

A-1  gives  the  equations  for  finding  these  asymptotic  roots. 
Sets  and  Z2  correspond  to  the  slow  modes  of  the  system, 
where  the  inodes  associated  with  the  roots  in  become 
uncontrollable,  and,  for  regular  plants,  the  modes  associat¬ 
ed  with  the  roots  in  ^2  become  unobservable  as  the  gain 
increases.  Set  Z^,  the  infinite  roots,  are  associated  with 
the  fast  modes  of  the  system  which  become  dominant  as  the 
gain  increases. 

The  roots  in  set  correspond  to  the  transmission 
zeros  of  the  system  which  are  not  altered  by  output  feed¬ 
back.  Since  infinite  gain  cannot  be  implemented,  and  is  not 
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desirable,  the  closed-loop  roots  of  the  system  tend  to 

migrate  toward  the  transmission  zeros.  This  may  adversely 

affect  the  system  stability  if  the  location  of  these  zeros 

Is  In  the  unstable  region.  Reference  12  gives  a  procedure 

for  locating  the  transmission  zeros  of  a  system. 

As  the  gain  Increases  the  system  output  responses 

become  Increasingly  decoupled.  The  asymptotic  closed-loop 

transfer  function  for  the  continuous  case  has  the  form 

,  .  90,  go, 

r(X)  -  diag  { - (A-20) 

For  the  discrete  case  the  form  Is 

r(X)'-diag{ - - 5^}  (A-21) 

where  the  0^  (1  ■  1 are  determined  by  the  weighting 
matrix,  jC,. 

The  computer  program  MULTI  Is  of  great  value  In 
reducing  the  time  required  to  achieve  a' satisfactory  design. 
The  MULTI  User's  Manual  (8)  describes  the  program  and  Its 
operation. 
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Appendix  Bi  Forming  G(0)  and  G (0) Matrices  for 
Calculating  Steady  State  Vectors 


Introduction 

This  additional  Improvement  to  MULTI  provides  the 
capability  to  calculate  the  G(0)  and  2(0)**^  matrices 
directly  from  A,  B,  and  C  matrices  of  the  plant.  Also, 
the  user  is  given  the  option: 

1.  to  enter  the  plant  steady  state  output  vector 
^  and  have  MULTI  automatically  calculate  the  steady  state 
Input  vector  u  using  the  G(0)  Inverse  matrix; 

2.  to  enter  the  plant  steady  state  Input  vector 
u  and  have  MULTI  automatically  calculate  the  steady  state 
output  vector  ^  using  the  G(0)  matrix. 

User' s  Guide 

First  the  A,  B,  and  C  matrices  of  the  plant  are 
entered  into  MULTI  using  Options  2  and  3  or  9.  Then  the 
user  selects  Option  6  and  the  following  prompt  appears: 

FOR  OPEN-LOOP  TRANSFER  FUNCTION  ENTER  1  > 

FOR  CLOSED-LOOP  TRANSFER  FUNCTION  ENTER  2  > 

FOR  G(0)  AND  G(0)  INVERSE  MATRICES  ENTER  3  > 


If  3  is  entered,  the  matrices  are  computed  and  printed 
like  the  example  which  follows: 


THE  G(0)  MATRIX  IS: 


-.3192E-01 

.8212E+02 

-.1039E-01 


.1453E-01 

-.3726E+02 

.3564E-02 


.8929E+00 

-.1225E+03 

.1213E-01 


THE  G(0)  INVERSE  MATRIX  IS: 


.1922E+01 

-.3467E+01 

.1183E+01 


-.3769E-01 

-.1114E+00 

.4658E-03 


-.3949B-I-03 

-.8705E+03 

.4607E-01 


Next  the  following  prompt  automatically  appears: 


TO  USE  THE  G(0)  INVERSE  MATRIX  OR  G(0)  MATRIX 
TO  CALCULATE  THE  STEADY  STATE  VALUE  OF  THE: 

PLANT  INPUT  VECTOR  U  ENTER  >1 
PLANT  OUTPUT  VECTOR  Y  ENTER  >2 
TERMINATE  CALCULATIONS  ENTER  >  3 


If  1  is  entered  the  following  prompt  appears: 


ENTER  THE  ”N''  ELEMENTS  OP  THE 
PLANT  OUTPUT  VECTOR  Y: 

For  example/  enter  1/1/1  and  the  following  appears; 


THE  PLANT  OUTPUT  VECTOR  Y 
JUST  ENTERED  IS; 


1. 

1. 

1. 


THE  PLANT  INPUT  VECTOR 
U  »  G(0)  INVERSE  *  Y  IS 


-394.71625 

-874.03421 

1.22980556 


th«n  the  following  prompt  automatically  roappeare  and  the 
prooeie  le  repeated  until  the  oaloulatlone  are  terminated 
by  entering  3: 


TO  USE  THE  G(0)  INVERSE  MATRIX  OR  G(0)  MATRIX 
TO  CALCULATE  THE  STEADY  STATE  VALUE  OF  THEi 

PLANT  INPUT  VECTOR  U  ENTER  1 

PLANT  OUTPUT  VECTOR  Y  ENTER  2 

TERMINATE  CALCULATIONS  ENTER  3 

Programmer* ■  Guide 

The  following  FORTRAN  77  code  Is  added  to  the 
MULTI  Option  6  Overlay.  Comments  explain  each  portion  of 
the  program. 


MINTt,  *m  6(0)  AN0  6(0)  INViRSE  HATRICEB  ENTER  3  >  ' 

READ!,  TFTYRE 

ctiMimtttMontMMHtttmtiitttttmtitMuttttitiMiututttittnit 

etIltIMtITHlS  BECTION  CALCUUTfil  m  6(0)  AND  6(0)  INVERSE  NATRtCESUtttt 
IF(TmRE.CQ.3)  THEN 
60  la  I>l,20 
IT(I)>0. 

CT(I)*0. 

142  CONTINUE 

CMIlirOLUNIINE  SnS  DUHHY  VARIA61I  L«N  (01HEN8IDN  OF  NATRII  A)tlttMttll 
(.•N 

CtMMFOUOSINE  COPIES  HATRIX  A  TO  OUNRY  RATRIX  VARIAME  ACLttttttttttttt 
DO  144  (■I.N 
N  143  0-l,N 

ACl(l,J)«A(IW) 

143  CONTINUE 

144  CONTINUE 

CltlllFOUONINS  CALCULATES  DENOMINATOR  POLYNOHIAL  OF  ClIBl-A]  INVERSE  I  S 
CAU  FNOFS  (ACL.L.L.DET) 

CtltltF0UINlIN6  CALCULATES  NUHERATOR  POLYNOHIAL  OF  CtCBHI  INVERSE  t  I 
BO  170  IPM.P 
SO  14T  Iiei.N 

so  146  1-1,N 

ST(1)*I(MN) 

CTdKdM) 


CM.L  CAOII  (CT,IT,6F,L» 

CttUtttFOLLOilINfi  CAICULATEB  ELERENTS  OP  MTRII  KOIttttttttMttttttttttl 
62IIP,lHI-6P(l)yOETU) 

169  CONTINUE 

170  CONTINUE 
PN1NTI,»  * 

PRINTI  '  ' 

PRINT VTNE  1(0)  IWTRIX  III* 

PRINT!,'  » 

PRINT!,'  ’ 

CALL  IMTPR(iZ,P,H) 

PRINT!,'  ' 

PRINT!,  'THE  1(0)  INVERIE  NRTRII  111' 

PRINT!,'  ' 

PRINT!,'  ' 

CIII!!t!!!POLLOHlNG  CALCUIATEI  NRTRII  6(0)  lNVER8St!tt!l!l!!l!l!!!!!t!l!! 
lOBTilO 
tR«10 

CALL  LIRV2P  (62,P,lA,i21HV,lB6T,MI(AREA,IIR) 

IP  (IIR.n.S!)  THEN 

PRINT!, 'ACCURACT  TEST  PAtLCO....fiATA  HAY  IE  SUSPECT' 

ElK  IP  (IIR.E0.t31)  THEN 

PRINT!, 'HATRII  TOO  ILL-CONDITIONED  FOR  ITERATIVE* 

PRINT!, 'IHPROVEHEKT  TO  IE  EFFECTIVE' 

SLII  IP  (IIR.ED.12T)  THEN 

PRINT!, 'tlNGULAR  AND  CANNOT  SC  INVERTED' 

ELK 

CAU  NATPR  (I2INV,P,H) 

PRINT!,'  • 

PRINT!,'  • 

171  PRINTI, 'TO  USE  THE  6(0)  INVERSE  HATRII  OR  6(0)  HATRII' 
PRIHTI,'TO  CALCULATE  THE  STEADY-STATE  VALUE  OF  THEl’ 

PRINT!,'  ' 

PRINT!, 'PLANT  INPUT  VECTOR  U  ENTER  >1' 

PRINT!, 'PLANT  OUTPUT  VECTOR  Y  ENTER  >2' 

PRINT!, 'TERHINATE  CALCULATIONS  ENTER  >3' 

REA0!,I1YE8 
IP  (IIYCS.E0.1)  THEN 
PRINT!,'  ' 

PRINTI, 'ENTER  THE  ',P,'  CDLUHN  ELCHENT8  OF  THE' 

PRINTI, 'PLANT  OUTPUT  VECTOR  Y' 
READI,(IYYY(I,J),J«t,t),I>l,P) 

PRINT!  '  ' 

PRINTI/THE  PLANT  OUTPUT  VECTOR  V 
PRINTI, 'JUST  ENTERED  ISi' 

PRINT!,'  ' 

PRINT!,'  ’ 
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DO  172  I-I.P 

PRINTt/ 

CORTINUE 

1A«10 

CALL  VHUlFF(6ZINV,m,P,P,J,IA,IA,UUU,IA,IERl 
PRIRTI,'  ’ 

PRlNTt  ’  * 

PRINTVTHE  PLANT  INPUT  VECTOR* 

PRINTt/U  •  6(0)  INVERSE  t  y  IS* 

PRINT!,’  ’ 

PRINT!,'  ’ 

DO  173  I*1,P 

PRINT!,*  ’,UIJU(1,1) 

CONTINUE 
PRINT!,’  ’ 

PRINT!,’  ’ 

60  TO  171 

ELSE  IP  (IIYES.Ee.Z)  THEN 
PRINT!  '  ' 

PRINT!! ’ENTER  THE  *,P.’  COLUHN  ELEHENT8  OF  THE’ 
.  PRINT!, ’PLANT  INPUT  VECTOR  U’ 
READ!,((UUU(1,J),>1,1),1«1,P) 

PRINT!  ’  ’ 

PRINT!*, ’THE  PLANT  INPUT  VECTOR  U*' 

PRINT!, ’JUST  ENTERED  IS* 

PRINT!,’  * 

PRINT!,’  * 

00  174  1*1,P 

PRINT!,’  *,UUU(I,t) 

CONTINUE 

lA-10 

CAU  VHULFF(&Z,UUU,P,P,J,IA,1A,YYY,IA,IER) 
PRINT!,’  ’ 

PRINT!  ’  ’ 

PRINT!! ’THE  PLANT  OUTPUT  VECTOR* 

PRINT!, ’Y  «  6(01  !  0  IS’ 

PRINT!,’  ’ 

PRINT!,’  ' 

DO  179  1M,P 

PRINT!,’  ’,YYY(I,1I 

CONTINUE 
PRINT!,’  ’ 

PRINT!,’  ' 

60  TO  171 

END  IF 
IF 

60  TO  S017 


Appendix  C:  Converting  Radiana  to  Degrees 

for  Plotting 

f 

Introduction 

This  addition  to  MULTI  provides  the  capability  to 
selectively  convert  the  vectors  and  states  of  the  simula¬ 
tion  output  from  radians  to  degrees  prior  to  plotting  the 
data  or  calculating  the  figures  of  merit.  This  is  useful 
because  the  equations  of  motion  for  aircraft  often  have 
the  variables  describing  angles  with  dimensional  units  in 
radians  and  control  surface  deflections  in  dimensional 
units  of  degrees. 

User' s  Guide 

Option  26  is  executed  to  generate  the  simulation 
data.  After  the  calculations  are  completed  the  following 
prompt  appears: 

CONVERT  RADIANS  TO  DEGREES?  1-YES,  O-NO 

If  0  is  entered  the  option  is  skipped.  If  1  is  entered, 
then  prompts  appear  as  shown  in  the  following  example: 

ENTER  #  OF  OUTPUTS  (Y  VECTOR)  TO  CONVERT 
?  2 

ENTER  INTEGERS  FOR  Y  VECTOR  CHOICES 
?  1,3 


ENTER  #  OF  INPUTS  (V  VECTOR)  TO  CONVERT 
?  2 

ENTER  INTEGERS  FOR  V  VECTOR  CHOICES 
?  1,3 

ENTER  #  OF  .STATES  TO  CONVERT 
?  3 

ENTER  INTEGERS  FOR  STATE  VECTOR  CHOICES 
?  1,3,4 


Programmer ' a  Guide 

The  following  FORTRAN  77  code  is  added  to  the 
Option  26  overlay.  Comments  explain  each  section  of  the 
code. 


PMKTt/CONVEftT  RAfilAMS  TO  DEBREES?  1-YES,  O-NO* 
REMt.IYES 
IF(1YES.IIE.1)  THEN 
BO  TO  0990 
HBE 

PR1NTI«'CNTER  I  OF  OUTFITS  (Y  VECTOR)  TO  CONVERT* 
REMI.NUNBER 
END  IF 

lF(liURBER.CB.O)  THEN 
BO  TO  0830 
asE 

PRINTS, 'ENTER  INTESERS  FOR  Y  VECTOR  CHOICES*. 

READS, (ICLHS(BN),BN>I.NUH6ER) 
n  8800  QN>1, NUMBER 
DO  8700  M.NT 

YP(I,IICLHS(gN)41>)«YPII,(ICLHS(QN)M)  >137.2930 
870*^  CONTINUE 

8800  CONTINUE 
END  IP 

8130  PRINTS, 'ENTER  t  OF  INPUTS  (V  VECTOR)  TO  CONVERT* 

READS, NUH8ER 
IF(NUnBER.Efl.0>  THEN 
60  TO  1900 

asE 

PRINTS, 'ENTER  INTEGERS  FOR  V  VECTOR  CHOICES' 

READS, IIClN3t0N),0N«l,NUN8ER) 

DO  8870  ON* I, NUMBER 
00  8860  1«1.NT 

VP(I,(]CLHS(8N)«l»*VP(I,liaM8l0N>M))S37.2938 
8860  CONTINUE 

8870  CONTINUE 
END  IP 


C-2 


8900  PRlNTt, ’ENTER  t  OF  STATES  TO  CONVERT’ 

READt.NUNBER 
1F(NUHBER,EB.0)  THEN 
60  TO  8990 
ELSE 

PRINTt, 'ENTER  1NTE6ERS  FDR  STATE  VECTOR  CHOICES’ 
READt,(tCLH8(0N),aN>t.NUNBER) 

00  6970  QNM.NUHBER 
DO  8760  IM.NT 

IP(I.(lCLNS(0N)M))>XP(I,(IClJIS(0N)il))}57.295B 

8960  CONTINUE 
8970  CONTINUE 
END  IF 

8990  CONTINUE 


Appendix  D:  Longitudinal  State  Space  Matrices 

Tables  in  this  appendix  list  the  linearized  longi¬ 
tudinal  state  space  matrices  of  the  X-29  aircraft  for  the 
following  flight  conditions: 

(1)  0.4  Mach  at  Sea  Level 

(2)  0.7  Mach  at  15000  Feet 

(3)  1.2  Mach  at  15000  Feet 


D~1 


TABLE  D-1 


LONGITUDINAL  STATE  SPACE  MATRICES 


FLIGHT  CONDITION:  0.4  Mach  at  Sea  Level 


.OOOOE+OO 

.3217E+02 

.OOOOE+00 

.OOOOE+00 


(Plant  Matrix) 


.OOOOE+00 

.1660Er01 

.3172E-03 

.2823E-03 


.OOOOE+00 

.4880E-01 

-.lOOlE+01 

.1012E+02 


(Control  Input  Matrix) 


.lOOOE+01 

.4273E+00 

.9880E-I-00 

'.4340E+00 


$  43 


.OOOOE+00 

'.4891E>01 

‘.9426E-03 

.9328E-01 


OOOOE+01 

OOOOE+00 

OOOOE+00 


•OOOOE+OO 

.2592E-01 

-.8021E-03 

-.3603E-01 


.OOOOE+OO 

.lOOOE+01 

.OOOOE+OO 


.OOOOE+OO 

.OOOOE+OO 

.lOOOE+01 


.OOOOE+OO 

.3779E+02 

•.6761E-02 

•.2679E+00 


.lOOOE+01 

.OOOOE+OO 

.OOOOE+OO 


Notes 

1.  States  (listed  in  order)  are  pitch  angle, 
change  in  forward  velocity,  angle-of-attack,  and  pitch 
rate . 

2.  Control  inputs  (listed  in  order)  are  can'ards 
strake  flaps,  and  thrust. 

3.  Outputs  (listed  in  order)  are  pitch  rate, 
change  in  forward  velocity,  and  angle  of  attack. 


TABLE  D-2 


LONGITUDINAL  STATE  SPACE  MATRICES 


FLIGHT 

CONDITION:  0.7  Mach  at  15000 

Feet 

A  (Plant  Matrix) 

.OOOOE+00 

-.3213E+02 

-.3243E-06 

.OOOOE+00 

.OOOOEtOO  .OOOOEtOO 

-.1199E-01  -.1718E+01 

-.1298E-03  -.1138E+01 

.29d6E-03  .1892E-f02 

.lOOOE+01 

-.3504E+00 

.9917E+00 

-.4928E+00 

B  (Control  Input  Matrix) 

.OOOOE-t-OO 

-.5895E-01 

-.1161E-02 

.1720E+00 

•OOOOE+OO 

.2621E-01 

-.7808E-03 

-.5630E-01 

. OOOOE+OO 
.2724E+02 
-.2099E-02 
-.1929E+00 

C  (Output  Matrix) 

.OOOOE+00 

.OOOOE+00 

.OOOOE+00 

.OOOOE+OO  .OOOOE+OO 
.lOOOE+01  .OOOOE+OO 
.OOOOE+OO  .lOOOB+01 

.lOOOE+01 

.OOOOE+OO 

.OOOOE+OO 

Notes 

1.  States  (listed  in  order)  are  pitch  angle, 
change  in  forward  velocity,  angle-of-attack,  and  pitch 
rate. 


2.  Control  inputs  (listed  in  order)  are  canards, 
strake  flaps,  and  thrust. 

3.  Outputs  (listed  in  order)  are  pitch  rate, 
change  in  forward  velocity,  and  angle  of  attack. 


TABLE  D-3 


LONGITUDINAL  STATE  SPACE  MATRICES 


PLIGHT  CONDITION: 

1.2  Mach  at  15000 

Feet 

A  (Plant  Matrix) 

.OOOOE+00 

.OOOOEtOO 

.OOOOE+OO 

.lOOOB+Ol 

3213Et02 

-.1083Et00 

-.2668E+02 

-.3814E+00 

-.1892E-06 

-.6824E-04 

-.2461E+01 

.9930E+00 

.OOOOEtOO 

-.7393E-02 

.1920E+02 

-.9340E+00 

B  (Control  Input  Matrix) 

.OOOOEtOO 

.OOOOE+00 

. OOOOE+OO 

-.2959E+00 

. 3583E-01 

.1079E+03 

3266E-02 

-.7448E-03 

3645E-02 

.  3880E-t-00 

-.1195E+00 

-.7638E+00 

C  (Output  Matrix) 

.OOOOE+00 

•OOOOE+OO 

.OOOOE+OO 

. lOOOE+01 

.OOOOE+00 

.lOOOE+01 

.OOOOE+OO 

. OOOOE+OO 

.OOOOE+00 

.OOOOE+OO 

.lOOOE+01 

.OOOOE+OO 

Notes 

1. 

States  (listed 

in  order)  are  pitch 

angle , 

change  in  forward  velocity. 

angle-of-attack,  and  pitch 

rate . 

2. 

Control  inputs 

(listed  in  order)  are  canards. 

strake  flaps 

,  and  thrust. 

3.  Outputs  (listed  In  order)  are  pitch  rate, 
change  in  forward  velocity,  and  angle  of  attack. 


Appendix  E:  Longitudinal  Design  Parameters i 

Controller  Matrices ,  and 
Figures  of  Merit 

Two  seta  of  tables  in  this  appendix  tabulate  the 
design  parameters/  controller  matrices/  and  figures  of 
merit  for  the  following  longitudinal  maneuvers  and  flight 
conditions: 

Direct  Climb 

1 .  0.4  Mach  at  Sea  Level 

2.  1.2  Mach  at  15000  Feet 

Vertical  Translation 

1.  0.7  Mach  at  15000  Feet 

2.  1.2  Mach  at  15000  Feet 

Pitch  Pointing 

1.  0.7  Mach  at  15000  Feet 

2,  1.2  Mach  at  15000  Feet 


TABLE  E-1 


DESIGN  PARAMETERS  AND  CONTROLLER  MATRICES 
FOR  LONGITUDINAL  CONTROLLERS 


FLIGHT  CONDITION:  0.4  Mach  at  Sea  Level 

MANEUVER:  Direct  Climb 

COMMAND  VECTOR  v*:  1.1,  0.0436,  3,  4.1 

0 ,  0  ^  0 ,  0 

0,  0,  0,  0 


Plant  Only 


Alpha 

Epsilon 

Sigma** 

K  *** 

-0 

1.  7 

.1253E+02 

-.8427E-02 

-.3313E+01 

0.714 

1.4 

1.2 

.1495E+02 

-.2592E+00 

-.8588E+01 

0.01 

.2647E-01 

. 3192E-01 

.1603E-02 

Plant  and  Computational 

Time  Delay 

Alpha 

Epsilon 

Sigma** 

Kq*** 

1.9 

.1400E-t-02 

-.8427E-02 

-  .6625E+01 

0.714 

1.4 

1.2 

.1670E+02 

-.2592E+00 

-.1718E+02 

0.02 

.2958E-01 

.3192E-01 

.3206E-02 

Plant  with 

Actuator  : 

Dynamic  s 

Alpha 

Epsilon 

S iqma  *  * 

V  ★** 

-0 

2.0 

.1474E+02 

-. 8427E-02 

-.3313E+01 

0.714 

1.4 

1.2 

.1758E+02 

-. 2592E+00 

-.8588E+01 

0.01 

. 3114E-01 

. 3192E-01 

.1603E-02 

Plant 

with  Actuator  Dynamics  and  Computational 

Time  Delay 

Alpha 

Epsilon 

Sigma** 

K  *** 

-0 

2.1 

.1548E+02 

-.7022E-02 

-.4969E+01 

.555 

1.8 

1.0 

.1846E+02 

-.2160E+00 

-.1288E+02 

0.015 

.3269E-01 

.2660E-01 

. 2405E-02 

Notes 

*^Each  pulse  entry  in  v  has  four  parts: 

1.  The  time  in  seconds  the  input  reaches  steady- state; 

2.  Steady-state  value  (angles  in  radians); 

3.  Time  the  input  leaves  steady- state;  and 

4.  Time  the  input  reaches  zero. 

**Diagonal  element  values  of  the  matrix. 

***Ko  *■  alpha  * 

E—  2 


TABLE  E-2 


DESIGN  PARAMETERS  AND  CONTROLLER  MATRICES 
FOR  LONGITUDINAL  CONTROLLERS 


FLIGHT  CONDITION:  1.2  Mach  at  15000  Feet 

MANEUVER:  Direct  Climb 

COMMAND  VECTOR  v*:  1.1,  0.0436,  3.0,  4.1 

0,  0/  0,  0 

0,  0,  0,  0 


Plant  Only 


Alpha 

Epsilon 

Sigma** 

V  *** 

-0 

1.7 

.1866E-«-01 

.1094E-02 

-.1760E+01 

0.7143 

1.4 

0.6 

.8219E+01 

-. 3208E-01 

-.5695E+01 

0.01 

.7845E-02 

.5574E-02 

-.2935E-02 

Plant  and  Computational 

Time  Delay 

Alpha 

Epsilon 

Sigma** 

Kq*** 

1.9 

.20a5E+01 

.1094E-02 

-.3520E+01 

0.7143 

1.4 

0.6  - 

■.9186E+01 

-. 3208E-01 

-.1139E+02 

0.02 

.8768E-02 

.5574E-02 

-.5870E-02 

Plant  with 

Actuator  Dynamics 

Alpha 

Epsilon 

Sigma** 

Kq*** 

• 

2.0 

.2195E+01 

.1094E-02 

-.1760E+01 

0.7143 

1.4 

0.6 

. 9669E+01 

-.3208E-01 

-.5695E+01 

0.01 

.9230E-0Z 

.5574E-02 

-.2935E-02 

Plant 

with  Actuator  Dynamics  and  Computational 

Time  Delay 

Alpha 

Epsilon 

Sigma** 

K  *** 

-0 

2.1 

.2304E+01 

.1094E-02 

-. 2639E+01 

.  555 

1.8 

0.6 

.1015E+02 

-.3207E-01 

-.8542E+01 

0.015 

.9690E-02 

.5574E-02 

-.4402E-02 

Notes 

^Each  pulse  entry  in  v  has  four  parts: 

1.  The  time  in  seconds  the  input  reaches  steady- state; 

2.  Steady-state  value  (angles  in  radians); 

3.  Time  the  input  leaves  steady- state ;  and 

4.  Time  the  input  reaches  zero. 

**Diagonal  element  values  of  the  matrix. 

***Kq  -  alpha  *  K^. 


TABLE  E-3 


DESIGN  PARAMETERS  AND  CONTROLLER  MATRICES 
FOR  LONGITUDINAL  CONTROLLERS 


FLIGHT  CONDITION:  0.7  Mach  at  15000  Feet 
MANEUVER:  Vertical  Translation 
COMMAND  VECTOR  v*:  0,  0,  0,  0 

0,  0,  0,  0 

1.1,  -0.0349,  50,  50 


Plant  Only 


Alpha 

Epsilon 

Sigma** 

-0 

.4545 

2.2 

1.9 

1.0 

0.015 

.7435E+01 

.1113E+02 

.2679E-01 

.5984E-02 
-.1079E+00 
. 3683E-01 

-.4229E+01 

-.1293E+02 

.3290E-02 

Plant  and  Computational  Time  Delay 

Alpha 

Epsilon 

Sigma** 

V  *** 

-0 

• 

.  4545 

2.2 

1.9 

1.0 

0.03 

.7435E+01 

.1113E+02 

.2679E-01 

.5984E-02 
-.1079E+00 
. 3683E-01 

-.84588+01 

-.2586E+02 

.6581E-02 

Plant  with 

Actuator  Dynamics 

Alpha 

Epsilon 

Sigma** 

V  a.** 

-0 

.4545 

2.2 

2.0 

1.0 

0.015 

.7826E+01 

.1171E+02 

.2821E-01 

.5984E-02 
-. 1079E+00 
. 3683E-01 

-.4229E+01 

-.1293E+02 

.3290E-02 

Plant 

with  Actuator  Dynamics  and  Computational 

Time  Delay 

Alpha 

Epsilon 

Sigma** 

If  *** 

-0 

.  4545 

2.2 

1.9 

1.5 

0.025 

.7435E+01 
.1113E+02 
. 2679E-01 

.8976E-02 

-.1618E+00 

.5524E-01 

-.7049E+01 
-.  215.5E+02 
. 5484E-02 

Notes 

^Each  pulse  entry  in  v  has  four  parts: 

1.  The  time  in  seconds  the  input  reaches  steady- state; 

2.  Steady-state  value  (angles  in  radians); 

3.  Time  the  input  leaves  steady- state;  and 

4.  Time  the  input  reaches  zero. 

**Diagonal  element  values  of  the  matrix. 

***Kq  *  alpha  *  K^. 


E-4 


TABLE  E-4 


DESIGN  PARAMETERS  AND  CONTROLLER  MATRICES 
FOR  LONGITUDINAL  CONTROLLERS 


PLIGHT  CONDITIONi  1.2  Mach  at  15000  Feet 
MANEUVER:  Vertical  Translation 
COMMAND  VECTOR  v*s  0,  0,  0,  0 

0,  0,  0;  0 

0,  1.1,  -0.0218,  50,  50 


Plant  Only 


Alpha 

Epsilon 

Sigma** 

V  *** 

-0 

1.8 

.1975E+01 

.1823E-02 

-.8799E+01 

0.4545 

1.8 

1.0 

-.8702E+01 

-.5346E-01 

-.  284  8E+02 

0.05 

.8307E-02 

.9291E-02 

-.1467E-01 

Plant  and  Computational  Time  Delay 


Alpha 

Epsilon 

Sigma** 

X  WWW 

-0 

1.8 

.1975E+01 

•  .1823E-02 

-.1760E+02 

0.4545 

1.8 

1.0 

.8702E+01 

-.5346E-01 

-. 5695E+02 

0.1 

.8307E-02 

.9291E-02 

-.2935E-01 

Plant  with 

Actuator  1 

oynamics 

Alpha 

Epsilon 

Sigma** 

K  *** 

-0 

1,8 

.1975E+01 

.1823E-02 

-.1056E+02 

0.4545 

1.8 

1.0 

.8702E+01 

-.5346E-01 

-. 3417E+02 

0.06 

.8307E-02 

. 9291E-02 

-.1761E-01 

Plant 

with  Actuator  Dynamics  and  Computational 

Time  Delay 

Alpha 

Epsilon 

Sigma** 

K  *** 

-0 

1.8 

.1975E+01 

.1823E-02 

-.1232E+02 

0..4545 

1.8 

1.0 

'.8702E+01 

-.  5346E-01 

-.3987E+02 

0.07 

.8307E-02 

.9291E-02 

-.2054E-01 

Notes 

^Each  pulse  entry  in  v  has  four  parts; 

1.  The  time  in  seconds  the  input  reaches  steady-state; 

2.  Steady-state  value  (angles  in  radians); 

3.  Time  the  input  leaves  steady- state;  and 

4.  Time  the  input  reaches  zero. 

’''*Diagonal  element  values  of  the  matrix. 

***K  »  alpha  *  K^, 


E-5 


TABLE  E-5 


DESIGN  PARAMETERS  AND  CONTROLLER  MATRICES 
FOR  LONGITUDINAL  CONTROLLERS 


PLIGHT 

'  CONDITION;  0.7  Mach  at  15000  Feet 

MANEUVER;  Pitch  Pointing 

COMMAND  VECTOR 

V*:  0^  0.8;  0.0349;  1.0 

,  1.8 

0,  0;  0;  0 

0,  0.8;  0.0349;  50; 

50 

Plant  Only 

Alpha 

Epsilon 

Sigma** 

Kq*** 

1.9 

.7432E+01  .5982E-02 

-.5637E+01 

0.357 

2.8 

1.0 

-.1112E+02  -.1079E+00 

-.1724E+02 

0.02 

.2678E-01  .3681E-01 

.4385E-02 

Plant  and  Computational  Time  Delay 

Alpha 

Epsilon 

Sigma** 

Kq*** 

1.9 

.7432E+01  .5982E-02 

-  .7046E+01 

0.357 

2.8' 

1.0 

-.1112E+02  -.1079E+00 

-.2154E+02 

0.025 

.2678E-01  .3681E-01 

.5482E-02 

Plant  with  Actuator  Dynamics 

Alpha 

Epsilon 

Sigma** 

Kq*** 

1.9 

.7432E+01  .5982E-02 

-.4227E+01 

0.357 

2.8 

1.0 

-.1112E+02  -.1079E+00 

-.1293E+02 

0.015 

.2678E-01  .3681E-01 

.3289E-02 

Plant 

with  Actuator  Dynamics  and  Computational  Time  Delay 

Alpha 

Epsilon 

Sigma** 

K  ■**  * 

-0 

1.7 

1 

.6649E+01  .5982E-02 

-.4227E+01 

0.357 

2.8 

1.0 

-.9951E+01  -.1079E+00 

-.1293E+02 

0.015 

.2396E-01  .3681E-01 

. 3289E-02 

Notes 

"^Each  pulse  entry  in  v  has  four  parts; 

1.  The  time  in  seconds  the  input  reaches  steady- state; 

2.  Steady-state  value  (angles  in  radians); 

3.  Time  the  input  leaves  steady- state;  and 

4.  Time  the  input  reaches  zero. 

**Diagonal  element  values  of  the  matrix. 

***JSo  =  alpha  *  Kj^. 


TABLE  E-6 


DESIGN  PARAMETERS  AND  CONTROLLER  MATRICES 
FOR  LONGITUDINAL  CONTROLLERS 


FLIGHT  CONDITION;  1.2  Mach  at  15000  Feet 

MANEUVER:  Pitch  Pointing 

COMMAND  VECTOR  V*!  0.8,  0.0218,  1.0,  1.8 

0,  0,  0,  0 

0.8,  0.0218,  50,  50 


Plant  Onli 


Alpha 
0.  357 

Alpha 


0.  357 


Epsilon  Sigma** 
1.9 

2.8  0.2 

0.07 


.2084E‘t-01  .3645E-03  -.1231E-I-02 

-.9182E+01  -.1069E-01  -.3985E+02 

.8765E-02  .1857E-02  -.2054E~01 


Plant  and  Computational  Time  Delay 


K 


Epsilon  Sigma** 


1.9 

2.8  0.2 

0.16 


.2084E+01  .3645B-03  -.2815E-I-02 

-.9182E+01  -.1069E-01  -.9109E+02 

.8765E-02  .1857E-02  -.4694E-01 


Plant  with  Actuator  Dynamics 

V  )ilr  A 

Alpha  Epsilon  Sigma**  _ ^0 _ 


1.9 

0.357  2.8  0.2  -.9182E+01  -.1069E-01  -.2846E^ 

0.05  .8765E-02  .1857E-02  -.1467E- 

Plant  with  Actuator  Dynamics  and  Computational  Time  Dela 


.2084E+01  .3645E-03  -.8796E+01 

-.9182E+01  -.1069E-01  -.2846E+02 

.8765E-02  .1857E-02  -.1467E-01 


A1 

pha 

Epsilon 

Sigma** 

-0 

1.  9 

. 2084E+01 

.3645E-03 

-.1231E+02 

0. 

357 

2.8 

0.2 

-.9182E+01 

-.1069E-01 

-. 3985E+02 

0.07 

.8765E-02 

.1857E-02 

-. 2054E-01 

Notes 

^Each  pulse  entry  in  v  has  four  parts; 

1.  The  time  in  seconds  the  input  reaches  steady-state 

2.  Steady-state  value  (angles  in  radians); 

3.  Time  the  input  leaves  steady- state;  and 

4.  Time  the  input  reaches  zero. 

^*Diagonal  element  values  of  the  matrix. 

***K„  =  alpha  *  K.  . 


TABLE  E-7  — Continued 


Plant  with  Actuator  Dynamics 


Output 

Peak  Value** 

Peak  Time 

Settling  Time 

Pitch  Rate 

2.6 

1.575 

*** 

Forward  Velocity 

“0.02 

2.8 

*** 

Angle  of  Attack 

1.34 

1.925 

*** 

P’^-Th  Angle 

7.7  . 

4.0 

7.8 

F.  jht  Path  Angle 

8.4 

5.2 

8.0 

Plant  with  Actuator  Dynamics  and  Computational  Time  Delay 


Oi’tgut 

Peak  Value** 

Peak  Time 

!  Settling  Time 

Pitch  Rate 

2.7 

1.225 

**  * 

F  rward  Velocity 

-0.038 

2.45 

*** 

Angle  of  Attack 

1.45 

1.925 

*  *  A 

P3-^ch  Angle 

7.7 

3.9 

4.0 

Fj.  .ght  Path  Angle 

8.4 

5.2 

8.2 

No  as 

1.  Figures  of  merit  for 

the  pitch 

angle  and  the 

flight  path  angle 

are  estimated 

from  the  response  plots. 

2.  Pitch 

rate ,  forward 

velocity,  and  angle  of 

attack  are  components  of  the  output  vector 

Z- 

3.  Pitch 

rate  and  angle 

of  attack 

are  commanded 

to  achieve  desired  responses  in  the  pitch  angle  and  the 
flight  path  angle  for  each  maneuver. 


*Coimnanded  inputs  in  order: 

(1)  pitch  rate  in  radians  per  second; 

(2)  forward  velocity  in  feet  per  second :  and 

(3)  angle  of  attack  in  radians. 

**Angles  are  converted  from  radians  to  d.egcees  nrior 
to  calculating  the  figures  of  merit. 

***Final  value  of  the  output  is  zero.  Settling  time 
undefined  using  definition  of  2  percent  of  final  value. 


TABLE  E-8 


LONGITUDINAL  FIGURES  OF  MERIT 

FLIGHT  CONDITION:  1*2  Mach  at  15000  Feet 

MANEUVER:  Direct  Clinto 

COMMAND  VECTOR  v:*  !•!/  0.0436/  3,  4.1 

0/  0,  0/  0 

0/  0/  0/  0 


Plant  Only 


Output 

Pitch  Rate 
Forward  Velocity 
Angle  of  Attack 
Pitch  Angle 
Flight  Path  Angle 


Plant 


Output 

Pitch  Rate 
Forward  Velocity 
Angle  of  Attack 
Pitch  Angle 
Flight  Path  Angle 


Peak  Value** 

Peak  Time 

Sett line  Time 

2.59 

1.4 

*** 

-0.04 

2.275 

*** 

0.756 

1.75 

*** 

7.6 

3.9 

3.9 

7.9 

4.4 

4.6 

and  Computational  Time  Delay 


Peak  Value** 

Peak  Time 

Settling  Time 

2.66 

1.225 

*  *  * 

-0.08 

2.1 

*  *  * 

0.78 

1.75 

*** 

7.7 

3.9 

4.0 

7.9 

4.2 

4.5 

E-10 


TABLE  E-8  — Continued 


Plant  with  Actuator  Dynamics 

Output 

Peak  Value** 

Peak  Time 

Settling  Time 

Pitch  Rate 

2.64 

1.4 

*** 

Forward  Velocity 

-0.039 

1.925 

Angle  of  Attack 

0.77 

1.575 

*** 

Pitch  Angle 

7.6 

3.9 

3.9 

Flight  Path  Angle 

7.8 

4.2 

5.5 

Plant  with  Actuator  Dynamlos  and  Computational  Time  Delay 


Output 

Peak  Value** 

Peak  Time 

Settling  Time 

Pitch  Rate 

2.79 

1.225 

iekh 

Forward  Velocity 

-0.069 

1.4 

*** 

Angle  of  Attack 

0.825 

1.575 

itifk 

Pitch  Angle 

7.7 

3.9 

4.0 

Flight  Path  Angle 

7.8 

4.2 

5.5 

Notes 


1.  Figures  of  merit  for  the  pitch  angle  and  the 
flight  path  angle  are  estimated  from  the  response  plots. 

2.  Pitch  rate,  forward  velocity,  and  angle  of 
attack  are  components  of  the  output  vector 

3.  Pitch  rate  and  angle  of  attack  are  commanded 
to  achieve  desired  responses  in  the  pitch  angle  and  the 
flight  path  angle  for  each  maneuver. 

*Commanded  inputs  in  order; 

(1)  pitch  rate  in  radians  per  second; 

(2)  forward  velocity  in  feet  per  second;  and 

(3)  angle  of  attack  in  radians. 

**Angles  are  converted  from  radians  to  degrees  prior 
to  calculating  the  figures  of  merit. 

***Final  value  of  the  output  is  zero.  Settling  time 
undefined  using  definition  of  2  percent  of  final  value. 


E-11 


TABLiS  E-9 


LONGITUDINAL  FIGURES  OF  MERIT 

FLIGHT  CONDITION:  0*7  Mach  at  15000  Feet 

MANEUVER;  Vertical  Translation 

COMMAND  VECTOR  v:*  0,  0,  0,  0 

“  0,  0,  0,  0 

1.1,  -0.0349,  50,  50 


Plant  Only 


Output 

Peak  Value** 

Peak  Tine 

Settling  Time 

Pitch  Rate 

-0.11 

1.575  • 

*** 

Forward  Velocity 

0.009 

1.925 

*** 

Angle  of  Attack 

-2.03 

4.025 

4 . 55 

Pitch  Angle 

-0.25 

3.6 

3.7 

Flight  Path  Angle 

1.8 

3 . 6 

3.7 

Plant 

and  Computational 

Time  Delay 

Output 

Peak  Value**  Peak  Time  Settling  Time 

Pitch  Rate 

-0.229 

1.575  *** 

Forward  Velocity 

0.003 

1.925  *** 

Angle  of  Attack 

-2.01 

3.85  2.975 

Pitch  Angle 

-0.4 

3.6  3.7 

Flight  Path  Angle 

1.7 

3.7  3.8 

E-12 


TABLE  E-9  — Continued 


Plant  with  Actuator  Dynamics 

Output 

Peak  Value** 

Peak  Time 

Settling  Time 

Pitch  Rate 

-0.18 

1.575 

*** 

Forward  Velocity 

.002 

1.75 

*** 

Angle  of  Attack 

-2,08 

3.5 

4.55 

Pitch  Angle 

-0.2 

3.5 

3.5 

Flight  Path  Angle 

1.9 

3.5 

3.5 

Plant  with  Actuator  Dynamics  and  Computational  Time  Delay 


Output 

Peak  Value** 

Peak  Time 

Settling  Time 

Pitch  Rate 

-0.394 

1.  575 

•k  liit 

Forward  Velocity 

0.004 

1.925 

*** 

Angle  of  Attack 

-2.007 

3.675 

2.975 

Pitch  Angle 

-0.4 

4.0 

4.1 

Flight  Path  Angle 

1.6 

3.6 

3.7 

Notes 


1.  Figures  of  merit  for  the  pitch  angle  and  the 
flight  path  angle  are  estimated  from  the  response  plots. 

2.  Pitch  rate/  forward  velocity/  and  angle  of 
attack  are  components  of  the  output  vector  y. 

3.  Pitch  rate  and  angle  of  attack  are  commanded 
to  achieve  desired  responses  in  the  pitoh  angle  and  the 
flight  path  angle  for  each  maneuver. 

^Commanded  inputs  in  order: 

(1)  pitch  rate  in  radians  ^  i;  second? 

(2)  forward  velocity  in  fee.  ner  second;  ar.d 

(3)  angle  of  attack  in  radian. 

**Anglea  are  converted  from  radians  to  .  prior 

to  calculating  the  figures  of  merit. 

**'*Final  value  of  the  output  is  zero.  Settling  time 
undefined  using  definition  of  2  percent  of  final  value. 


E-13 


TABLE  E-10 


LONGITUDINAL  FIGURES  OF  MERIT 

FLIGHT  CONDITION;  1*2  Mach  at  15000  Feet 

MANEUVER;  Vertical  Translation 

COMMAND  VECTOR  v;*  0,  0,  0,  0 

■“  0,  0,  0,  0 

1.1,  -0.0218,  50,  50 


Output 

Pitch  Rate 
Forward  Velocity 
Angle  of  Attack 
Pitch  Angle 
Flight  Path  Angle 


Output 

Pitch  Rate 
Forward  Velocity 
Angle  of  Attack 
Pitch  Angle 
i''  '.ight  Path  Angle 


Peak  Value** 

Peak  Time 

Settlincr  Time 

-0.09 

1.225 

*  *  lit 

0.004 

1.225 

**4t 

-1.25 

15.9 

5.575 

0.16 

3.0 

3.0 

1,19 

3.2 

3.2 

and  Computational  Time  Delay  ' 

Peak  Value** 

Peak  Time 

Settlinq  Time 

-0.19 

1.225 

*  ** 

0.10 

1.225 

*  ♦  * 

-1.25 

15.9 

5.9 

-0.30 

3,1 

3.1 

.95 

3.1 

3.1 

TABLE  E-10 — Continued 


Plant  with  Actuator  Dynamics 

Output 

Peak  Value** 

Peak  Time 

Settling  Time 

Pitch  Rate 

0.424 

0.175 

*  *  * 

Forward  Velocity 

0.009 

1.4 

Angle  of  Attack 

-1.25 

15.9 

6.1 

Pitch  Angle 

-0.2 

2.1 

2.3 

Flight  Path  Angle 

1.9 

1.6 

1.8 

Plant  with  Actuator  Dynamics  and  Computational  Time  Delay 


Output 

Peak  Value** 

Peak  Time 

Settling  Time 

Pitch  Rate 

-0.68 

1.4 

*  ** 

Forward  velocity 

0.18 

1.4 

*** 

Angle  of  Attack 

1.25 

15.7 

6.0 

Pitch  Angle 

-0.3 

2.0 

2.1 

Flight  Path  Angle 

1.0 

1.75 

3.0 

Notes 

1.  Figures  of  merit  for  the  pitch  angle  and  the 
flight  path  angle  are  estimated  from  the  response  plots. 


2,  Pitch  rate,  forward  velocity,  and  angle  of 
attack  are  components  of  the  output  vector 

3.  Pitch  rate  and  angle  of  attack  are  commanded 
to  achieve  desired  responses  in  the  pitch  angle  and  the 
flight  path  angle  for  each  maneuver. 

*Commanded  inputs  in  order: 

(1)  pitch  rate  in  radians  per  second; 

(2)  forward  velocity  in  feet  per  second;  and 

(3)  angle  of  attack  in  radians. 

** Angles  are  converted  from  radians  to  degrees  prior 
to  calculating  the  figures  of  merit. 

***Final  value  of  the  output  is  zero.  Settling  time 
undefined  using  definition  of  2  percent  of  final  value. 
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TABLE  E-11 


LONGITUDINAL  FIGURES  OF  MERIT 


PLIGHT  CONDITION:  0-7  Mach  at  15000  Feet 

MANEUVER:  Pitch  Pointing 

COMMAND  VECTOR  v:*  0,  0.8,  0.0349,  1.0,  1.8 

0 ,  0  ,  0 ,  0 

0,  0.8,  0.0349,  50,  51 


Plant  Only 


Output 

Peak  Value** 

Peak  Time 

Settling  Time 

Pitch  Rate 

2.16 

0.875 

•khii 

Forward  Velocity 

-0.01 

1.225 

•kitii 

Angle  of  Attack 

2.025 

1.4 

4.025 

Pitch  Angle 

2.2 

1.9 

2.0 

Flight  Path  Angle 

0.5 

2.3 

4.0 

Plant 

and  Computational  Time 

Delay 

Output 

Peak  Value** 

Peak  Time  Settling  Time 

Pitch  Rate 

2.36 

1.05 

*** 

Forward  Velocity 

-0.02 

1.225 

*** 

Angle  of  Attack 

2.01 

5.775 

4.6 

Pitch  Angle 

2.3 

1.95 

2.0 

Plight  Path  Angle 

0.7 

2.5 

5.0 

TABLE  E-11 — Continued 


Plant  with  Actuator  Dynamics 

Output 

Peak  Value"* 

Peak  Time 

Settling  Time 

Pitch  Rate 

2.28 

1.05 

4t  *  4r 

Forward  Velocity 

-0.01 

1.05 

*** 

Angle  of  Attack 

2.03 

5.075 

4.75 

Pitch  Angle 

2.1 

1 . 95 

2.5 

Flight  Path  Angle 

0.4 

2.3 

4.3 

Plant  with  Actuator  Dynamics  and 

Computational  Time  Delay 

Output 

Peak  Value"* 

Peak  Time 

Settling  Time 

Pitch  Rate 

2.42 

1.05 

*** 

Forward  velocity 

-0.023 

1.05 

h-kh 

Angle  of  Attack 

2.002 

7.7 

5.6 

Pitch  Angle 

2.0 

1.95 

5.5 

Flight  Path  Angle 

0.9 

2.3 

5.6 

Notes 

1.  Figures  of  merit  for  the  pitch  angle  and  the 
flight  path  angle  are  estimated  from  the  response  plots. 

2.  Pitch  rate,  forward  velocity,  and  angle  of 
attack  are  components  of  the  output  vector 

3.  Pitch  rate  and  angle  of  attack  are  commanded 
to  achieve  desired  responses  in  the  pitch  angle  and  the 
flight  path  angle  for  each  maneuver. 

''’Commanded  inputs  in  order: 

(1)  pitch  rate  in  radians  per  second; 

(2)  forward  velocity  in  feet  per  second;  and 

(3)  angle  of  attack  in  radians. 

’"’''Angles  are  converted  from  radians  to  degrees  prior 
to  calculating  the  figures  of  merit. 

""’''Final  value  of  the  output  is  zero.  Settling  time 
undefined  using  definition  of  2  percent  of  final  value. 


TABLE  E-12 


LONGITUDINAL  FIGURES  OF  MERIT 


FLIGHT  CONDITIONS  1.2  Mach  at  15000  Feet 


MANEUVERS  Pitch  Pointing 


COMMAND 

VECTOR  vs*  0.8,  0.0218, 
0,  0,  0,  0 
0.8,  0.0218, 

1.0,  1. 

50,  50 

8 

Output 

Plant  Only  ' 

Peak  Value**  Peak  Time 

Settling  Time 

Pitch  Rate 

1.36 

0.875 

*** 

Forward  Velocity 

-0.05 

1.05 

*** 

Angle  of  Attack 

1.25 

15.9 

2.975 

Pitch  Angle 

1.35 

1.  9 

1.9 

Flight  Path  Angle 

-0.40 

1.0 

2.5 

Plant 

and  Computational  Time  Delay 

Output 

Peak  Value** 

Peak  Time 

Settling  Time 

Pitch  Rate 

1.47 

1.05 

**  * 

Forward  Velocity 

-0.12 

1.225 

•kit  it 

Angle  of  Attack 

1.253 

1.225 

2.8 

Pitch  Angle 

1.4 

1.9 

2.0 

Flight  Path  Angle 

-0.4 

1.0 

2.5 

TABLE  E-12 — Continued 


Output 


Plant  with  Actuator  Dynamics 

Peak  Value**  Peak  Time  Settlino  Time 


Pitch  Rate 
Forward  Velocity 
Angle  of  Attack 
Pitch  Angle 
Flight  Path  Angle 


1.7 

1.05 

it  -kit 

-0.085 

1.05 

it  Mil 

1.28 

1.225 

3.15 

1.25 

1.9 

3.0 

-0.4 

1.0 

2.5 

Plant  with  Actuator  Dynamics  and  Computational  Time  Dela 
Output  Peak  Value**  Peak  Time  Settlina  Time 


Pitch  Rate 
Forward  Velocity 
Angle  of  Attack 
Pitch  Angle 
Flight  Path  Angle 


1.0 
1.225 
1.2 
2.5 
1.9 


1.  Figures  of  merit  for  the  pitch  angle  and  the 
flight  path  angle  are  estimated  from  the  response  plots. 

2.  Pitch  rate,  forward  velocity,  and  angle  of 
athack  are  components  of  the  output  vector  y. 

3.  Pitch  rate  and  angle  of  attack  are  commanded 
to  achieve  desired  responses  in  the  pitch  angle  and  the 
flijht  path  angle  for  each  maneuver. 

*Commanded  inputs  in  order; 

(1)  pitch  rate  in  radians  per  second; 

(2)  forward  velocity  in  feet  per  second;  and 

(1)  angle  of  attack  in  radians. 

**Ancjl«s  arn  converted  from  radians  to  degrees  prior 
to  calculatir  j  i',;'’  figures  of  merit. 

***Final  value  of  the  output  is  zero.  Settling  time 
undefincJ  using  definition  of  2  percent  of  final  value. 
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Appendix  F:  Longitudinal  Controller  Respon se  Plots 

This  appendix  includes  longitudinal  response  plots 
for  the  X-29  closed- loop  flight  control  system  incorporating 
a  controller  designed  with  the  Porter  method.  Response  plots 
are  provided  for  the  following  maneuvers  and  flight  condi¬ 
tions  : 


Direct  Climb 


(1) 

0.4 

Mach 

at 

Sea  Level 

(2) 

1.2 

Mach 

at 

15000  Feet 

Vertical  Translation 

(1) 

0.7 

Mach 

at 

15000  Feet 

(2) 

1.2 

Mach 

at 

15000  Feet 

Pitch  Pointing 

(1) 

0.7 

Mach 

at 

15000  Feet 

(2)  1.2  Mach  at  15000  Feet 
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Appendix  G:  Lateral  State  Space  Matrices 


Tables  in  this  appendix  list  the  linearized 
lateral  state  space  matrices  of  the  X-29  aircraft  for  the 
following  flight  conditions: 

1.  0.4  Mach  at  Sea  Level 

2 .  0.9  Mach  at  Sea  Level 

3.  1.2  Mach  at  15,000  Feet 

4.  0.9  Mach  at  50,000  Feet 


TABLE  G-1 


LATERAL  STATE  SPACE  MATRICES 


FLIGHT  CONDITION: 

0.4  Mach  at  Sea 

Level 

A  (Plant  Matrix) 

.OOOOE-I-OO 

.7181E-01 

.OOOOE+OO 

.OOOOEtOO 

.OOOOE+OO 

-.2052E+00 

-.9003E+01 

.5079E+01 

.lOOOE+01 

.7971E-01 

-.2637E+01 

“.8816E-01 

.7996E-01 

-.9968E+00 

.1124E+01 

-.1061E+00 

B  (Control 

Input  Matrix) 

.OOOOE+OO 

.7953E-03 

.9747E+00 

.6165E-01 

.OOOOE+OO 

.8768E-03 

.1503E+00 

-.4984E-01 

, 

C  (Output  Matrix) 

.lOOOE+01 

•OOOOE+OO 

.OOOOE+OO 

.lOOOE+01 

.OOOOE+OO 

.OOOOE+OO 

.OOOOE+OO 

.OOOOE+OO 

Notes 

1. 

slip  angle ^ 

States  (listed 
roll  rate ,  and 

in  order)  are  bank  angle,  side- 
yaw  rate. 

2. 

Control  inputs 

(listed  in  order) 

are  flaperon 

and  rudder. 


3.  Outputs  (listed  In  order)  are  bank  angle  and 
sideslip  angle. 


TABLE  G-2 


LATERAL  STATE  SPACE  MATRICES 


PLIGHT 

CONDITION: 

0 . 9  Mach  at  Sea 

Level 

A  (Plant  Matrix) 

.OOOOE+00 

.3119E-01 

.OOOOE+00 

.OOOOE+00 

.OOOOE+OO 

-.4664E-I-00 

-.6582E-»-02 

.1864E+02 

.lOOOE+01 

.4810E-01 

-.8976E+01 

-.2631E+00 

.4385E-01 

-.9990E+00 

.2229E-)-01 

-.1606E+00 

B  (Control 

Input  Matrix) 

.OOOOE-fOO 
-.2465E-02 
■ i2784E+01 
.1805E+00 

.OOOOE-fOO 

.1395E-02 

.60O7E-f0O 

-.1832E+00 

C  (Output  Matrix) 

.lOOOE+01 

.OOOOE+00 

.OOOOE+00 

.lOOOE+01 

.OOOOE+00 

.OOOOE+00 

.OOOOE+00 

.OOOOE-fOO 

Notes 


1.  States  (listed  In  order)  are  bank  angle,  side¬ 
slip  angle,  roll  rate,  and  yaw  rate. 

2.  Control  inputs  (listed  in  order)  are  flaperon 
and  rudder. 

3.  Outputs  (listed  in  order)  are  bank  angle  and 
sideslip  angle. 


Appendix  H:  Lateral  Controller  Response  Plots 

This  appendix  Includes  lateral  time  history'  plots 
for  the  X-29  closed-loop  flight  control  system  Incorporating 
a  controller  designed  with  the  Porter  method.  The  response 
plots  are  presented  for  the  following  maneuvers  and  flight 
conditions: 

Coordinaced  Turn 

1.  0.4  Mach  at  Sea  Level 

2.  0.9  Mach  at  Sea  Level 

3.  1.2  Mach  at  15,000  Feet 

4.  0.9  Mach  at  50,000  Feet 

Beta  Pointing 

1.  0.4  Mach  at  Sea  Level 

2.  0.9  Mach  at  Sea  Level 

3.  1.2  Mach  at  15,000  Feet 

4.  0.9  Mach  at  50,000  Feet, 

The  lateral  design  parameters,  controller  matrices, 
and  figures  of  merit  for  these  plots  are  located  in  Tables 
6-2  through  6-5  of  Chapter  6. 


33  03  8  03  »  33 

SllONB  <il  IS-lOiS  ONB  VNBS 


COnriRNOtO  BLTfl  RMCIL  I  OLD 

^.CO  2  CO  4. CO  6.CC  8. CO 


ROI  I  (IMD  raw  RqTL5 


I  03'*-  00  a-  o: 

H»A  ONb  1)021 


00  Oi  D0'7l 

(oia  j 


Appendix  I;  Other  Design  Approaches 


As  discussed  in  Chapter  7,  Conclusions  and  Recom¬ 
mendations,  an  unsuccessful  attempt  was  made  to  control 
acceleration  first  by  implementing  the  actuator  dynamics 
into  the  state  equations,  thus  eliminating  the  necessity 
for  a  feedforward  D  matrix.  The  result  was  transmission 
zeros  in  the  right  half  s-plane.  For  the  benefit  of  future 
work  in  this  area  this  appendix  expands  on  the  results  by 
providing  a  sample  of  the  actual  transmission  zeros  calcu¬ 
lated  on  TOTAL  using  the  Macro  program  written  by  Barfield 
(1).  The  same  results  were  confirmed  using  the  ZERO  pro¬ 
gram  written  by  Lewis  >1)  .  For  the  flight  condition  of 
0.7  Mach  at  15,000  feet,  the  following  longitudinal  equa¬ 
tions  of  motion  and  resulting  transmission  zoros  are  listed 
as  follows: 


I-l 


-.0002986  18.92  -.4928  .1720  -.05630  -.1929 

0  00'  -20.0  0  0 

0  0  0  0  -5.0  0 


Transmission  Zero  Locations 


-.36E-05 

-.16E-03 

-290.9 

+15.16 


For  the  lateral  equations  of  motion  at  0.7  Mach  at  15/000 
feet/  the  equations  and  resulting  transmission  zeros  are 
listed  below: 


• 

p 

-2.989 

-17.93 

0 

• 

& 

.0571 

-.2210 

.04334 

• 

<t) 

1 

0 

0 

■ 

.  . 

k 

-.09651 

7.836 

0 

• 

ft 

A 

A 

1.090  1.493  .2670 

-.9984  -.0008361  .0009071 

.05710  0  0 

* 

-.09504  .09608  -.08655 

R 

0  -10.0  0 

w 


Transmission  Zero  Locations 
-0.002476 

-2.048  .25  0  0 

for  M  * 

-4.0  000 

+1.952 


-0.002475 
-1.0 
-1.065 
+0.9385 

The  transmission  zeros  lie  even  farther  in  the  right  half 
s-plane  if  the  equations  are  changed  so  the  output  ^ 
vector  is  bank  angle  4)  and  sideslip  angle  3*  For  other 
measurement  matrix  values  and  locations ^  the  results  are 
similarly  unsatisfactory. 
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is  used  to  develop  and  refine  the  control  laws.  MULTI  also 
simulates  the  complete  closed-loop  control  system  and 
generates  appropriate  time  response  plots  for  analysis,)? 


'Aircraft  dynamics  for  several  points  in  the  flight 
envelope  are  represented  by  linearized  state  space  equations 
obtained  from  NASA  Dryden,  an  agency  responsible  for  the 
development  and  testing  of  the  X-29.  Decoupled  longi¬ 
tudinal  and  lateral^equations  are  used  to  design  separate 
longitudinal  and  lateral  controllers. 


Control  lav?s  are  developed  to  stabilize  the  aircraft 
and  perform  longitudinal  maneuvers  (direct  climb,  vertical 
translation,  and  beta  pointing)  at  three  different  fligiit 
conditions,  with  and  without  first  order-actuator  dynamics 
and  computational  time  delay  added  to  the  simulation.  The 
responses  are  com, pared. Although  some  degradation  of  per¬ 
formance  is  observed  wher.  the  actuator  dynamics  and  time 
delay  are  added,  the  results  show  fast  and  well-behaved 
control  of  the  aircraft /motion  parameters  with  quick  and 
sm.ooth  tracking  of  pilot  input  commands.  Similarly, 
lateral  control  lav?s ^re  developed  for  the  coordinated 
turn  and  beta  pointing  maneuvers  at  five  different  flight 
conditions  with  c^rfiparable  results. 
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